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SECTION 1. INTRODUCTION 


After the National Aeronautics and Space Administration (NASA) 
Office of Space Science (OSS) announcement of opportunities for partici- 
pation in the Apollo Telescope Mount (ATM) solar physics mission on Sky- 
lab in January 1971 and the receipt of resultant proposals, an ad hoc 
committee of solar scientists from throughout the United States selected 
a number of efforts to be pursued. These programs involved some X ray- 
instrumented rocket flights as well as the ground-based observational 
efforts that the Marshall Space Flight Center (MSFC) was asked to manage. 
The ropket flight tasks were managed directly from NASA headquarters. 

Because of its scientific nature, the Skylab Ground-Based Astronomy 
Program (SGAP) was assigned to MSFC’s Space Sciences Laboratory, where the 
author was project manager, Dr, Mona Hagyard was project scientist, and 
individual task scientists (see Table 1.1) were responsible for the se- 
lected efforts as contracting officer’s representatives (COR). 

Nine tasks were funded, with the participating astronomical obser- 
vatories and their associated projects as follows: 

A. The University of Hawaii constructed a photoelectric differen- 
tial coronal photometer for an observational investigation of coronal 
active regions. This Instrument measures simultaneously the intensities 
of several visible coronal lines suitably chosen for their diagnostic 
properties. Physical properties are studied to determine from this 
primary data the rates of energy loss and gain from the active regions 
and the effects of flare events on the corona. The Principal Investigator 
(PI) was Dr. John Jefferies. 

B. Lockheed Missiles and Space Company (LMSC) operated, at Kitt 
Peak Observatory’s McMath solar telescope, a system called spectra- 
spectrohellography that is capable of mapping the physical parameters of 
the solar atmosphere. In addition to the Kitt Peak solar telescopes and 
vacuum spectrograph, a wide-exit aperture and a specially constructed 
movie camera capable of rapid film advance were used to obtain spectral 
maps of regions of the solar disk with high spatial resolution (to one- 
half arc second) . A high-speed microdensitometer-computer system was 
added to allow rapid analysis of the spectral data to obtain solar velo- 
city and magnetic field maps of regions of interest. Dr. Alan Title was 
the Principal Investigator. 

C. The National Bureau of Standards (NBS) was to upgrade calibra- 
tion capabilities in support of ATM-related measurements as a part of the 
SGAP program in the following areas: 

1. A hydrogen arc source of known radiant flux for the cali- 
bration of spectrometric-detector systems over the region of 500 to 3700 
X was developed. 



INSTITUTION 

INVESTIGATOR 

DESCRIPTION 

University of Hawaii 

Jefferies 

Construction & operation during ATM of photo- 
electric differential coronal photometer in 
various simultaneous visible lines. 

Lockheed Missiles and 
Space Company 

Title 

Addition of microdensitometer/ computer system 
to S^Hg & modeling of solar data. 

National Bureau of Standards 

Wiese - 
Madden 

Develop H2 arc UV calib. source (500-3700A) . 

Study contamination effects on photocathode 
response In VUV (500-1500) . 



Develop sjrnchrotron calib, capability (600 - 
200 A). 

California Institute 
of Technology 

Zirin 

Construct equatorial mount and vacuum system, 

& Install 65 cm FVU photoheliograph at Big Bear. 



Purchase and install at Big Bear, filter & IR 
vidicon for 10830 A He & 3933 A Ca-K line 
f iltergrams. 

Lockheed Missiles and 
Space Company 

Mart in- 
Ramsey 

Observations of D3 line of Hel (5876 A) 
(photographic) , 



Filter studies In Call (8542 A). 

University of California 
at San Diego 

Hudson 

Construct beam switching radiometer for the 
60” Mt. Lemon telescope for meas, at 1 mm - 
1 u & observations during ATM, 

Lockheed Missiles and 
Space Company 

Kulander 

Theoretical study of He emission in visible 
(Hel-p3(?5876 A) and UV (Hel Q 584 & 537 A) . 

Uttar Pradesh State 
Observation (India) 

Pande 

Equipment procurement & loan to allow 
observations for molecular line studies. 

Johns Hopkins 
University 

Gotwols 

Observation & analysis of 565-1000 MHz spectra 
from APL 60* dish. 


CONTRACTING OFFICER 
REPRESENTATIVE 

J. McGuire 

B. Duncan 
R. Linton 

E« Reichmann 

E* Miller 

W, Fountain 

Dr« Hagyard 
H. Miller 

C . Baugher 
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2. A study was conducted to determine the effects on photo- 
cathodes caused by the removal or addition of monolayers of contaminants 
in vacuum; the wavelength region of interest was 500 to 1500 A, 

3. A capability was established for radiometric calibration 
by utilization of the NBS synchrotron facility. Windowless diodes were 
developed as transfer standards in this spectral range. Drs. W. Wiese 
and R. Madden were the Principal Investigators. 

D. The California Institute of Technology installed the 65-cm 
FVU (functional verification unit) photoheliograph at its observatory at 
Big Bear Lake in California. The installation of the new vacuum photo- 
heliograph provides Big Bear with a large-aperture telescope capable of 
fully exploiting the site’s excellent observing conditions. The new unit 
was used in conjunction with ATM for extremely high resolution studies of 
active regions . Provisions were made for use of various cameras and ^ 
detector systems as well as feeding a vacuum spectrograph at the Coude 
focus. Filtergrams were made in lines extending from the 3933 A calcium 
K-line to the 10830 A helium line. Dr. Harold Zirin was the Principal 
Investigator. 

E. The Lockheed Solar Observatory at Rye Canyon prepared both of 
its spar telescope systems for observations during the ATM mission. On 
one spar, studies in the line of Hel were to direct attention to obser- 
ving solar flares and transient events during periods of high disk activity 
and to prominence observations during periods of low disk activity. On the 
other spar, a telescope was fitted with a filter for high resolution photo- 
graphic studies in Call at 8542 A. This line is believed to originate at a 
higher level in the solar atmosphere than the H- a line and should be 
extremely valuable in relating the ATM X ray and XUV data to filtergrams 
and spectroheliograms taken at wavelengths originating at lower levels in 
the chromosphere. Drs. S. Martin and H. Ramsey were the Principal Inves- 
tigators . 


F. The University of California at San Diego built a beam-switching 
photometer system for infrared observations of solar flares. The observing 
site is at Mt. Lemmon, Arizona, with observations to extend from the 700-y 
region down to possibly Ip, utilizing the 60-inch Cassegrainian telescope. 
These infrared continuum observations should define the solar spectrum 
well enough to allow differentiating among the major theories of white- 
light flare emission. Dr. Hugh Hudson was the Principal Investigator. 

G. The LMSC at Palo Alto, California, conducted a theoretical 
study of helium visible and ultraviolet emissions from solar active regions. 
Calculations of the statistical equilibrium populations for a 19-level del 
and a 10-level Hell ion were performed. Results of these studies should 
allow the interpretation of line observations in terms of electron density 
and temperature profiles for the active regions. Dr. J.L. Kulander was 

the Principal Investigator. 
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H. The Uttar Pradesh State Observatory In India was to study 
dissociation and excitation equilibria of various molecules in the 
photosphere^ spots, and faculae. Detection equipment was supplied to 
India on loan for an observational program using their existing horizontal 
solar telescope and its associated spectrograph, with the goal of improv- 
ing available models of the solar atmosphere parameters which influence 
these studies. Dr. M. C. Pande was the Principal Investigator. No final 
report was required from this effort, although publications in the scien- 
tific journals are anticipated. 

I. The Applied Physics Laboratory (APL) of Johns Hopkins Univer- 
sity performed, during the ATM mission, an observational program of high 
time-resolution spectral observations of solar radio bursts in the 500- 

to 1000-MHz range. With the 0.01 sec time resolution available from their 
spectrograph and 60— ft diameter antenna, it should be possible to determine 
whether a finite frequency drift rate is present in the radio bursts. 

This, in turn, would affect models of the emission mechanism and electron 
densities in the source region. Mr. B. Gotwols was the Principal Inves- 
tigator . 


The nature of the SGAP program was to provide corollary data for 
postmission analyses in conjunction with the ATM. As such, it complements 
those programs, such as the National Oceanic and Atmospheric Administration 
(NOAA) effort, which were designed to provide the real-time, or near-real- 
time, data required by the ATM mission. 

The ATM Instruments were designed to take advantage of orbital 
altitude to obtain data in the X ray and ultraviolet portions of the 
spectrum where atmospheric absorption limits observations from the earth. 
The SGAP observations extended the spectral coverage of events observed 
by the ATM by utilizing atmospheric windows in the visible, microwave, 
and infrared portions of the spectrum. Figure 1.1 shows the spectral 
ranges covered in relation to each other. 
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SECTION 2. UNIVERSITY OF HAWAII MULTICHANNEL CORONAL SPECTROPHOTOMETER 


The Principal Investigator for this experiment was Dr. John T. 
Jefferies, with Drs. Frank Orrall, Jack Zlrker, and Don Landman as co- 
Investlgators . The University of Hawaii instrument measured forbidden 
line (magnetic dipole) radiation in the visible and near infrared spec~ 
tral regions from some of the same ions that emit the permitted line 
radiation observed by the ATM instruments. A number of prominence line 
intensities with temporal and spatial resolutions comparable to that of 
the ATM equipment are measured. Selected line intensities, together with 
neighboring sky and continuum background, were digitized and recorded on 
magnetic disk while the spectrophotometer scanned the image of a coronal 
region. In this way, simultaneous, monochromatic, high spatial and 
temporal resolution coronal region images were obtained. The data, re- 
duced in real-time with a PDP 11/45 computer system, were then available 
for immediate analysis. Certain forbidden line ratios are sensitive to 
electron density changes over relatively narrow temperature ranges while 
other line ratios are insensitive to density and may be used to fix the 
temperature or composition. Diagnostic methods utilizing line ratios 
have been developed and applied by the experimenters. The temperature, 
density, and composition of the active regions can thus be determined, 
thereby providing inputs for physical models of the active regions as 
functions of time and activity. These physical models can be used to 
determine the rates of energy loss by radiation and conduction, and the 
coronal energy balance can be studied. 


The apparatus was designed to record simultaneously many 0.33p to 
l.lp line intensities in order to eliminate systematic errors in obtain- 
ing line intensity ratios resulting from seeing and sky transparency 
fluctuations during a measurement. The spatial resolution obtainable by 
the instrument is on the order of 5 arc seconds. The instrument's rela- 
tively high temporal resolution capability enables following closely the 
time variations of physical properties of active regions. 

^^8ure 2.1 is an overall view of the instrument. It is roughly 
25 cm in diameter by 4.4 m long and is mounted on one face of the 3.6-m 
spar at the Mees Solar Observatory (Figure 2.2). The glmbal mounting 
allows programmed spatial scanning within 3 with the stepping motor 
translation drives designed to provide 1.66 arc second/step motions. 
Scanning is under computer control so that rectangular, circular, etc., 
rasters can be programmed and used as desired* 


The optical system, shown schematically in Figure 2.3, consists of 
a coronagraph telescope, relmaging optics, component wheel, spectro- 

shutter /sun-positioner. The detector is a thermoelec trlcally 
cooled silicon vldicon camera tube placed at the spectrometer's output 
ocal plane, located behind the spectrometer. After exposure, the dura- 
tion of which can be arbitrarily selected, the target is read out in a 
slow-scan mode (approximately 2 seconds/frame), with the video signals 
amplified, digitized, and routed via CAMAC to the computer. Any pixel 
set (up to 236 x 384 at 12 blts/plxel) can be selected by the program 
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Figure 2.2 Mees Solar Observatory 
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Figure 2.3 Optical System 



for processing. The data acquisition and control system (DACS) is config- 
ured around a PDP 11/45 computer with 24K words of memory. Peripherals 
include teletype, line printer, DEC tape, and disk units. All other inter- 
facing, such as CRT display, motor controls, etc., is via a CAMAC system to 
take advantage of data-way standardization. The DACS control functions 
Include; telescope pointing; vidicon exposure, read-out, and erasure; 
pixel selection and data transfer; and on-line reduction with associated 
CRT displays. A more complete instrumentation description is contained in 
Appendix A. 
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SECTION 3. LOCKHEED SPECTRA-SPECTROHELIOGRAPH 


Dr. Alan Title was the Principal Investigator for this effort 
involving spectra-spectrohellography (s2hG), an extension of spectro- 
heliography (SGH) . In both S^HG and SHG, the solar image is moved with 
respect to the entrance slit of a spectrograph and recorded on film which 
Is simultaneously being moved with respect to the exit slit. The most 
Important difference is that S2 gh records a wide spectral region and SHG 
records only a narrow band of wavelengths. Following is a brief descrip- 
tion of the system's opera t ion , taken largely from Dr. Title's original 
experiment proposal. ' 

Ordinarily, spectroheliograms are made by fixing both the film 
and the sun's image and moving the spectrograph. Exposure is through 
a narrow exit slit. The result is a photograph of the sun with variations 
in density representing variations in the intensity of solar radiation 
(integrated over a small fixed wavelength band). It is difficult to 
determine from a single spectroheliogram whether these changes in density 
are due to Doppler shifts, Zeeman splittings, or variations in continuum 
intensity or line strength. Various techniques using combinations of 
spectroheliograms have been developed in an attempt to distinguish 
physical effects. The S^HG system allows simultaneous collection of the 
data needed for these various studies. 

The S^HG spectrograph (Figure 3.1) has a wide exit aperture and a 
digitally controlled, stepping motor driven movie camera capable of rapidly 
advancing the film. In recording spectra, the spectrograph remains sta- 
tionary and the sun's image is allowed to drift past the entrance slit. 

A single frame of the film is exposed while the sun moves a fixed distance 
with respect to the entrance slit. The film is then rapidly advanced and 
the next frame exposed. Thus, S^hG spectra are recorded in a series of 
frames on the film, each frame containing the entire spectral profile 
along a line on the sun's disk (Figure 3.2). Since successive frames con- 
tain spectral profiles for adjacent solar regions, the film contains the 
spectral line profile for each position on a two-dimensional grid on the 
solar disk. By inserting quarter or half-wave plates and a polarization- 
dependent beam splitter into the system, one obtains spectral profiles in 
polarization pairs, as shown in the figure. 

If a line displays a Zeeman effect and the spectra have been recor- 
ded in right and left circular polarizations, longitudinal magnetic field 
maps can be constructed. If Zeeman sensitive lines are photographed in 
several pairs of linear polarizations, then transverse magnetic field maps 
can be constructed. S^HG magnetic field maps yield resolutions as high 
as one-half arc second under favorable seeing conditions. 

Velocities to about 7 meters per second absolute can be measured 
by comparing the spectra against a standard line inserted into the spectro- 
graph by a laser or atomic absorption cell. Operating in the velocity 
field mode, the S^hG observations cover a region of 1 x 2 arc minutes in 
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Figure 3.1 S^HG Data Collection System 



Figure 3.2 Typical s2rG Solar Spectra 
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less than 15 seconds. This enables studying the velocity field develop- 
ment in regions of substantial physical extent. Using beam splitters 
and folding optics allows simultaneous recording of spectra of two or 
more lines formed at different heights for studies of changes in phase 
and amplitude with height. 

2 

The S^HG cameras were mounted on Instruments at the McMath Solar 
. Tower at Kitt Peak National Observatory in Arizona. At Kltt Peak the 
large spectrograph has a dispersion of 0.14 A/mm at A5250 (5th order), 
so that 20vi corresponds to 0.00276 t at X5250. A shift of 20 microns 
corresponds to 71 gauss or to 157 meters /second. The location of the 
peaks can be reliably determined to a micron by comparison with adjacent 
points, so that fields could be measured to an accuracy of 3.5 gauss and 
velocity of 7 meters/ second. 

Due to a fortunate coincidence, fifth order Fe A5250.2 (5 x 5250.2 
= 26251.0) and fourth order H-alpha X6562.8 (4 x 6562.8 = 26251.2) are 
nearly coincident in the spectrograph focal plane. Hence, a dichroic 
beam splitter that separates red and green light and the normal polariza- 
tion beam splitter used for magnetic field studies make it possible to 
simultaneously study the spectral region surrounding H-alpha and the 
magnetic field in X5250. When S2 hG is operating in the H-alpha mode, it 
is possible to conveniently study a spectral region 3.3 X wide, centered 
on H— alpha, which allows following velocities in H— alpha for + 75 fcra/ 
second. ” 


The S2 hg film is scanned by a Photometric Data Systems Model lOlOG 
microdensitometer whose output is encoded digitally and fed to a PDP 11/45 
computer. The computer also provides control for the scanning functions. 
The microdensitometer is capable of scanning some 6000 points per second, 
with positional encoding to an accuracy of one micron. 

Appendix B is Section I of Dr. Title’s final report, containing 
discussions of theory and results. Sections II and III of the report are 
not included herein, but they contain operation manuals and programming 
Information for the microdensitometer and computer system, should outside 
users be interested. 
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Figure 3.3 S^HG Data Reduction System 



SECTION 4, NBS RADIOMETRIC CALIBRATION OF ATM EXPERIMENTS 


The National Bureau of Standards participated in the SGAP program, 
although not for ground-based observations per se, but rather to enhance 
their capabilities to provide calibration services to the ATM and its 
associated rocket programs in the ultraviolet wavelength region ♦ 

The NBS effort was divided into three distinct phases as follows: 

A. Continued development work on a hydrogen arc plasma source of 
known radiant flux for the calibration of spectrometrlc-detector systems 
over the wavelength region of 500 to 3700 A* 

B. The second part of the program was a study of the effect of 
ultra-high vacuum and particular contaminants on the sensitivity of 
photocathodes in the windowless region of the vacuum ultraviolet (VUV) . 
These investigations of the cause of observed changes in windowless 
photomultiplier response in space help establish a greater integrity of 
calibration from the laboratory to a space environment* 

C. The third part of the program was the development of cali- 
bration capability at the NBS Synchrotron Ultraviolet Radiation Facility 
(SURF) for the purpose of extending radiometric calibrations to shorter 
wavelengths* For certain applications the calculabllity of this source 
can be utilized for direct calibrations* Also, work was done toward the 
development of diode detectors as transfer standards for the wavelength 
region shorter than 600 A. 

The H 2 arc is capable of producing accurate continuous emission 
in the VUV at reasonable temperatures. The known frequency distribution 
of the continuum emitted from a pure hydrogen plasma, with the eraissivlty 
exactly calculable, allows calibrations using the continuum intensities 
at any wavelength, with accuracies dependent on the accuracy of the plasma 
diagnostics* Extensive numerical calculations of the continuous emission 
of hydrogen from 950 A to 3650 X have been performed by NBS under the 
assumption that a hydrogen arc plasma is in local thermal equilibrium 
(LTE) and at atmospheric pressure* The accuracy of the calculations is 
estimated to be within 2 percent. The main uncertainty originates from 
the high density plasma corrections. 

The arc was initially operated at 13,000 ‘’K, being an excellent 
standard down to about 1600 A where molecular Lyman band emission begins* 
Atomic Jine radiation is also present from a few Lyman lines in the 1216 X 
to 950 A region. Modification of the arc for operation at higher tempera- 
tures up to 25,000 allowed an extension of the calibration range down 
to 1100 A in that the molecular line emission became insignificant due to 
dissociations. 



The arc is wall-stabilized by a stack of electrically insulated, 
water-cooled disks with an axial bore of a few millimeters. Ring elec- 
trodes allow emission to be observed axially. Provision is also made 
for side-on observations. For wavelengths below 1100 X, the LiF window 
is removed, requiring a differential pumping system. Measurements are 
then possible down to about 500 A, although Inaccuracy increases as wave- 
length decreases. 

Utilizing the H 2 arc, calibration should be reliable to an ac- 
curacy of about 15 percent from 1650 A through 6000 A. Inaccuracy will 
increase as wavelength decreases to about 25 or 30 percent at 500 A wave- 
length. 

The photocathode contamination study’s objective was to understand 
the effects on the performance of photocathodes caused either by removal 
or addition of surface monolayers of contaminating substances. The tests 

vacua of 10“^® torr, over the wavelength region 500 to 
1500 A, and include the energy analysis of the photo— emitted electrons for 
diagnostic purposes, A number of contaminating gases were used with the 
element tungsten as a VUV calibration detector on the photocathode. Tung- 
sten was used since its physical properties are well known. Induced con- 
taminants included the common gases 0,, H„, N„, etc., and such things as 
methanol, aldehydes, acids, esters and ketones. 

Surfaces were prepared in the atomically— clean state and, subse- 
quently, exposed to controlled amounts of known contaminating gases 
measuring concurrently the changes in photoelectron yield. The main 
constituents evolved during the outgassing of space plastics were of 
relatively high molecular weights. However, contamination effects of 
these species can probably be partially understood by a study of these 
related lighter hydrocarbons. During these experiments, the energies of 
the photoelectrons were measured to determine both the cause of the 
changes in yield which are observed and the origin of the photoelectrons 
(e.g., the photocathode proper or the contaminating layer). 


Results of these studies have application to contamination of such 
pUotocathodes used as photomultipliers for the far (windowless) ultraviolet 
when exposed to outgassing of various space vehicle sources and to the 
"clean-up" of such surfaces as the exposure to the hard vacuum of space 
continues. 

NBS has a unique facility which allows the sample photocathode to 
be exposed to far UV radiation in the windowless region without destroying 
the required hard vacuum. It consists of a gas discharge light source 
connected through a two— stage differential pumping unit to a high vacuum- 
normal -incidence VUV monochromator. This monochromator, of stainless steel 
with metal 0-rings, has a 10-lnch diffusion pump. The bakable experimental 
chamber fastens to the exit slit of the monochromator and, with a high- 
speed pump and LN 2 baffle, maintains 1 x 10“^® torr with the light source 
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running. A sample manipulator allows two samples to be studied simul- 
taneously and independently heated for cleaning. An electron energy 
analyzer la a part of the main experimental unit. 


The third part of the NBS program was to establish a capability 
or the radiometric calibration of ATM experiments from 600 1 to 200 X 
by utilization of the NBS Synchrotron Ultraviolet Radiation Facility 
^ (SURF). The approach was two-pronged: First, calibration procedures 

were developed using synchrotron radiation as a calculable standard 
source of flux for experiments where this approach is feasible; and 
second, synchrotron radiation was utilized as a continuum light source 
with a monochromator and an absolute detector for the calibration of 
transfer standard diode detectors. 


u 4 ®5^^J^^'^tron radiation is a strong source of continuum radiation 
which is highly collimated (f /360) , linearly polarized and of calculable 
ntenslty, assuming that the number of electrons circulating (electron 
current) is known. This effort Involved determination of the electron 

a measure of the proportional radiation intensity, and careful 
definition of the photon beam collection geometry. Also Involved was the 
fitting of a test port with a grazing incidence monochromator which is 
rotatable about the optic axis— required by the polarization of the beam 

nd optical systems to be tested— and the associated test fixtures and 
electronics . 


The high collimation of synchrotron radiation requires optical 
angular aperture to be scanned. These problems may 
restrict the instruments for which this calibration technique is 
applicable. In the second approach SURF acts as the continuum source of 
ThI the calibration of transfer standard detectors, 

for til T"’ photodiodes, one with a window 


for 1 Iv *^'^7 : Piiocoaioaes, one with a window 

for the wavelength region of 1200 ^ up to 2400^X. and a windowless one 

it 1200 A down to below 600 A. This effort involved 


extending from about 

development of a transfer standard to extend calibration to shorter 
wavelengths. Work concentrated on AI 2 O- photocathodes because of their 
comparative stability of photoelectrol Jleld with temperature. The 

testing is a specially constructed double- 

lenEt^^ad^2^ transparent to this short wave- 

radiation, hence, the oxide is used to eliminate natural oxida- 
tion and resultant change of photoelectric yield with time. 
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SECTION 5. CAL TECH VACUUM PHOTOHELIOGRAPH 


The California Institue of Technology had a twofold effort in the 
SGAP program; mounting the 26-inch photoheliograph functional verification 
unit, including its own specially constructed vacuum tank onto the telescope 
mount yoke at Big Bear Lake Observatory, and the construction and placing 
into operation of a special new birefringent filter for observations in the 
10830 A line of helium. Dr. Harold Zirin was the Principal Investigator 
for this effort, and Richard Proutt was the experiment manager. 

Among the necessary efforts in the study of Solar Physics is the 
observation and understanding of successively smaller structural detail on 
the sun's surface. The Big Bear Lake site was selected Initially for its 
excellent "seeing" conditions and the fact that some of the best resolution 
work to date had come from that site. Installation of the 26-lnch photo- 
heliograph was intended to take advantage of this "seeing" by providing 
good optics with a large aperture. 

The 26-inch photoheliograph is a Gregorian telescope which was 
expressly designed for future space flights. The instrument (Figure 5.1) 
was built as a functional verification unit with special attention to 
reducing the effect of solar heating on the various mirrors. The para- 
boloidal primary mirror is coated with a special high reflectance silver 
coating which greatly reduces thermal distortion of the shape of the pri- 
mary mirror. An image of the sun is formed at a diagonal heat stop mirror 
which rejects all the solar energy except that from a small field approx- 
imately 5 arc minutes on the side which is only 3.3 percent of the sun’s 
area. This image is then enlarged by an ellipsoidal secondary mirror and 
reflected by two diagonal mirrors which direct the image to a location 
where it may be studied (Figure 5.2). The entire instrument, up to the 
cameras, is enclosed in a vacuum tank with a vacuum entrance window and 
exit windows so that all effects due to air currents inside the telescope 
structure may be eliminated. The vacuum tank was mounted in the existing 
fork mounting at the Big Bear Solar Observatory, along with two other 
vacuum telescopes, a 10-inch refractor and an 8.6-inch refractor, which 
simulated the ATM hydrogen-alpha telescope. The beam from the 26-inch 
telescope may either be utilized directly or may be diverted through a pair 
of mirrors to the Coude focus of the telescope, where a powerful Littrow 
spectrograph enables spectrograms to be made of various selected features 
on the sun. Alternate observations using narrow-band birefringent filters 
with spectrograph observations and observations using a video magnetograph 
allow the determination of magnetic fields in near real time. Another 
instrument used with the 26 -inch system is the universal birefringent 
filter which enables observations in any spectral line from 4000 to 7000 
angstroms. These systems obtain long-term observations of the various 
features that were observed from the ATM so that their background, history, 
and fine structure can be better understood. The 26 -inch photoheliograph 
is theoretically capable of resolution down to 0.2 arc seconds at H-alpha. 


17 



SECONDARY MIRROR 
ASSEMBLY 


INTEGRATING STRUCTURE 
ASSEMBLY 


PRIMARY MIRROR 
ASSEMBLY 


Figure 5.1 26 -Inch Photoheliograph Assembly 







Figure 5*2 Optical Layout 

An automatic alignment system has been provided which is designed to keep 
the telescope in good alignment at all times despite changes in tempera-- 
ture and the like. Many of the apparent larger scale phenomena that are 
observed with lower resolution originate in smaller features* In partic- 
ular, solar granulation is made up of structures only a few hundred kilo- 
meters across, yet these are thought to be one of the main mechanisms by 
which the chromosphere and corona of the sun are heated* Similarly, solar 
spicules, which are only a few hundred kilometers across, extend upward 
into the atmosphere and are thought to play an important role in the 
transition between the chromosphere and the corona. Inside sunspots many 
features which are near the limit of resolution appear to play an impor- 
tant role in the growth of magnetic fields and the occurrence of flares* 
Therefore, high resolution observations of what is occurring inside these 
features give physical Insight into what is happening and complement and 
extend the ATM data even though the features that were observed by ATM 
scientific experiments were, in general, more coarse than those which are 
resolved by the system discussed here. 

o 

A 10830 A birefringent, narrow band filter was especially designed 
and constructed for observations at this infrared wavelength, which is the 
only helium line observable on the disk of the sun. Since helium Is 
difficult to excite, it is seen only where the temperature is high. Thus, 
it provides an opportunity to observe the chromosphere on the disk without 
interference from photospheric radiation below, since the photosphere is 


19 


cooler. Because of problems with the filter, no simultaneous observations 
could be made with ATM, although postmission observations will be made and 
reported. These observations are complementary to the extreme ultraviolet 
(XUV) data from the ATM experiments in that they also look at the hotter 
part of the solar atmosphere. The instrumental setup for 19830 I observa- 
tions includes, in addition to the filter and telescope, a C0« cooled RCA 
infrared image converter tube and a 35mm camera for recording^the image. 

Appendix C is the Cal Tech final report describing installation 
and initial results of the photoheliograph at Big Bear Observatory. 



Figure 5.3 Big Bear Lake Observatory 
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SECTION 6 * LOCKHEED FILTER STUDIES IN Hel AND Call AT RYE CANYON 


The Lockheed SGAP effort included research programs involving disk 
^nd limb observations in the D^ line of Hel and high resolution observa- 
tions in the Infrared Call line at 8542 X, The two programs had comple- 
mentary objectives. Drs. Sara Martin and Harvey Ramsey were the Principal 
Investigators* 

During periods of high solar activity on the disk, the Hel observa- 
tions were directed toward photographing solar flares and transient events* 
Effort was made to take these D^ observations simultaneously with ATM 
spectra and spectroheliograms at X ray and XUV wavelengths. During periods 
of low disk activity, the effort was concentrated on acquiring prominence 
observations for correlation with ATM coronagraph observations* 

The Infrared Call observations were high resolution, time-lapse 
photographs of only selected active centers or comparable sized areas of 
the solar surface* These obseivations were taken during the ATM observ- 
ing periods (primarily the second manned mission) as well as several days 
after the ATM observations in order to acquire a reasonable chronology of 
the active regions for which X ray and XUV spectroheliograms were taken 
from space. The Call and Hel filtergrams taken in conjunction with the 
ATM X ray and XUV spectroheliograms permit study of the origin of features 
at different levels in the solar atmosphere. 

2 The Hel program consisted of building and testing a filter for the 
D 3 (2 P 3 D) neutral helium line at 5876 K and placing it in operation on 
one of two solar spar telescope systems at the Rye Canyon Observatory 
(Figure 6*1). The filter was mounted on the 28-inch diameter, 12—foot long 
spar which contains four optical systems, each occupying one quadrant of a 
cross fo^ed by radial struts around a 5-inch center tube which contains 
the guiding optics. The ^narrow band 0.4 A filter is used with a 7 A block- 
ing filter. A second 7 A blocking filter was mounted in a full disk occult- 
ii ^8 system for studies of limb prominences. The narrow band filter permits 
referencing of oppositely polarized transmission peaks of 0.4 A separated 
by 0.8 A which allows referencing centerline D 3 and continuum at + 0.8 t 
or determining Doppler differences at + 0.4 A subtraction techniques. The 
line is seen to appear in either emission or absorption against the photo- 
sphere during flares; the strength relative to the continuum is used to 
infer electron density in the emitting volume. 

Mounted <gn the other solar spar system at Rye Canyon was the optical 
system and 0.25 A tunable filter for the infrared Call line high resolution 
observations at 8542 A. This solar spar system is slightly larger (14 -foot 
length and 36-inch diameter) but otherwise very similar to the other spar. 
The Call line is believed to originate at a higher level in the solar 
atmosphere than H-alpha. Therefore, these observations should prove 
extremely valuable in relating the ATM X ray and XUV data to H-alpha filter- 
grams and to spectroheliograms taken at wavelengths originating at lower 
levels . 
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Figure 6.1 Lockheed Solar Spars 


TUis task, provided the photography of Call features with excellent 
spatial resolution and permitted the recording of all except the most 
rapid changes in solar flares and accompanying transient events. The 
optical system included, in addition to the 0.25 X birefringent filter 
v^th polarizing beam splitters, a 10-lnch objective lens and a broadband 
blocking filter mounted on the track which served as the optical bench, 
as well as the camera system for data taking on infrared film. A visual 
monitoring system consisting of a conventional television camera and 
monitor was provided for operator convenience’. 

Instrumentation and observations are described in Appendix D. 
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SECTION 7. UNIVERSITY OF CALIFORNIA AT SAN DIEGO (UCSD) SOLAR 

INFRARED OBSERVATIONS 


The objective of the UCSD observations was to explore the last 
remaining unknown region of the solar flare electromagnetic spectrum^ 
the infrared wavelengths between Ivim and 1 mm. Flares appear prominently 
at almost all other wavelengths, and historically the development of ob- 
servational techniques for each wavelength has revealed phenomena of major 
Importance. Indeed, the X ray flare observations within the last decade 
have completely changed the accepted picture of solar flare energetics. 

Atmospheric opacity limits ground-based observations to certain 
regions of the infrared spectrum. Out to about 25u one finds numerous 
transparent "windows**; between 25p and lOOy, however, the sun would not 
be detectable even at excellent mountaintop sites, Longward of lOOp, 
atmospheric transmission windows are narrow and require narrow-band 
detection to exploit them fully. In the present series of observations 
the 25p and 300 to 1000 p regions were explored. Observations were primarily 
from Mt. Lemmon, Arizona, although a two-week run was scheduled at Mauna 
Kea, Hawaii. 

Infrared detection generally Involves **sky cancellation** to help 
remove the effects of emissions from the atmosphere and various parts of 
the telescope. With the 60-inch Mt- Lemmon telescope at submllliraeter 
wavelengths, angular resolution is good enough to permit two beams on 
the solar disk. Chopping between them cancels both the background signal 
and the quiet solar signal, so that a source of transient activity In one 
beam (e.g., a flare) can be recognized by Its time profile. At the longer 
wavelength end of the Infrared spectrum, phenomena similar to those in the 
millimeter-wave region of radio observations were expected. In this region 
both impulse and gradual bursts are found, arising apparently from distinct 
emission mechanisms. The nature of the emission mechanism of the impulsive 
component remains unclear, and the submillimeter observations should help 
in understanding it. In the 1 to 30pm range, a most interesting phenomenon 
is the **white-light** flare. Understanding the white-light flare phenomenon 
potentially has very great Importance, partly because of Its close relation- 
ship with the energy release in the **flash phase** of the flare, and also 
because the visible continuum may also represent a large fraction of the 
flare energy. Infrared observations provide the crucial test to decide 
among different emission mechanisms. Solar flare Infrared continuum might 
arise from any layer of the solar atmosphere above the level of the temp- 
erature minimum, below which the optical depth becomes large. This region 
Includes the site of H— ot flare emission. Possible emission processes would 
include white-light continuum, thermal continuum from the H-ot source, ther- 
mal continuum from hot regions, non— thermal Bremsstrahlung, and synchro- 
tron emission, A great deal of the justification for this effort was the 
possibility, however slight, for obtaining data on a white light flare in 
order to choose between these mechanisms. Although this did not material- 
ize, numerous subflares were observed, and an Important new phenomenon was 
discovered: temperature fluctuations associated with the 300-second oscil- 

lations of the solar atmosphere. Appendix E describes this discovery. 
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Figure 7.1 UCSD Telescope at Mt. Lemmon, Arizona 
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SECTION 8. LMSC STUDIES OF HELIUM EMISSION IN VISIBLE AND UV 


This was a theoretical study by Dr. J.L. Kulander of Lockheed, the 
objective of which was to develop a better understanding of active solar 
region phenomena such as flares, prominences, and filaments through 
observations of visible and ultraviolet eralSBlons from helium in the 
regions of interest. 

The study provided the framework for the Interpretation of and 
correlation between the various emission and absorption features found 
in active solar phenomena. A detailed calculation of the statistical 
equilibrium populations for a 19 level Hel atom and 10 level Hell ion 
was conducted. The results are used to^develop relationships between 
the visible D3 line emission (and 4686 A line of Hell) and Important UV 
resonance line emissions such as the 584 A and 537 X lines of Hel for a 
wide variety of temperatures and densities as would be anticipated in 
solar active regions. The specific interpretation of visible Hel Do line 
observations' provided estimates of electron density and temperature in 
solar active regions. 

A theoretical study of the visible and UV line radiation of Hel 
atoms and Hell ions from flare and flare associated regions was made. 
Primary emphasis was placed upon the Interrelations between the Hel - Do 
line (5876 A) and the first two resonance lines of Hel (584, 537 I). This 
emphasis resulted from the planned program to obtain D3 flltergrams on a 
patrol basis at the Lockheed Rye Canyon Observatory during the ATM mission. 
The calculations include other visible and UV lines and continue. The 
purpose of the program was to obtain a theoretical understanding which can 
be used to interpret the simultaneous UV measurements from the ATM and the 
ground-based visible line measurements. 

A code was developed to calculate the radiation fields and popu- 
lations of the individual terms in Hel and II through principal quantum 
number n = 4. A plane-parallel model was used, with solutions analyzed 
for electron densities IQlO to IQl^ cm’3 and electron temperature from 

c Optical depths and some line intensities are presented 

for a 1000 kilometer thick layer. 

Appendix F is the final report for this study effort. 
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SECTION 9. MOLECULAR LINE OBSERVATIONS AT UTTAR PRADESH, INDIA 


The Uttar Pradesh State Observatory (UPSO) In India is doing a 
study of dissociation and excitation equilibria of various molecules 
in the photosphere, spots, and faculae, using data obtained during the 
Sky lab mission. Detection equipment was supplied to India on loan for 
an observational program at the Manora Peak facility using their existing 
horizontal solar telescope and its associated spectrograph with a view 
toward improving available models of the solar atmosphere’s parameters 
which influence these studies. Dr. M.C. Pande was the Principal Inves- 
tigator . 

The solar group at UPSO is investigating the following types of 
problems concerning the sun during the ATM mission: 

A. A comparative study of the dissociation and excitation equi- 
libria of various molecules in the photosphere, spots and faculae. 
Theoretical studies concerning these have previously been undertaken 
based on various available models of these features. The present effort 
is based on an observational program using a spectrograph in conjunction 
with the horizontal telescope. 

B. Center-to-llmb variations of the atomic and molecular lines 
in the spectra of photospheric and active regions and their intercompar- 
Ison. 


The equipment complement of the horizontal solar telescope and its 
associated spectrograph includes: 

A. Image-forming optics, Including a 46-cm coelostat and a 27-cm, 
f/66 skew Cassegrain telescope capable of forming an 18-cm image of the 
sun. 


B. A spectrograph employing a 127 x 203 mm^ plane reflection 
grating with 600 lines/mm blazed at 2.5)i and a 25-cra, f/36 camera giving 
a dispersion of 1.2 i/mm in the first order. In this program spectral 
orders up to the fifth order were used. The recording facilities were 
provided for photoelectric spectrophotometry (with scanning facilities 
in the fourth order at speeds varying between 1 X in 768 seconds to 120 
X/minute) and photography. Simultaneous spectrophotometry in different 
spectral regions, not exceeding four in number, is possible, in the 
region 3000 to 9000 A. 

Instrumentation which was placed on loan to Uttar Pradesh for 
observation in the ATM time frame includes: 

A. Refrigerated photomultiplier tube chambers to reduce dark 
current noise in the detection. 



B« PbS detectors to extend the spectral coverage Into the 
Infrared. 

C. An Image Intenslfier tube to allow prominence and weak 
emission line studies. 

D. Film magazine adaptable to the image tube, film, and associated 
electronics . 

Although no report to NASA was required for this equipment loan, 
publications in the open literature are expected covering the Sky lab and 
subsequent time frame observations. The basic Skylab program was center- 
to-llmb variations of Mg trlplett, Na, Di and D 2 lines, C 2 bands, MgH, 
and CH 4300 A region on each clear day for photospherlc and spot spectra. 
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SECTION 10. JOHNS HOPKINS RADIO BURST SPECTRAL OBSERVATIONS 


The Applied Physics Laboratory (APL) of Johns Hopkins University 
performed an observational program during the ATM mission of high time 
resolution spectral observations of solar radio bursts in the 500 to 
1000 MHz range. Bruce Gotwols was the Principal Investigator. With the 
high time resolution available from their spectrograph and 60-foot dia- 
meter antenna, the observational goal was to determine the presence or 
absence of frequency drift in solar radio bursts and associated impli- 
cations on models of the emission mechanisms and electron densities in 
the source region. Observations during ATM were principally a continua- 
tion of earlier work; hence, no instrumentation development was required 
for this project. 

A principal goal of the 565 to 1000 MHz spectral observations was 
to utilize the extremely high 0.01 second time resolution of the APL 
spectrum analyzer to determine if, in fact, there is a finite frequency 
drift rate in the radio bursts from the sun. A finite drift rate would 
strongly imply an emission mechanism operating at or near the local 
plasma or gyro frequency. Electrons travelling radially outward from 
the photosphere with velocities of appreciable fractions of the speed 
of light should traverse the distance from the 1000 MHz to 500 MHz plasma 
levels on the order of a tenth of a second, thus allowing an unambiguous 
conclusion as to the presence or absence of a finite frequency drift rate. 

It is also possible to obtain a time history of the activity at 
any selected frequency within the passband of the spectrum analyzer. 

From the lifetime of a drifting burst at a particular frequency, the 
surrounding plasma temperature may be determined, and from the frequency 
drift rate the density gradient can be estimated. Thus, data are 
provided for modeling the physical parameters of the corona -chromosphere 
transition region where at least some of these bursts are thought to 
originate. 

The antenna is a 60-foot diameter paraboloid mounted on an X-Y 
mount with computer generated pointing control. Focal length is 25-feet, 
at which a cavity backed spiral, with right-hand circular polarization is 
located. 

The radio frequency section of the receiver is antenna mounted, 
behind the paraboloid, as is the intermediate frequency amplifier. The 
remaining portions of the spectrograph are located in an instrumentation 
trailer adjacent to the antenna. Figure 10.1 is a block diagram of the 
complete spectrum analyzer. The minimum detectable signal is approx- 
imately 20 X 10~22 watts/m2/Hz. 

The data recording was accomplished by letting the log amplifier 
output modulate the intensity of the spot on a high resolution cathode 
ray tube. Simultaneously the spot is driven laterally in synchronization 
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Figure 10.1 Block Diagram of the 565-1000 MHz High Time Resolution Solar Spectrum Analyzer 
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with the local oscillator tuning of the spectrum analyzer and photo- 
graphed with a continuously moving 35mm camera film. Figure 10.2 is a 
typical recording. The spectrum analyzer output was also recorded on 
magnetic tape in analog form. Appendix G discusses results of this task. 



Figure 10.2 Solar Radio Burst Spectrum 
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Figure 10.3 APL Sixty Foot Radio Burst Antenna 
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INTRODUCTICN 


The nxulti-chsrmel coronal spectropliotcmeter system described in this 
report was conceived as part of the 1973 SKYLAB mission ground-based solar 
observational support. Much of the, Apollo Telescope t-tount instrunrjnt 
package aboard SrCYIAB was configured especially to investigate coronal 
physics via high spatial and terporal resolution observations of the line 
radiation in the XUV and EUV spectral regions. Our apparatus was designed 
to ccnplenent that v;ork by measurd.ng, in coordination with the SKYLAB 
scientific schedule, the corresponding forbidden line radiation in tlie 
visible and near-visible spectral regions, together with a number of prcmi- 
nence line intensities, with resolutions ccnparable to that of the ATIi 
equipment. 

The instrument was in op>eration on a daily basis throughout the mission. 
Ha«;ever, due principally to limiting signal-to-noxse capability, sparse 
coronal activity and spotty weather conditions, our observationsal program 
was only partially successful. By far, our rrost scientifically significant 
results were obtained during the second manned mission. V'Jhen appropriate, 
supplementary observations were made using our other existing instrumentation; 
10-inch coronagraph-coude spectrograph and Ha Zeiss system, Ha flare patrol 
camera, K-line telescope, Ha and green-iine dual coronagraphs , and polar imeter- 
photameter , 

Two papers describing the multi-channel coronal spectrophotcmeter system 
have been presented; one was published in Astroncmical Observations \'ith 
Television-T\fpe Sensors , University of British Columbia (1973) , and the other 
was delivered at the Solar Physics Division^ meeting, of the American Astrono- 
mical Society, 9-11 January 1974 (paper P4) . 

Anticimting the analysis of the coronal observations, we also carried 
cut neiv Fe'^i2 excitation equilibrium calculations based on our improved ground 
configuration proton iitpact excitation rates. Three p^jers describing this work 
have been published in Solar Physics . 

The succeeding sections to this report describe, in turn, the multi- 
channel coronal spectrophotometer system's design and operation, the obser- 
vations made \'dth the system, and the proton ijipacf excitation calculations . 

The Appendix contains copies of the scientific papers resulting from the work 
performed under this contract. 

Photographs of the spar-mounted instrument (with the two-dimensional 
vidicon detector unit) and the data-aoquisition and control systen are shasn 
in Figures 1 and 2, respectively. 
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II. INSTFUMEOT DESIQJ M© OPEIWriON 


The projected scientific programs constrained the instrunent system 
design in several directions: 

(a) To have the facility of measuring simultaneously as many coronal and 
prcminence emission lines as possible^ the apparatus had to be sensitive 
and achronatic over a very wide spectral range. In particular, we desired 
the capability of measuring simultaneously the three Fe+12 lines: AX3388, 

10747 and 10798. For obtaining line intensity ratios, the feature of 
simultaneity e] iminates any systonatic errors resulting frcrti sJq^ transpa” 
rency and seeing fluctuations throughout a measurar^t; 

(b) To be able to derive meaningful physical models fron the data, the instru- 
ment required a spatial resolution sufficient to resolve the intrinsic fila- 
mentary structure of coronal active regions (typical diameters >5^^) ; 

(c) To enable close time-variation observations of the physical properties of 
active rcjgions, the instrument was designed with as high a temporal resolution 
capability as possible; and 

(d) To maximize scientific throughput a data-aoquisition and control system 
with on-line reduction and display facilities was required. To allow for the 
possible coordinated operation of several solar telescope systems it was 
desired that the DACS have multiprograirining capability. 

The follaving subsections describe in turn the various apparatus subsystem 
design characteristics and functioning. 

A. MECHANICAL CCNFIGURATIO^S 

Figure 3 shews an overall imchanical view of the instrument. It is -v25 cm 
in diameter by M.4 m long, constructed principally of aluminum, and mounted 
on one face of the 3.6-m spar at the Mees Solar Observatory. It consists suc- 
cessively of an on-axis coronagraph telescope, a spectrometer, and a silicon- 
vidii^n camera tube detector unit. The gimbal mounting permits 5/6 ^/step 
spatial scanning via corputer driven stepping motors. The limits of travel 
^ ^ from the sun center. Various component motions (dust cover, shutter, 
diffraction grating, optical coiponent viieel and objective lens focus) are 
driven by computer-controlled synchror»ous motor (cf . Figure 9) * 
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B. OPnCTiL CONFIGURATION 


The cptical system, shavn schematically in Figure 4, is conprised of an 
achromatic coronagraph telescope, reimaging lens, optical ooirponent viieel, 
calibration light sources, shutter/sun-positioner and spectrareter. 

(1) Coronagraph Telescope 

As indicated in Figure 4, the coronagraph telescope is used on-axis so 
that viien observing the corona, the disk's image is occulted by the side. of 
the pinhole asserrbly. Scattered light in the field is nunimized by apprcpriate 
baffling in the telescope tubes. 

I 

The objective lens is a three-element (FK-5AZFS-N4/F-7) oiled (DC 510, 
50CS) apochrcmat v.dth a 15.2-cm v.’or3ctng aperture and 220-on nominal back focal 
length (i.e., f/14.5). It is color corrected (with no aspheric surfaces) for 
out of focus" to better than tv’/o Rayleigh limits throughout the spectral range 
0. 3n<X<l.ly. 

The pinhole acts as a field lens to "form an inage" of the 0^ aperture 
on the Lyot stop located just in front of the reimaging lens. (Only*” the Op 
apertaire diffraction image was masked out; it did not appear necessary to step 
the images frem Oi surface reflections.} The most effective pinhole diameter 
corresponds to OO^'; diffraction from smaller pinholes increases the stray 
light in the field, v/hile the larger the pinhole, the less close in to the 
solar limb one can observe. The reason for the latter circumstance is that 
the instrument's aperture stop is the spectrometer entrance slit and photo- 
^heric light passed by the pinhole is scattered into the field by the re- 
imaging optics. A set of interchangeable pinholes with assorted diameters 
vere fabricated; most ebservations were made with the 30'" and 40" ones. 

A value of 'v30x10“6t was estimated for the 0]^ scatter using a bench 
test after initial assembly of the lens, when presurrably it was optimally clean. 
In actual operation, the estimated scattering after average use v;as about 
tidce the above value. Hiis result was obtained by conparing multi-channel 
coronal spectrophotareter ebservatiens of the green coronal anission line 
Fe+13 5303 a with those made simultaneously by the coronal polarimeter- 
photemeter (also noun ted on the spar) . During one cleaning operation (on 
23 June 1973) , the oil voided out beto;een the lens elements. ^In the process 
of reasseirbly a relatively deep scratch developed in the rear junction and for 
subsequent ebservations a smell step v/as placed at the position of the Lyot 
stop to mask out the image of the scratch. With, this stop in place, the 
scattering figure did not appear to be degraded over tlie earlier value. 

(2) Reimaging Lens 

This lens is used to relay the coronal light to the spectremeter entrance 
slit and, as menticned above, supports the Lyot stop on its face. It is a 
triplet (BK- 6/fluorite lBK-6) achromat designed to produce with minimal 
aberration an image v/ith 5" equivalent to 25y, a scale determined principally 
fcy the vidicon read beam diameter. 
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(3) Optical CCnroonent V?heel 

Ihis rotatable v4ieel is used to position under progratmable control a 
selection of assorted optical ocnponents just in front of the spectrorreter 
entrance slit. These ocqponents, equally spaced in eight slots around the 
wheel, include polarizers, spectral filters, neutral density filters used 
for system intensity calibration via the intensity at the center of .the sun's 
disk, and a srrall mirror to direct li^t from the calibration lairps through 
the spectrometer. 

(4) Calibration Light Sources 

Provision was made for the concurrent mounting of two types of spectral 
calibration lanps; pencil disdiarge and hollow cathode, and either type can 
be selected by means of a flippable mirror. These lanps are used for system 
spectral calibration. 

(5) Shutte r/Sun-Positioner 

The coiputer controlled shutter is situated right behind the spectroneter 
entrance slit. The shutter's front faoe is polished and angled so that when 
the shutter cones on-axis, light is deflected onto a photo-diode. The anpli- 
fied signal is read by the cerputer through a 12-bit ADC (cf. Figure 7) and 
is used for prograimvable sun-centering of tlie telescope axis by monitoring the 
gross change in light when crossing the limb. 

(6) Spectronster 

The spectraneter, which is '''66 cm long, is shewn sdioratically in Figure 5. 
The collinetor and camera mirrors are f/8 off-axis paraboloidal sections v;ith 
focal lengths of -vei on. The diffraction grating (75 groove/rmi, 6^28' blaze 
angle, 3.0y 1®^ order Littra.-; blaze wavelength) produces dispersion in the 
plane of the figure while the prism/mirror separates the overlapping grating 
orders by providing dispersion perpendicular to this plane. The spectrum is 
thereby presented to the vidicon target in the form of bands corresponding 
to the Iw orders in an echelle-type format. Using orders 3 to 9 and remaining 
near the blaze, the entire spectral range of interest '\,0.3u-l.lu can be covered, 
with much overlap, on the active target area with three grating orientations. 
The spectrometer has a low dispersion ('Ua/25u in 6th order (vslOOA)) so that 
with the target formatting used (■^33yx33y pixels) only line intensity and not 
line shape information can be obtained. As noted above, the spectroneter ^ 
entrance slit defines the aperture stop and for all data were fixed at 5"x20"; 
with respect to the sun the slit orientation remains fixed, independent of 
telescope pointing. 

C. DETECTOR 

T\vo versions of the detector unit were used during the research period. 

The active component in each case was a silicon-vidicon camera tube, chosen 
because of its extremely wide spectral response ('^0.3y-l.lu) and relatively 
hi^ quantum efficiency over this range. RCA 4532 tube with filament light 
shield w^^s used throughout. The t\±>e's target was situated in the spectraneter 
output focal plane. Bia earlier less versatile version v;as used throughout 
the first feco manned fligiits; tire final unit was in operation subsequently. 

Vfe describe the latter unit first. 
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liie final detec±or assorbly is shov/n in outline fom in Figiare 6. The Unit is 
tlierrnoelectrically cooled to '^-6° C to reduce leakage current and enable ex- 
tended exjxjsure tiires . {The fan mounted on the back of the unit v/as disconnected 
because it introduced sane noise into the video signal. This resulted in a 
yidioon operating terperature '^9^ c warmer than we had anticipated with correspond- 
^gly larger leakage currents.) An evacuated double quartz windo-v'' is mounted 
m front of the tube to prevent frosting on the tube's entrance v;indo;-/. In 
^^fstion, light integration times for v.oak emission lines in the corona are 
limited to '^iSsec by leakage current and scattered li^it in the field. 

A block diagram of the detector control electronics is shown in Figure 7. The 
vi(^con control unit at the ri^it supplies all the camera tube ooerating voltages 
(filament, grids, forms, etc.) as v,ell as the control pulses {clock, frame and 
line retrace) . The vidicon is operated in a slov'/-scan mode, the scan rate 
being determined by the video digitization rate together v/itli oonputer house- 
keeping time. ^ Horizontal cuid vertical deflection currents are generated via 
a crystal oscillator driving scalers connected to DACs. The fast axis of the 
scan consists of 256 pi>:els covering '''8.9 mm on the target in a direction 
^P^^^brately perpendicular to the spectral orders. The sla.'/ scan axis consists 
of 384 lines ooveriiig '''13.0 mm. V7ith our 12-bit ADC, a scan rate of 20 psec/pixel 
IS appropriate so that this arrangement therefore produces a raster scan over 
'^116 mm of target area witJi a resolution of 98,304 pixels in '''2.1 sec (includ- 
ing retrace time) . Provision is made for the PDP-11/45 computer to read the 
deflection scalers so that it can check the read beam position at any time. 

The read beam may also be prograirrrably blanked or unblanked at the beginning 
of any frame. 

The video output signal is amplified, digitized and input to the PDP-11/45 
via a CAMAC input register. To enable real-time monitoring of the system, the 
video is both displayed line by line on a CKT and routed into a Hughes scan 
converter for regular TV display. 

As a figure of merit for the system, we have found that the^ limiting factor in 
^stecting v/eak spectral lines with the detector as it- new stands is the presence 
of electrical noise in the video signal leading to unit signal-to— noise ratio 
for an ''■5 x10"°Iq green-line intensity witli an H2 sec integration time. This 
figure should be capable of imuch improvement, since we do not even approach the 
detector preamplifier noise specification limits. 

ttitil the above detector system was reaefy , we configured a preliminary system 
around an SSR, Inc. Optical Multiciiannel Analyzer. This system was similar to 
the final one, differing principally in tlie target read out format and rate. 

Vhth the Or-A the tar<^t is divided electronically into 500, ''25y vbde channels 
and upon readout, Vr^iich is at the rate of 32.8 msec/ffaire, the signal in each 
^-^bannel^is added together. Since this device provides only a one-dimensional 
array Oi channels for data taking, wiiile the spectremeter oresents the scectrum 
in a ti-»-diirensional format, the camera tube was oriented v/ith the diannels 
along the grating dispersion direction and order separation acconplished 
by physically marking the target to acoorrmodate the desired set of spectral 
features. Tnin copper masks were fabricated by photo-etching holes in the 
appropriate places using standard electronic circuit board lay-out techniques, 
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and TOunted up against the vidioon faceplate. The OMA console unit contains 
axi tne neoessa^ po.ver, scanning, leakage and background subtraction, 
arrplifying, digitizing and iremory electronics to operate the vidioon. Various 
features {e.g. , frame blanking) v/ere prograrmvably controlled. Provision \iras 
made for local cooling of the VD.dioon target by '''10-15° C by enabling the 
nounting of an ann ular themoelectric unit on top of the mask; hoisflever^ this 
was not done during observations since the specially ordered thermoelectric 
module was not available until after our final detector was ready. 

There ^ several naj or disadvantages of this system when compared to the bv-o- 
dimsnsional one: (1) different observing programs require ch^ging physical 

masks; (2) spectral lines of interest falling on the same channels in different 
orders get sumred over when reading out the vidioon; (3) adding in the masked 
out portion of the target signal during channel read out to the very much 
^ller unmasked portion decades severely the signal-to-noise ratio; and (4) 
^tensity ^libration conplications due to the Very slight moving of the spec- 
trum relative to the mask as the telescope changes position through the day. 

preliminary system had in practice a noise figure comoarable 
to ^at quoted above for the h-r>-dimensional version, in principle the latter 
system shool-d be at least an order of magnitude better. 

D. DATA-AOQUISmON AI® OOTROL SYSTEM 

The DACS is oonfig\ared about a Digital Equipment Corporation PDP-11/45 coiTiDuter 
and IS equipped v/ith 4 OK of core memory, a teletype, two disk cartridge drives 
^ a dual DEG tape drive. All major user defined peripherals, including a 
CRT poyit plot display, character generator, light pen and Printer 100 char/ 

^c printer are interfaced via a CAI'FVC standard datavray controlled by a Borer 
type 1533 crate coitroller. The oonputer configuration is sham in Ficjure 8. 

^ displayed sdiematically in Figure 9, all electro-mechanical instrurrmt 
factions are controlled by the PDP-11/45 via CAMAC interfaces to stepping 
and s^chronous motor drivers. Botli detector systems were interfaced* to the 
ccnputer via Oy-IAC modules and the arrangement for the b'.X)-dirrensional vidioon 
mit is given in Figure 7, (In these Figures; I(o)R = input (output) register; 
- prograimed pulse generator; and 4BD = fourfold busy-done.) 

An enormous amount of softoare was varitten to inplemant the various system 
aspects. To effect the functions of calibration and telescope raster scan- 
lyng, a large cissortment of CA^S\C handlers were assembled into a coherent, 
^teractive system of operational subrontines. Using these basic assembly 
routines, a v/ioe v^iety of FORTRAN sofb^are was developed for DACS purposes, 
includ^g coloration, telescope scanning, sun-centering, focus, etc., programs. 

oonplete operating systems had to be prepared for each of the to*;o 
detector packages. 

All sof b^^are was run vnder the PDP-H/45 Disk Operating System Monitor and 
^ we had originally planned, under the RSX 11-D Real-Time Executive 
Operating System. Th° latter system enables a foreground/background nulti- 
programraing capability. DEC w’as simply not able to deliver the RSX 11-D in 
twe f or in^^raOon into the system. In fact, DEC's delay of several 

in their initial delivery of the PDP-11/45 seriously undermined our 
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development timetable. In practice, \ve have had an ongoing difficulty with 
the DOS; the monitor "bombs" frequently and erratically, causing loss of 
system parameters (telesccpe position, etc.) and data; DEC has not been able 
to solve this problem. 

While the flexibility deriving frcan the standardisation and modularity inherent 
to the CA^7vC scheme has proved extremely valuable in interfacing the .variety 
of components to the PDP-11/45, we did run into a unnecessary problems 
particularly ;^^ith regard to implementing the ta'r>-diirensional detector unit, as 
outlined in Quarterly Status Report #8. In addition, the Printer-100 printer 
has had a some^diat spotty madianical history; viien it is do'/.n we are forced 
to fall bade on TTY output, greatly slo^^ing davn data acquisition (i.e., tele- 
scope raster scanning) speed as well as inconveniencing software development. 
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III. OBSERVATIONS 


/ 


Ihrou^iout the report period we have used the muLLti-channel coronal spectro- 
photcjneter system* v/ith both detector asseirf^lies* to perform a variety of 
prominenoe and coronal observations. In addition to tlie usual prcminence 
emissicn lines, the most significant lines that v.-e have cbserved are Fe+9 
6374 A, Fe+10 7892 A, Fe+13 5303 A and He 10830 A. Vflienever possU^le, we 
coordinated our observations with those of SKVLAB as relayed via VMVH. 

tost of the coronal cbservations v;ere made either with the 5303 A line alone 
or with the 5303 A and 6374 A lines sirnultaneously; the signal-to-noise ratio 
for the 7892 A line v;as generally too poor to give much more than ^ indica- 
tion of its presence. V’e were not successful in detecting tlie Fe+12 10747 A 
and 10798 A lines, even after careful searching in regions where there \<jas 
substantial 5303 A and 6374 A emission? the system sensitivity at these wave- 
lengths, ha^ever, is clearly evidenced by the relative ease with vhich the 
10830 A line is observed. 

Ihe great majority of our scientifically meaningful observations were obtained 
during the second manned mission. During the first mission there v;ere no 
particularly relevant events for us to look at in relation to the SKVL^ d>- 
servati-onal programs and our work was principally of a diagnostic and instru- 
mental shakedc/ai nature. The uns>pected delay in making operational the ne^>? 
two-dimensional vidioon detector system, coi^led with generally poor v^eather 
throughout the third irannad mission, precluded any scientificcilly significant 
output during that period. 

The next four figures depict data obtained using the CMA adapted detector sys- 
tem? similar cbservations have been performed with the too-dimensional system. 
In Figure 10, we show a typical prcsrunence spectrum and illustrate the sensi- 
tivity to He 10830 A. A RG-10 hi^-pass colored glass filter inserted in 
front of the reimaging lens to very strongly attenuate the Ha line, rendering 
its intensity corrparable to that of the 10830 A line. The relatively long 
signal averaging in tliis case resulted from a combination of using the RG-10 
filter and the desire to bring out the He 6678 A line. Note that tlie line 
shapes are instrumental . 

Figure 11a shcs'/s a spectrum that includes the three coronal lines 5303 A, 

6374. A and 7892 A, tocetiier wdtli Ha. Hie gross pedastral structure is due to 
scattered light delimited by the target mask slots and some\-iiat' obscures the 
6374 A line which falls on a slot edge in this airrangement. In Figure 11b, 
we shew the 5303 A and 6374 A signals \-/ith a different mask/grating position- 
ing. Here orders 5 and 6 aire both allowed to pass through the target mask 
over the spectral region indicated. (VJiile this circumstance of overlapping 
orders reduces the signal-to-noise ratio in the derived line intensities, 
the close lateral proximity of the two lines made it necessary to record 


44 



sinu^tanc 20 us line intensities in this nvanner. IJiis situation is ccyrpletely 
opiated by tlie tov'o-dimensional detector system.) Signals are shoivTi for two 
different coronal positions. Again, the line sh^s are instroirental. 

Figures 12 and 13 sho;-/ the results of a tv-ro-dimensional coronal scan simulta- 
neous in 5303 A and 6374 A. In Figure 12, the data are plotted as a (>,0 
raster v.itJi tlie numbers denoting the eibsolute emission line intensities, 
c^ibrated via the measured intensity at the ceiiter of the sun's disk, and 
witli isophote ooiitours dravm in? in Figure 13, the same data are redrai-?n as 
curves of line intensity versus 0 at each of the various values of p, 

(Note: p = R/Rp = distance from sun center in solar radii, and 0 = geocentric 

position angle m^as^^ed counterclodwise from north.) As described above, 
tte spa^al resolution, determined by the spectrometer entrance slit, vras 
5 X 20 ' . The ^ telescope is mounted on the spar so that the orientation of 
this slit has its long side tangent to the solar’ limb at e = 135®; the orienta- 
tion is independent of telescope pointing. i 



IV. CALOJIATIONS 


A principal function of the multi-channel coronal spectrcphotorreter was to 
observe forbidden line intensities frcm coronal ions. Scare ions especially, 
such as Fe+12 v;ith three suclr observable transitions (XX3388, 10747, and 
10790), are particularly suitable for coronal region diagnostics, provided 
the line intensity ratios can be reliably related back to the ionic excita- 
tion distribution. The various proton collisional excitation rate constants 
are significant elenents in this reduction and since many of them had not 
been calculated xvith sufficient detail previously, we embarked upon a program 
of doing so for the most inportant cases. Three papers reporting on this" 
work — t^^/o describing proton excitation calculations in Fe+l2 and Fe+13, and 
the third applying these results to detemining the Fe+12 excitation equilibriurc- 
were publislied in Solar Physics and are reproduced in the Appendix to tiiis 
rqxjrt. 



FIGURE CAPTIONS 


1. Maiti-diannel coronal spectrcphotoreter rciounted on the Mees Solar 
Observatory 3.6-m spar. 

2. Data-aoquisition and control system. 

3. Overall mechanical viev/ of the instrument. 

4. Scheiretic dravrLng of the optical system. 

5. ScheitBtic drawing of the spectrometer. 

6 . T^vo-dirrensional silicon-vidicon detector asseitbly. 

7. Block diagram of the -bwo-dimensionai detector unit and sun-centering 
photo-diode electronics and system integration. 

8 . Block diagram of the oorrputer configuration. 

9. Block diagram of the electromeclianical instrument control electronics 
and system integration. 

10. Data of 24 June 1973, 2:05 tff; p = 1.04, e =^262°; 3000 frame x 32.8 

msec/frame; RG-10 filter; Aperture: 5" x 20". 

11. (a) Spectrum showing the three coronal lines 5303 A, 6374 A and 7892 A, 

together with Ha. Data of 5 September 1973, 21:34 uT; p =,1.04, 

0 = 296*^; 1000 frame x 32.8 msec/frame; Aperture: 5" x 20". 

ft)) Spectrum shoving the 5303 A and 6374 A lines with the target slot/ 
spectrum configuration used for simultaneous scan observations. 

Data of 11 Septeirb^ 1973, ''-21:25 uT; 500 frame x 32.8 msec/frame; 
Aperture: 5" x 20". 

12. Simoltaneous 5303 A and 6374 A emission line scan. Data of 18-19 
Septeriber 1974, 23:19-00:01 uT; the nunbers denote absolute line 
intensity in units of erg an“^ sec“l ster~l; Ala^s = jd5%, AI^^l - 
jd5 erg cm"^ sec~l ster“i; = 956.7". 

*. 

13. Sane data as for Figure 12. 
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Figure 3. 

Overall niechaniccLL via^? of the instrument. 
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Figure 4. 

Schemtic drawing of the optical system. 
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Figure 5. 

Schematic drai-djig of the spectrometer. 
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Block diagram of the two-diinensional detector 
electronics and system integration. 


tinit and sun-centering photo-diode 
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Block diagrazTi of the cxix^puter configuration^ 
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Fig-ore 11. 

(a) Spectrum shewing the three ccrannl lines 5303 A, 6374 A ejnd 7892 A, togetlier with 

Ha. Data of 5 September 1973, 21:34 uT; p = 1.04, 0 = 296“; 1000 frame x 32.8 
msec/frarc; Aperture: 5" x 20”. 

(b) Spectrum shaving tlie 5303 A and 6374 A lines with the target slot/spcctrim configu- 

ration used for simultaneous scan obscrv'ations . Uam of 11 Scptomlxir 1973, 0-21:25 
uT; 500 frair: x 32.8 m-sec/frame; Aperture: 5" x 20". 
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contract are reprinted in this Section. They include, in order: 


• 1. D. A. I^nclran, F. Q. Orrcill, and R. Zare, "Astronomical Observations 
witlT Television-'I^pe Sensors," University of British Columbia (1973), 
p. 255. 

2. p. A. Lfndrran, F. Q. Orrall, £ind R. Zane, Solar Physics Division meet- 
ing of tl::3 /'rcrican AstroJicmical Society, 9-11 Jan. 1974, paper P4 
(rJastract only) . 

3. D. A. Landman, Solar Physics 30 , 371 (1973) . 

4. G. D. Finn and D. A. Landman, Solar Pliysics 30, 381 (1973). 

5. D. A. Landrren, So3.gr Physics 31, 81 (1973) . 
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A. SUMMARY 

The basic objective of the contract was to create a system capable 
of producing maps of the magnetic field straight from spectra. The theory of 
the extraction of magnetic field information is contained in Section I in 
Part B on measurement of the magnetic field by fourier transform techniques. 
Part C contains contour maps of a high gradient magnetic field region. 

Section II is an operators manual, program description, and the Fortran 
coding for the implementation of the measurement procedures. Section III is 
an operators manual for microdensitoraeter which also contains the machine 
coding for the control computer. 

Data on the magnetic field was taken at Kitt Peak during the first 

and second manned Sky lab missions. Due to construction of a new solar 

magnetic field measuring facility at Kitt Peak maintenance of the main Kitt 

Peak Solar telescope suffered. Unfortunately during the period of reduced 

maintenance the main heliostat drive developed severe shake problems in light 

(■ 

to moderate winds. The amplitude of the oscillation of the image was often 
50 arc seconds and could exceed an arc minute. Because the oscillations were 
an appreciable fraction of a sun spot diameter there was little point to 
observing during the third manned mission. 

In addition to the data taken at Kitt Peak, Lockheed Observatory 
operated its multislit spectrograph during the Skylab mission. On 5 September 
1973 at l8:00 UT there was a class IB flare that was well covered by the 
Skylab instruments as well as the Lockheed multislit system. Because of the 
basic similarity of the multislit and the spectroheliograph data it was 
reasonably straightforward to produce densitometer tracings of the flare 
region. A description of the flare reduction procedure is contained in 
Part D of Section I. 

The densitometer traces of the flare represent the first ground 
based record of the hydrogen alpha spectrum with such compicLe temporal and 
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spatial resolution. In order to demonstrate and properly illustrate the 
spatial and temperal development of the Ha flare a short movie has been 
produced. The movie includes the original multislit data and the tracings 
of the flare region. Part E of Section I is a description of the movie. 

In addition to the solar useage of the densitometer we have actively 
encouraged use Of the system by outside users and have made the operating 
system of the densitometer generally available. The manufacturer of the 
densitometer now distributes program TRACE as contained in Section III as 
the normal operating system for desk based PDP-11 systems. Two systems, 
one at JPL and the other at Nice Observatory in France, are in operation. 

The National Science Foundation has funded the University of California, 
Berkeley, to essentially duplicate the Lockheed system. The Naval Observa- 
tory Station at Flagstaff is also in the process of acquiring a duplicate of 
the Lockheed system. 

Students and staff of the Universities of California and Southern 
California and the Center for Astrophysics have used the system. One of the 
users was able to trace an image tube echelle spectra. The automatic tracing 
of the echelle plate would have been impossible on any other densitometer. 

The special curve following routines that adjoin program TRACE 3 which 
make echelle tracing possible will be published shortly. When the routines 
are in final form they will be supplied to the TRACE user's group. 
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B. MEASUREMENT OF MAGNETIC FIELDS BY FOURIER TRANSFORM TECHNIQUES 

1. INTRODUCTION 

Although there are currently a number of instruments that measure 
the longitudinal component of the solar magnetic field, there are very few 
systems that attempt to measure the vector field. The measurement of the 
vector field can be especially difficult in and around sunspots. We shall 
present in this paper a method for the measurement of the vector field in 
high field regions. The method is based on the fourier transform, proper- 
ties of circularly and linearly polarized spectra arising from simple 
Zeeman triplets. 

Beckers has briefly noted that Fourier spectroscopy and the 
resulting fourier transformed profiles are useful for directly determining 
field properites. In the discussions below we will amplify on the advan- 
tages and some of the shortcomings of the analysis of the transformed 
profiles. It will be demonstrated that when Zeeman splitting is on the 
order of half of the full width at half maximum (FWHM) of the basic line 
profile, then the magnitude of the field, the inclination of the field, 
the azimuthal angle of the field, and the velocity can be determined in- 
dependent of the shape of the line profile, if the line profile obeys 
certain reasonable assumptions. In the sections on analysis of the errors 
in the fourier transform method, what constitutes a reasonable line profile 
will be discussed. 

The methodology presented here has been successfully implemented. 

The basic data are pairs of spectra in right and left circular polarized 

light and three pairs of orthogonal linear polarizations. The data 

acquisition system is called a spectra-spectroheliograph and has been 

2 

discussed in some detail by Title and Andelin . Data are microdcnsito- 
metered and digitized using a filter densitometering system described by 
Schoolman . The digitized data are organized, reduced to absolute inten- 
sities, and analyzed with a set of programs described by November^, In this 
paper we shall not present actual magnetograms. The magnetograms will be 
discussed in a series of papers on the vector field in sunspots. 
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2 . 


FOURIER TRANSFORM PROPERTIES 


Circular Polarizations 


For a normal absorption Zeeman triple in unpolarized light, the 
profile has the form: 

p = M(X) - A’(f(X-X^+ A) + h(X-X^-A)) - B g(X-X^), (l) 

where is the central wavelength of the line in the absence of a field; 
M(\) is the continuum intensity; f,h and g are the individual profiles 
of the displaced and undisplaced components; A is the Zeeman splitting; 
and A' and B are parameters that depend on the angle of the field and 
the properties of the atmosphere. 


In the region of a line, ve shall assume a constant continuum in- 
tensity. Then, for convenience, we drop the constant and treat the profile 
as an emission profile. Further, we shall assume: 

f(X-X^) = g(X-X^)= h(X-X^), (2) 

and that the profiles are symmetric: 


f(x-X ) * f(X -X). 
o o 


(3) 


Using assumptions (2) and (3) and dropping the continuum intensity, the 

Zeeman profile observed in unpolarized and right and left circularly 
polarized light can be written: 


where 


p(X-X^) = A’[f(X-X^+ a) + i)] + B f(\-X,^) 

= (o) ‘■(’“V -"(a) ® 

A' = A + C. 


( 4 ) 

(5) 


68 



The fourier transforms of eq.uations ( 4 ) and (5) are; 

p(t) = [ 2 (A+C) cos At + B] f (t) (6) 

^/rCp\^^^ " C(A+C)cosAt + B] f (t) (+)i(A-C) sin A + f (t), (7) 

\lcp) 

where : 

f(t)=J' f(u) du. (8) 

-CO 

The tilde indicates a fourier transformed function. A catalog of circu- 
larly polarized profiles and the real and Imaginary parts of their fourier 
transform is contained in Appendix I. 

If the fourier transform of the basic profile (equation 8) does not have 
any zero crossings^ and the field is not purely transverse (i.e., Aj^C), then the 
first zero of the imaginary part of the fourier transform of either circu- 
larly polarized profile yields the separation of the Zeeman components. 

That is, at: 


and 


At^ = jr, 
sin At^ = 0 
A = 


( 9 ) 


Since delta is directly proportional to the magnitude of the field strength, 
the zero crossing of the imaginary part directly determines the field strength. 
The fourier transforms of common profile functions, Lorentzian, Gaussian and 
Voight profiles, have no zero crossings. 

If (A+C) is greater in magnitude than B, then the real part of 
equation (7) or equation (6) can be used to obtain the field strength. 

Since (A+C) is greater than B, the first zero occurs when: 


At^ = Jt/2 + e, 


( 10 ) 
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and the second zero when; 



where £ is smaller than rr/2. 


Hence, 


L = 


git 

(tg + t^) 


( 11 ) 


( 12 ) 


In addition, the fourier transform of the circular polarized 
profiles readily yields the inclination of the field. The slope at the 
origin of the imaginary part of either circular polarization transform 
(equation 7) is; 


- (A-C) A f (o), (13) 

while the value of the real part of the transform at the origin is; 

= (A+B+C) f (o). (14) 


The ratio of 0^ to is then independent of the profile shape; 

V^o = (A+S+C)" (15) 

Since delta is known from the value of the zero crossing, a quanity SS may 

be defined such that 

SS = CjC A 
1 0 


(A-C) 

- . . ■ a 

' (A+Bi<;) 


( 16 ) 


The value of At^/g is rt/2, hence the cosine at At^/g is zero while 
the sine is unity. Therefore, the ratio of the real and imaginary parts 
of equation (7) yields; 
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Again, using the assumption that (A+C) is greater than B, the first zero of 
the real part of the transform occurs -vriien: 

cosCAtg) = . (18) 

Combining (I 7 ) and (I 8 ), we obtain: 



For most models of line formation, the quantity is directly related to 
the cosine of the angle of inclination. has the additional advantage 
that it is independent of the central component. 


When the Scares' relations are valid. 


A = 1/4 (1 ± cos y)^ 

(20) 

B = 1/2 sin^ Y 

(21) 

C = 1/4 (1 ± cos y)^, 

(22) 


where y is the Inclination of the field to the line of sight. 

Using Seeires' relation, result (I 6 ) is just: 

SS = cos Y* (23) 

Besides yielding the magnitude and inclination of the field, the sum of the 

right and left circular polarizations determine the Doppler shift. The 

fourier transform of the sum (the unpolarized Zeeman profile), is a 

symmetric function about the central wavelength of the profile. Hence 

its transform is real. Therefore, the fourier transform with respect to 

any other wavelength X must be of the form; 

s 


Pg (t) = p (t). 

(24) 

where 


s = X * 

s 0 

( 25 ) 
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Since p(t) Is real, the inverse tangent of the ratio of the imaginary and 
real parts of Pg("t) yields the offset, s; 

s = tan‘^ [Imag (p^Ct]) / Real(p^[t]) }. ( 26 ) 

The velocity shift with respect to the undisplaced wavelength k is; 

o 

v=scA^, (27) 

where c is the velocity of light. 

Linear Polarizations 


Using the same assumptions as used in equation (^) for circularly 
polarized light, Zeeman profiles observed through a linear polarizer at 
angle phi with the projection of B in the azimuthal plane has the form 

% \ ‘‘I * ( 28 ) 

where and are functions of gamma and phi. The fourier transform 
of equation (28) is: 

P (t) = [2E cos At + F ] ?(t). 

<p cp cp ' 

In the case, that Unno’ s relations” hold: 

^cp “ (3(A(1 + cos^ Y “ sln^ Y sin 2tp) 

^ cp = 0 / 2(1 + sin 2cp) sin^ Y ) 

where cc and ^ are functions of the absorption parameters that describe the 
line and of the angle Y. For weak lines, a and 0 are unity- 


(29) 

(30) 

(31) 
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It is useful to analyze the sura and difference of orthogonal pairs 
of linearly polarized profiles. Using relations ( 3o) and (3l), the sum 
and difference transforms are; 

Pg jp (t) = [a(l + cos^ y) cos ^t + P sin^ y] ^ 

^ 2 ^ 2 

Pq jp (t) = [-a(sin Y sin 2:p) cos At + p sin ^ sin 2 <p] 

f(€). 

A catalog of linear polarization profiles and their sum and difference 
fourier transforms is contained in Appendix II. 

If P is less than or equal to ot, the sura transform will be 
zero for the values of t symmetric about At = n. That is, at; 

At^ = n - e , 
and At^ = « + €» 
the sum transform is zero. Hence; 

i = jt/(t3 + t^^). (34) 

Then, given the value of A, the sum and difference transforms can be 
compared at t^ such that Atg is it/2. Then; 


(32) 

(33) 


% cp j = sin 2 
cp ht^=-n/2 


9 


(35) 


Hence, a pair of linear polarizations can yield the azimuthal angle of the 
field. Because of the possibility that the field is 0 or 90° (parallel or 
perpendicular) to the analyzer, it is useful to analyze several orthogonal 
pairs. Also, since the sum profile is a symmetric function, the velocity 
shift can also be obtained with a pair of linear polarizations. 
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3 . 


ACCWftCY AM) TiTKIIATIONS OF T HE FOURIER TRANSFORM MRTHOn 
Total Field Strength 


From the discussion above, the magnitude of the magnetic field, 
the angle to the line of sight, the azimuthal projection, and the line 
of sight velocity ar*e readily available from the fourier transforms of 
the line profiles. However, as with any measurement method, there are pro 
blems that occur because of both random and systematic errors. The 
fourier transform method is remarkably insensitive to some classes of 
error or noise and sensitive to others. 


The measurement of and hence the total field strength, |b|, is 
straight foreward. It depends only on the zero crossing of the imaginary part 
of the circular polarization transform. The accuracy of the zero crossing 
technique increases as the magnetic splitting increases with respect to the 
width of the basic profile. The fundamental reason for the improvement in 
accuracy with large splitting to width ratio is that the larger the splitting, 
the lower the spatial frequency at which the zero crossing occurs, while 
the narrower the profile, the higher in spatial frequency its transform has 
significant amplitude. Because of its fundamental importance, the ratio of 
splitting to full width at half maximum (fwHM) shall be defined as: 


Q = A/PWHM. 

^ 30 j 

The amplitude of the imaginary part of a circular polarization transform is 
proportional to cosine gamma. Thus, the zero crossing determination will 
also depend on cosine gamma. 

In order to get some measure of the effectiveness of the trans- 
form technique, a program was written to create artificial profiles that 
could then be subjected to various sytematic and random effects that 


74 



simulate some solar and measurement problems. To check the validity of 
the zero crossing procedure, Zeeman profiles were constructed using the 
Sears relations and Gaussian and Lorenzian profiles. The profiles so 
constructed were subjected to random noise proportional to the maximum 
signal value of the undisturbed profile. These synthetic profiles 
Were then transformed and analyzed for the first zero crossing, re- 
peated evaluation of the same profile subjected to random noise, it was 
possible to calculate the standard deviaH on of the zero crossing versus 
percentage noise in the profile. The procedure was carried out for a 
variety of gammas and splitting-to-width ratios for both Gaussian and 
Lorentzian line shapes. 


Upon completion of the analysis, it was found that, if the error 
in the zero crossing was normalized by cosine gamma, the percent error 
of the zero crossing was a function only of the ratio of the splitting 
to the width of the line profile, so that for a noise n(P„j^) in the 
profile: 




n(P ) 
max 


max 


N cos y. 


(37) 


where Wg is the noise reduction and n(B) is the standard deviation in the 
measurement of B. Shown in Figure 1 is the noise reduction factor, the ratio 
of the noise in the profile to the noise in the zero crossing, versus Q,, the 
ratio of the splitting to the full width at half maximum (fwHM) for Gaussian and 
Lorentzian profiles. From the figure it is seen that when the field is vertical 
at a magnetic splitting equal to FWHM, the noise in the measured field is reduced 
by a factor of 12 (Gaussian) from the noise in the profile. For a 60° 
inclination of the field, the reduction factor is 6 and for 75 °, it is 3 . 


In order to measure [ B[ , the position of the unpolarized line 
center must be known accurately because an error in the center position 
will be reflected as an error measured in the splitting. For a pure 
longitudinal field, an error in the line center position will cause an 
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error of equal magnitude In the splitting for a single circular polariza- 
tion. However, the centering error will cause an equal hut opposite error 
in the splitting inferred from the opposite circular polarization. Thus, 
at least for longitudinal field, the mean splitting obtained from the sum of 
the right and left circular profiles will have the correct value. 

For other than pure longitudinal fields, the centering error is 
somewhat more complex. From equation (24) in the presence of a centering 
error, s, and a splitting. A, the condition for the zero of the imaginaiy 
parts of the transform is: 

A sint(s + a) + C cint(s - A) + B sints « 0. (38) 

Since, in the absence of a centering error, the zero of the imaginary 
part occurs when: 


At^. = It , 

in the presence of a small error condition (58) occurs when 

At^ = jt + Ct, (39) 

where Q is just the splitting error. 


Substituting (39) into (38) yields the relation between the 
centering and splitting error: 


A sint(s + C) + C sint(s - C). " B sints «= 0. 


(40) 


If the errors are small, then; 


^ (A + C - B) 


(W) 
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Using the Seares relations, equation (ho) becomes: 

C = ±(cos y)s. (42) 

From results(4l) and (42), the splitting error averages to zero when the 
splittings obtained from the right. and left circular polarizations are 
averaged. Further when Seares relations hold, the error in |b| caused 
by a centering error is diminished by the cosine of the inclination. 


One of the assumptions of the fourier transform technique is that 
all three of the Zeeman components have the same profile shape. However, 
if the undisplaced profile is symmetric it may differ from the shape 
*^f the displaced components without affectxng the value of the zero 

crossing because the imaginary part of the circular transform 
is free of all profile components that are symmetric about the profile 
center. The lack of dependence on the central component is very useful 
because it means scattered photospheric light or the existence of molecu- 
IS'^ lines centered on the profile do not affect the value of |b| . 

Another transform assumption is that the profile shape is symmetric 
about the undisplaced center. There are at least two physical conditions which 
can cause line profile asymmetry - magnetic field gradients and velocity 
field gradients . Magnetic field gradients cause mirror asymmetry in the 
displaced profiles. That is, the profile displaced to high wavelengths is 
the mirror imge of the component shifted to shorter wavelengths. In the 
case of the mirror asymmetry , the profile can be considered to be made of 
a sum of profiles that are shifted by differing amounts. Since the 
fourier transform procedure is a linear process, the zero crossing will 
reflect a weighted average magnetic splitting. In the case of a velocity 
gradient asymmetry, all three profiles are asymmetric in the same direc- 
tion. The first effect of a velocity gradient will be an error in the 
center wavelength of the sum line profile. As discussed above, if 
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simultaneous profiles are obtained in right and left circular polariza- 
tions, the velocity error will increase the field estimated from one pro- 
file and decrease the field from the other. The order of the error will 
be the same, so that the average of the right and left circular polariza- 
tion fields will be a good estimate of [b] and the difference will be a 
measure of the velocity gradient error. 

Even if the line profiles are symmetric and well centered, the 
fourier transform can still yield erroneous results for |b| if there 
exist photometry errors. To get an idea of the magnitude of this error, 
profiles of the form: 

pj(X) a 

were analyzed. It is clear that for pure longitudinal fields that 
the value of the zero crossing is unaffected by 6 . However, as the field 
inclination increases, the effect of non-zero 6 on the zero crossing increases. 
On the other hand, as Q, increases, the effect of non zero 6 should decrease, 
since the profile components overlap region decreases. Shown in Figure 2 
are plots of percent error in the zero crossing versus inclination of 
field for 6 = -.05 and -.1 for Q = . 5 . Shown in Figure 3 are plots of 
percent error in the zero crossing verus Q for the same values of 6 for 
a field inclination of 8 ^ . Negative values of 6 cause an increase 
in the splitting values while positive values cause a decrease of splitt- 
ing, For values of 6 < 1«08| the magnitude of the crossing error is 
nearly independent of sign of 6 . Note that the values of Q and Y used in 
Figures 2 and 5 respectively were chosen to illustrate maximal error sensitivity. 

In practice, it should be possible to correct the photometry so 
that the error in 6 is less than |. 05 l. 

THE INCLINATIOM OF THF! FTFT.D 


Once the value of |b| is known, the inclination can be determined 
as indicated by the series of results (I 7 ), (I 8 ) and (I 9 ). Result (I 9 ) 
yields cosine y independent of the value of the central component and 
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only requires Unno' s relations to hold. However, since result (19) is 
independent of the central component, it cannot be expected to be of much 
value for small values of the ratio A/C, For inclinations of less than 
^5 , A/C < .03. But by 60 , A/C = .11. Thus, it should not be surprising 
if the accuracy of cosine y markedly deteriorates at inclinations of less 
than 45°. 

For angles of inclination less than , the slope at the origin of the 
imaginary part of the transform can be utilized. However, not only does this 
method require that Scares relation hold, but perhaps more importantly, intro- 
duces dependence on the strength of the possibly contaminated central component. 
Another method of obtaining cosine Y is to take advantage of the properties of 
linear polarization transforms, since these depend on sine gamma. 

The method used to obtain result (I9) can be improved upon. Since 
the value of A is known, the ratio of the real to the imaginary part of 
the transform can be multiplied by sin A t to create a function; 


RS(t) = 


(A+C) costA + B 

Ta^^c) 


(44) 


The integral of the product of RS(At) and the first two Legendre polynomi- 
als properly normalized over the range in At for 0 to tt yield 

J RS(At) d(cosAt) (^5) 

J* RS(t) d(cosAt) 

Since the Lengendre polynomials are orthogonal functions, the integrals of 
RS(At) times higher order polynomials yield information on the degree of 
asymmetry and or differences between the central and displaced components in 
the line profiles. 


The value of the integral approach of equation (46) is that a 
significant portion of the transformed function is used rather than the 
single point which is used for result (I9). However, the formation of 
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RS(zit) entails dividing out a denominator that can take on zero values. By 
use of Gaussian quadrature, the points at which RS(At) can become large can be 
avoided without loss of numerical accuracy. When the correlation coefficient 
of the error in |b| and the error in cosine Y are evaluated from profiles sub- 
jected to random noise, for inclination greater than h<f and Q > . 5 , the 
correlation coefficient is less than . 01 . 


Using the random synthetic profile program, the noise in cosine 
gamma was evaluated. As with the error in 1 b|, the noise in cosine y 

is reduced by a factor normalized by cosine gamm. Thus, for an amount 

of noise n(P ); 

max' 


and 


N V) ^ 

Y cos y 


n(P ) 
max' 


max 


cos y, 


n(cos y) 


n(P 


max 


) 


p 

max 



(47) 

m 


Thus, the error in cosine gamma Is independent of gamma. Shown in Figure 

4 is a plot of versus Q. Relations (47) and ( 48 ) hold for angles greater 

than 45 . For angles less than 45°, methods for determining cosine y 
from result ( 46 ) fail. 


Since it is not until slightly under 15° that the ratio: 

a/b < . 03 , 

it can be expected that result ( 2 J>) is useful between 15*^ and 45° for 

determining cosine y. Numerical experiment In fact shous that the slope 

of the iraglnary part of the transform Is useful for Inclination less 
than 45°. 


The polynomial measurement of cosine gamma is only weakly dependent 
on differences in shape between the central and displaced components, and 
errors introduced by such differences ore manifested in the coefficients of the 
IS er order Legendre polynomials in the expansion of BSCdt). Also, as with 
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|b| errors in centering the line cause equal and opposite effects on the 
value of cosine gamina measured* As a consequence^ the average value of 
cosine gamma will be a good estimate even in the presence of a centering 
error. 


Photometry errors cause a somewhat different effect on cosine 
gamma than on the measurement of |b|. As shown in Figure 5^ the percent 
error in cosine gamma versus inclination of field only slowly increases 
with angle, while as shown in Figure 6, the error in cosine gamma versus 
Q. does not decrease with Q as does the error in | b| , but rather increases. 
Further, as seen from Figures 5 Q^nd 6, the size of the error in cosine Y 
is only a slowly varying function of both gamma and Q. Shown in Figure 7 
is a plot of the error in gamma versus gamma for a ten percent error in 
cosine Y. Figure 7 d,eomonstrates that estimates of the inclination of 
larger angles are almost certainly correct while the accuracy of angles 
less than 40*^ is extremely sensitive to photometry. 
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5. DISCUSSION 

From the arguments above there are two major advantages of the 
fourier transform method. The first is that the magnitude of the field 
can be measured independent of the shape of the profile. The second is 
that for large fields, the magnitude of the field can be measured with an 
accuracy which is a factor of 10 to 25 better than that of the basic line 
profile. Since measurement accuracy depends on the square of the number of 
photons counted, the factor of 10 to 25 increase in accuracy represents a 
factor 100 to 625 in observing time required - 


The main disadvantage to the method is a possible systematic photo- 
metry error. At present, we can measure the field in sunspots relatively 
to an accuracy of about one percent, with a possible systematic error of 
5 percent. Cosine Y can be measured relatively to 3 to 4 percent and 
absolutely to 8 to 10 percent* 

The fourier transform method does not require a great deal of com- 
puter time. This is not because a fast fourier transform method is 
used, but rather because a relatively few fourier components need be 
evaluated to find the zero crossing. Usually less than ten fourier com- 
ponents must be evaluated. For Just a few transforms it is sufficient 
to determine the sines and cosines recursively. 
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FIGURE CAPTIONS 


Figure 1 


Figure 2 


Figure 5 


Figure 4 


Figure 5 


Figure 6 


Figure 7 


« 


Noise reduction in zero crossing from noise in profile data 
versus the ratio of split ting /FWHM. Curve (a) is for Gaussian 
profiles while (b) is for Lorenzian profiles. 

Percent error in zero crossing versus inclination of the field 
for photometric errors of size 6, for Q s .5. 

Percent error in zero crossing versus splitting FWHM for several 
photometric errors for a field inclination of 82.82® (cos Y = . 125 ). 

Noise reduction in the measurement of cosine gamma versus 
spl i tt i ng/FWHM . 

Percent error in cosine gamma versus gamma for photometric errors 
of size 6 for Q = 1.5. 

Percent error in cosine gamma versus splitting/FWHM for photo- 
metric errors of size 6 for a field inclination of 82.82®. 

Percent error in angle versus angle for a 10 percent error in 
cosine gamma. 
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NOISE REDUCTION 





Figure 2 





Figure 5 
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APPENDIX I 


Shown in Figures A1 through a 6 are plots of Gaussian based Zeeman 
profiles in a single circular polarization and the corresponding real and 
imaginary parts of the profiles' fourier transform. The set of figures 
include profiles with inclinations to the vertical from O*’ to 75° in 15° 
steps of . 25 . Scares' relations were used for the relative amplitudes of 
the three Zeeman components. The profiles are normalized so that peak 
amplitude at zero inclination is unity. The transforms are normalized so 
that the amplitude of the real part of the transform is unity at zero 
spatial frequency. On an individual figure which covers a single field 
inclination each row contains a profile and the corresponding real and 
imaginary transforms. Successive rows on a given figure are for suc“ 
cessively higher ratios of Q, the ratio of splitting to FWHM, 


92 



Figure A1 

CIRCULAR polarizations 
FIELD INCLINATION O 

PROFILES TRANSFORMS 


REAL 


IMAGINARY 




Figure A3 


CIRCULAR POLARIZATIONS 

field inclination 30 

PROFILES transforms 


REAL 


imaginary 



ZO\ 


Figure 

CIRCULAR POLARIZATIONS 

FIELD INCLINATION 4S 


PROFIL ES 


transforms 


REAL 


IMAGINARY 


.5 







96 



Figure A5 
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Figure A6 
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APPENDIX II 


Shown in Figures AJ through A15 are plots of Gaussian Zeeman profiles 
in pairs of orthogonal linear polarizations and the fourier transforms of 
the sum and differences of the pairs. The profiles and their transforms 
are organized differently from the circular case and a smaller fraction of 
the possible angular and splitting combinations are shown. Included are 
inclinations from the vertical from 50 to 90° in 50° steps, orientations 
with respect to the azimuthal projection of the field from to 90° in 
15 steps at three splitting to FWHM ratios. Each figure represents a 
single inclination from the vertical and a single splitting to FWHM ratio. 
The first two columns contain the profile at ± 45° to the azimuthal pro- 
jection of the field. In the third column of the first row is the transform 
of the sum profile. The difference profile is zero and is not shown. 

Further since the sura profile is the same for all azimuthal orientations 
it is only shown in row one. All successive rows contain profiles in 
orthogonal polarizations at the indicated angles to the azimuthal projection 
of the field and the fourier transform of the difference profile. 
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Figure A7 
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Figure A9 
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Figure AlO 
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Figure All 
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Figure A12 
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Figure AIJ 
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Figure A15 
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C. CONTOUR MAPS OF A HIGH GRADIENT REGION 

In the first paper on the Spectra-Spectroheliograph a sequence of 
original data was shown that apparently exhibited field reversal at a field 
strength of over 2000 gauss. In the course of all our observations this 
region presented the most complex magnetic structure. It was therefore 
chosen as the final test for the data reduction procedures. 

Shown in Figures 8 a and b are contour maps of the magnetic field 
strength (jBj), and the continuum intensity for the high gradient region. 

'■ f 

Corresponding points in 8 a and b as well as 9 ^^ b^ and c, represent identical 
points on the surface of the sun. The field strength is shown in hundreds of 
gauss. The continuum intensity is in arbitrary units j however ^ lower numbers 
indicate lower intensities. The regions with intensities below 50 are 
umbral. Figures 9^^ b, and c show contour maps of the longitudinal component 
of the field fe cos v), the inclination of the field (y)j and the line of 
sight velocity (v). The longitudinal component of the field is given in 
hundreds of gauss. The inclination of the field is given in degrees from 
the vertical. The velocity is in kilometer per second, positive downward. 

Comparison of Figures 8& and 9a shows that the maximal field strength 
occurs where the longitudinal field is zero. For convenience the contour of 
zero longitudinal field is shown on the field strength map (8a) and is 
marked by zeros at the outside of the map boundaries. It is clear from ,8a, 
that there exist gradients of the order of 6 OOO gauss/arc second in the 
magnitude of field in the neighborhood of the zero longitudinal field contour. 
The maximum gradient in the longitudinal component of the field is approxi- 
mately 700 gauss/arc second. 

Comparison of Figures 9a and 9 b shows that, at least for this region, 
the contours of the longitudinal fjeld inside, the spot are dominated by the 
inclination of the field. Figure 9 c indicates that the velocity toward the 
surface apparently increases as the field inclination becomes more vertical. 
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The data reduction for the contour maps was performed completely 
automatically. The film was traced using the roster features of program 
TRACE in less than five minutes. Uie data on the computer was reorganized 
under a program called PASS2 (see Section 2) in a running time of 8 minutes. 
Then the data for the maps were created by PASS5 in slightly under 20 
minutes. 


FICSURE CAPTIONS 


Figure^ 8. Contour map of Sunspot High Gradient Region, a - Magnitude of 
the field in hundreds of gauss, b - Continuum intensity in 
arbitrary units. 

Figure 9* Contour map of Sunspot High Gradient Region, a - Longitudinal 
component of the field (b cos Y). b - Inclination of field to 
line of sight (Y). c - Line of sight velocity (positive 
downward) . 
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D. 


DATA REDUCTION PROCESS 


As an illustrative example, let us consider the reduction of the HD: 
spectra of the flare of 5 September 1975 as recorded by Lockheed’s multi - 
slit spectrograph (mSS), This event was seen and studied by Skylab. A 
sample of the MSS data is shown in Figure 10. Bie MSS differs from a conven- 
tional spectrograph in that many entrance slits are placed at the solar image . 
A 7 Angstrom blocking filter centered on HX is introduced into the beam to 
prevent the spectra from the several slits from overlapping. 

The spectra were first digitized with the PDS microdensitometer. 

(The densitometer, the associated PDP-11 computer, and the control software 
which facilitates the use of the densitometer, are described in detail else- 
where in this report.) A total of ^5 frames (taken at 15 second intervals) 
were traced. This covered the period from just before the onset of the flare 
to well into the declining phase. On each frame, 70 positions along the slit 
were traced, each separated by about one arc second. A calibrated step wedge 
was also traced. The d.ata were in the computer's disk memory. 

The processing of the numerical data proceeds in three steps. First, 
the raw digital data are converted to relative intensity values, using the 
step wedge results and the known characteristic curve of the film. Next, the 
profile of the 7 A blocking filter was divided out, thereby recreating the 
original solar spectral profile. Finally, the results are displayed and 
compared in several ways. 

One obvious display is the comparison of a profile to the pre -flare 
profile at the same location. This can be done both for successive positions 
along the slit and for successive times at a given slit location. An example 
is shown in Figure 11. The drawings are made with a computer-controlled CRT. 
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We have found that a particularly useful display is produced by 
stacking the profiles diagonally one behind the other. This gives a perspec- 
tive view of the flare spectrum along the slit. Examples from selected frames 
are shown in Figure 12. Moreover, when these perspective views are projected 
cinemagraphically, the spectral dynamics of the entire event become clear. 

This display is included in the l6 mm film which is appended to t>n.s report. 


FIGURE CAPTIONS 

Figure 10. Multislit spectrograph in ttn for the flare of 5 September I 9 V 5 . 

Figure 11> Flare profile (solid line) compared to the preflare profile 
(dotted line). 

Figure IE Perspective view of the flare spectrum along the slit. Each 
frame shows the spatial distribution. 
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THE FLARE OF 5 SEPTEMBER 1973 ~ A FIIM 


The accompanying l6 mm film shows the solar flare of 5 September 1975 
(a Skylab event) as recorded and measured by the iK^ckheed Solar Oh.-'ervatory . 
The film consists of six sequences, each covering approxim-ately tVie sati'' 
period of time, but each showing the event in a different way. Only wh»^‘n all 
of these aspects, plus others not detected at Lockheed, are taken into 
account can a coherent picture of the flare process be assembled. 

Sequence 1. This shows a direct view of the flare as seen through a filter 
centered on ir? (the Balmcr-a line of neutral hydrogen at 6563 A). The picture 
were taken by the slit monitor camera of the multi ~slit spectrograph. Tiic 
vertical black lines are the entrance slits of the spectrograph. 


Sequence 2 . liiis shows a series of HD; spectra as recorded by the multi -silt 
spectrograph. A 7 A filter, centered on W( is placed in tfie optical path. 
Thus, each entrance slit produces a segment of spectrum about 7 A wide. The 
slits are sufficiently separated that there is only a small overlap between 
adjacent segments. The vertical dark lines seen in the pictures are ea'-h 
the core of the Itx line coming from a single slit. As the sequence progn ;!ses 
the bright flare can be seen in emission (bright) in the center of H3, Thf' 
rapj-.ily moving surges appear in absorption (dark) and are Doppler-.uhifteu 

uwu,, from the core of the line. Blue (shorter wave'length) ir; to IcCt, /•fU tn 
ill!' fight. 

The flare is seen through a narrow filter centered 1.;':' A to > ho 
blu- (short Wavelength) side of the center of Rx. Thus, f<?atures which am 
risin.r from the sun and are thereby Doppler -shifted toward the blue will 
’ipp-ar prominently. 

y:-:ii uencc ,_ 4 . This is similar to Sequence 5, except tliai, the filter sits on 
III rrd .-.idc ot Hd:, so that fa].lirig features appear prominently. 
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Sequence 5 . This is a repeat of parts of sequences 1 (slit monitor) and 
2 (multi-slit spectra). The arrow indicates the slit along which the flare 
spectrum has been digitized. In the multi-slit spectra sequence, note that 
the bright flare emission occurs at two slightly separated positions along 
the slit. 

Sequence 6 . This shows, in animation, the spectral dynamics of the flare. 
The data were digitized with the PDS microdensitoraeter, scanning perpendicu- 
lar to the slit. After reduction, the individual traces along the slit were 
drawn on a CRT stacked one "behind" the other, so as to give a perspective 
view of the spectral intensity along the slit. The sequence is divided into 
two parts, each taken from one of the two bright emission features indicated 
in Sequence 5- Note that the center of the line goes into smooth emission 
in the first part, while a "self -reversal" absorption profile persists at 
line center in the second part. 
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FOREWORD 
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in support of the Apollo telescope mount was conducted for 
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NAS 8-28028. This program was performed under the technical 
direction of; 
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Code S&E - SSL-TE 

George C. Marshall Space Flight Center 
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This program was performed with the assistance of: 
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Haskell Shapiro, Engineering Consultant, 

Corona del Mar, California 
Spectra Optics, La Crescenta, California 
and others too numerous to mention without whose assistance this 
project could not have been carried out. 
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ABSTRACT 


This program has resulted in the installation of a 65 cm 
Photoheliograph Functional Verification Unit (FVU) in a major 
solar observatory. The program also has led to development 
of a 10" vacuum refractor telescope and a 10830 A birefringent 
filter for obtaining Helium line observations. The telescope 
cluster, which also includes a previously developed 8.6" 
Singer-Link vacuum refractor and a Hasselblad white light tele- 
scope, is being used in a daily program of solar observation. 
The 10" and 8.6" vacuum refractors were used to obtain high 
resolution ground based observation in support of the Apollo 
Telescope Mount. Results were facsimile transmitted directly 
from the Big Bear Solar Observatory to the Lyndon B. Johnson 
Space Center at Houston, Texas during the course of the ATM 
mission. 



I Introduction and Historical Background 


1 . Big Bear Solar Observatory (BB50) Description and 
Capabilities 

In full operation at the start of the SGAP program, Big Bear 
Solar Observatory (BBSO) , located in Big Bear Lake in 
California's San Bernardino Mountains, possessed two excellent 
10" refracting telescopes, an 8.6" Singer Link vacuum refractor 
used to obtain Skylab astronaut training data, a small Hassel- 
blad white light full disc telescope, a spectrograph fed by a 
10" telescope through Coude* mirror and a video magnetograph. 

The telescope drive system which is of the declination/right 
ascension design possessed a principal guider with offset 
control affecting all telescopes except the 8.6" Singer Link 
which possessed its own guider enabling it to always point 
at the sun's center. The BBSO dome and dome shutter are 
independently guided. BBSO possessed a variety of 1/4 A and 
1/2 A H~a filters and a 0.3 A calcium K-line filter. During 
the course of the SGAP program a Universal Birefringent Filter 
tunable for wavelengths between 4000 and 7000 A was delivered 
and its checkout begun. 

2 • 65 cm Photoheliograph Functional Verification Unit (FVU) 

Description and Status 

A 65 cm photoheliograph development program was conducted by 
JPL and Caltech under contracts sponsored by NASA headquarters 
during the time period from 1966 through 1973. The 65 cm photo- 
heliograph flight hardware, which was initially proposed for 
an ATM-B mission that never materialized, was not built but a 
Functional Verification Unit (FVU) was designed and fabricated 
complete with flight quality CER-VIT optical elements. The 
FVU is a 65 cm aperture Gregorian design with an Invar, low 
thermal coefficient of expansion, truss structure used as a 
base for mounting of the optical elements. Its 250 cm focal 
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length f/3.85, primary mirror is coupled with an ellipsoidal 
secondary mirror possessing a 13 times magnification ratio 
which yields an effective telescope focal length of 3250 cm. 

Due to design and fabrication constraints, of the 65 cm primary 
mirror aperture approximately 2.5 cm at the outer edge is 
masked off and lost due to the beveled rim and slightly turned 
down edge which remained after polishing was completed. 

Therefore, although the FVU effective local ratio was designed 
to be f/50, it is closer to f/55 in actual operation. The 
theoretical diffraction limit for a primary mirror of this 
dimension is approximately 0.2 arc seconds for high contrast 
targets. Due to problems encountered in polishing and mounting 
the primary mirror, its resolution is somewhat reduced from the 
theoretical limit since LUPI- grams obtained during tests indicate 
that it possessed a surface good only to approximately 1/5 to 
1/6 wave peak to peak. A diffraction limited primary mirror 
should possess a surface figure good to from 1/8 to 1/10 wave 
peak to peak. 

3. Installation of the 65 cm Photoheliograph FVU at BBSO 

It became clear quite early in the 65 cm FVU program that Big 
Bear possessed seeing sufficiently good that imagery nearing the 
FVU's resolution limit might be obtained frequently during the 
summer months and occasionally although' less frequently in the 
winter months. It was also apparent that the integration of the 
65 cm FVU into the telescope cluster at BBSO would represent 
bhe least costly and most timely means for obtaining operational 
experience with a telescope of this design. As it turned out, 
interference testing of the FVU in full simulation of the sun, 
relative to total energy and collimation angle, and vacuum 
conditions of space would have been more costly than mounting 
the FVU at BBSO and using it as a ground based telescope where 
meaningful solar data could be obtained as well. 
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It was therefore proposed to MSPC as part of the SGAP program 
to install the 65 cm FVU on a modified mount at BBSO, providing 
as a part of this new installation one of the existing BBSO 
10" refractors, the existing 8.6" Singer Link vacuum refractor, 
and the Hasselblad white light telescope. The new design was to 
provide a Coude' capability so that the 65 cm FVU could be used . 
to feed the spectrograph on command. In addition a Helixim 

O 

10830 A birefringent filter was proposed for development. A 
Proposal was submitted, revised and resubmitted to MSFC and a 
contract to accomplish this installation prior to the launch of 
the Skylab ATM mission was awarded to Caltech in early 1972. 

The following sections describe the work that was accomplished 
under that contract, 

II Scope of Work Accomplished 

The scope of work that was accomplished is described in the 
following sections divided on the basis of major subsystems: 

i* 65 cm Photoheliograph FVU Subsystems 

In order to use the 65 cm FVU at BBSO the following modifica- 
tions and newly developed equipment items were required. 

A. A vacuum tank of sufficient size to enclose the entire 
telescope less camera rack and provide an operating 
environment at several millimeters of mercury pressure 

or less. In association with the vacuxom tank, an aperture 
window was required possessing sufficient quality in 
transmission that the 65 cm FVU image quality would not 
be compromised . 

B. A new vacuum guide telescope was required to provide 
pointing accuracy and stability of better than one arc 
second . 
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C. An aperture shutter was required to protect the 65 cm 
PVU structural elements in the case where inadvertent 
mispointing might have caused the primary mirror's 
image of the sun to fall on portions of the secondary 
assembly not protected by the heat stop mirror. 

D. As time progressed it became obvious that an extension 
sunshade was also necessary to protect the large aperture 
window from high thermal inputs at its edge which were 
conducted from the front of the vacuum tank which was 
exposed to the direct sun. 

After competitive procurements Ball Brothers Research Corporation 
(BBRC ) ! Boulder, Colorado was selected to construct the vacuum 
tank and Muffoletto Optical Company, Baltimore, Maryland was 
awarded a contract to provide the entrance aperture window. 
Selection in both cases was based on price, familiarity with the 
existing 65 cm FVU and its interface requirements, and the 
historical reputation of both contractors for meeting specifica- 
tions. The window was fabricated of Corning BSC 52 glass possess- 
ing extremely low inhomogeneity. To aid in testing, a wedge 
was purposely ground and polished between the window surfaces 
so that light reflected from each of the surfaces could be 
caused to produce a regular interference pattern that could in 
turn be used to evaluate the transmissive properties of the 
window once the surface quality was known. As it turned out, 
pi^oducing the window was much more difficult and time consuming 
than had ever been imagined. Its late delivery eventually 
caused the entire program to slip so that data from the newly 
modified telescopes was available only during the third of the 
astronaut visits to Skylab/ATM, However, data from existing 
telescopes was collected during each of the two prior crews 
visits so no operational time was lost during the first two 
manned portions of the Skylab/ATM flight program. 
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The new guide telescope was mounted within the 42" diameter 
vacuum tank on the FVU structure to minimize effects of any 
mechanical misalignments between telescope and vacuum tank. 

This necessitated provision of a separate guide telescope 
aperture window in the vacuum tank. 

Valving and a high capacity pump had previously been procured by 
JPL. It was tested, repaired and installed at BBSO with the 
pump exhaust vented through the concrete dome walls, A 
molecular sieve was procured to prevent any backstreaming of 
contaminants in case of failure of the pump with the valve 
open. A long corrugated stainless steel flexible tube was 
procured and installed on hangers near the pump for ready 
connection to the vacuum tank when pumping is necessary. 

The vacuum system holds exceptionally well. Pumpdown is 
accomplished no more frequently than once a week or so, as 
long as the tank isn't deliberately vented for maintenance. 

Vacuum within the specified range is achieved within less than 
one-half hour after pumping starts. All vacuum system 
specifications were exceeded and the system works very 
satisfactorily . 

On the 65 cm FVU itself, both the second diagonal mirror assembly 
and instrumentation rack assembly were removed from the 
integrating structure. The second diagonal mirror assembly 
has been remounted in a hat section on the sun end of the 
42" vacuum tank where it redirects the beam to the coude ' . 

The second diagonal mirror has been provided with two axis 
motor driven tilt capability for remote alignment and eventual 
image motion compensation drive. The instrumentation rack 
was redesigned in three sections and these are mounted on the 
Underside of the rear of the 42" vacuum tank directly rearward 
of the Coude’. Flip mirror and beamsplitter assemblies have 
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been mounted to feed any of three optical benches with provision 
for other benches if desired. Typically the Zeiss Universal 
Birefringent Filter is placed on the east bench and an H-a 
filter on the west bench. The center bench has less clearance 
but has been used for a motor driven Nikon camera and could be 
used for a TV monitor. 

2 . 10” Vacuum Refractor 

Although the objective and field lenses for the 10" telescope 
already existed in the original mount, a new structure had to be 
constructed to position and hold these elements. The decision was 
made to make the 10" refractor also a vacuum telescope using the 
objective as the entrance aperture. It was also recommended that 
independent guiding and offset drive be provided in a gimballed 
mount so that the 10" vacuum refractor could be made to point 
at the opposing limb of the sun if the main mount and 65 cm 
FVU were trained on one limb of the sun. The amount of travel 
required was twice that required for the Singer Link telescope. 

In addition a sunshield and split beam optical bench were 
provided to improve image quality and utilization respectively. 
Furthermore, an achromatizing lens was introduced so that image 
size identical to that obtained in H— a could be obtained from 
the calcium K-line filter. The following equipment was provided: 

A. A vacuum tank of sufficient size to hold and position 
the objective lens, field lenses, and achromatizing 
lens complete with suitable access ports for maintenance 
and adjustments. 

B. A sunshield and dust cover for the sun end of the vacuum 
tank . 

C. An independent vacuum guide telescope to be strapped to 
the 10" refractor provided with a translating sensor unit 
to be used for offset pointing. 
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D. Micrometer /motor drive units for guiding and providing 
offset pointing capability. Worm drive has provided 
improved stability under external loads whether 
disturbances were due to wind or observer operation. 

E. A vacuum gualified motor driven field lens and adjust~ 
able field stop assembly combined on one optical bench. 

F . A manually retractable achromatizing doublet lens 
designed by Dr. Ira Bowen of The Hale Observatories 
staff to provide identical magnifications for H-a and 
Calcium K-line images., 

G. A dual optical bench with beamsplitter assembly that 
can accept spectral or neutral density beamsplitters 
for simultaneous recording by two camera systems. 

H* Valving and a vacuum gauge for use with the high 

capacity rotary pump used on the 42” vacuum chamber. 

All of the above elements wer«= built at Caltech with the exception 
of the achromatizing lens. Tl.e system was integrated without 
difficulty, holds vacuum almo.;;t indefinitely and in short has been 
remarkably trouble free. It las been operational since first 
pumpdown at BBSO with the exception of time required to install 
a flip mirror replacing the b«.:amsplitter feeding the side by 
side optical benches. The flip mirror reduces efficiency losses 
and other side effects generrvted by spectral and neutral density 
beamsplitters. The 10” vacuijm refractor is usually provided with 
an H-a filter on one bench and the Leighton Videomagnetograph 
on the straight through bench. Its guider system appears to 
somewhat damp out guiding disturbances and mispointing 
occasionally seen in the 65 cm PVU guider. 
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3* 8.6” Singer Link Vacuum Refractor 


The 8.6" Singer-Link Vacuum Refractor was originally piggy- 
backed on the first telescope installation at BBSO. It has always 
provided independent guiding and offset pointing so that it can 
constantly guide on the sun's center. Accordingly, only minor 
modifications were necessary to gimbal arid mount the Singer 
Link in position at the side of the 10" refractor and under 
the declination axis. The vidicon camera and beamsplitter 
were rotated 120° to avoid interference with the 10" telescope 
and the BBSO mounting pier. 

4. Integrating Structure 

integrating structure was reguired to mate the 42" vacuum 
tank containing the 65 cm FVU, the 10" vacuum refractor and 
the 8.6" Singer Link vacuum refractor with the existing BBSO 
pier and right ascension and declination drive assemblies 
and axis. Fabrication alignment of this integrating structure 
was critical if it were to accept the existing declination 
axis and properly align the 65 cm FVU anditsCoude' system. 

The Coude’ system is of necessity an important part of the 
integrating structure. The Coude mirrors are self aligning 
for changes in declination. They are passively driven by 
mechanical coupling to the telescope base near the existing 
spectrograph. The Coude' system, because of its location, is 
extremely difficult to gain access to. It operates in vacuum 
which makes alignment access even more difficult than it would 
be just due to structural reguirements. The system works but 
will eventually receive added modifications to improve align- 
ment accessability and to add additional alignment adjustments. 

Calculations of telescope pointing extremes and seasonal 
variation of declination coupled with telescope configuration 
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showed that the northeast and northwest surfaces of the BBSO 

srid mounting' plate assembly needed to be shaved to provide 
adeguate clearance in normal operation. A concrete sawing firm 
was contracted with to accomplish this task after the old tele- 
scope was dismantled and prior to installation of the new 
telescope cluster. Adequate clearance is now available for 
midsummer operations encompassing any of the allowable 
telescope pointing attitudes. 

5 • Hasselblad White Light Telescope and Systems Accessories 

The Hasselblad White Light Telescope which is mounted on two 
separate optical benches was mated to the 42" vacuum tank 
containing the 65 cm FVU. A broadband 500 A bandpass green 
■^^ker f erence coating was applied to the objective lens in hopes 
that a greater degree of protection would be afforded the 
various neutral density and spectral filters placed in the 
beam near the camera. The shutter was retained. Even with 
the interference coating, experience has shown that the heat 
load on these thin filters is too great if left in the beam for 
very long. As long as the shutter is used for protection, no 
breakage occurs. Therefore, while the interference coating 
provides additional protection for everything in the image field 
of the objective, use of the objective shutter is still required 
between exposures to prevent inadvertent damage. 


Optical finders including a rifle sight have been added to the 
42" vacuum tank and the 10" telescope to aid in initial pointing 
and boresighting. 

Racks have been provided on the 42". vacuum tank for mounting 
lead weights necessary for telescope balance when instrumentation 
is moved or replaced by instrumentation of different mass or eg. 

O 

6* 10830 A Birefringent Filter 

o 

A 10830 A helium line birefringent filter was procured under 
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terms of this contract. The filter was originally specified to 

O O 

provide half height bandwidth at 10830 A not exceeding 0.6 A. 

O 

It was specified to be tunable over a spectral range of 20 A 

o . 

centered at 10830 A . Spectral nonun if ormity at any point in 

O 

the 30min clear aperture was not to exceed + 0.1 A. Each filter 
cell was to be capable of independent adjustment so that the 
filter might be fine tuned at its equilibrium temperature. The 
filter was to reject all side bands over the actinic range of 
the silicon matrix television camera tube. All air to glass 

O 

surfaces were to be anti reflected for 10830 A and exhibit 
1/10 wave flatness. All internal optical surfaces were to be 
oil contacted. The filter was to incorporate a contrast element 
whose inner order separation is 5/7 the 32n element giving a 
bandpass of 7/5 times that of the smallest element. The filter 
was to utilize two concentric' heating shells, the outer shell to 
provide the inner shell with a stable temperature hold within 
which to operate. The outer shell was to be thermostated at 
40‘*c + 0.03**c. The filter was to operate on 24 volts D-C. 

A revision later changed the half height bandwidth specifications 

o o 

at 10830 A to not exceeding 0.3A . 

As a result of competitive procurement Douglas Martin doing 
business as Spectra Optics of La Crescenta, California was 
awarded a contract to construct this filter on April 21, 1972. 

On March 5, 1973 Spectra Optics delivered the completed filter 
for acceptance testing. Utilizing the Hale Laboratory spectro- 
graph, it was determined that the filter met specifications 
and it was stored awaiting completion of the telescope cluster 
at BBSO. VThen it became possible to check out the 1083 0 A 
filter at BBSO using the spectrograph there, no transmission could 
be observed at all using the silicon vidicon camera. The filter 
was returned to Spectra Optics for examination. It was found 
that the polaroids' had all deteriorated and had become 
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individually opaque. Whether this was a reaction with the oil 
immersion bath or simply deterioration due to extended storage 
at 46 °c was not determined but the latter reason is suspected. 

Polaroid supplied Spectra Optics with an improved substrate 
material which was tested at elevated temperatures for a period 
of months with no detectable change in transmission. With new 
polaroids made from the new substrate material the filter has 
been reassembled. It is now undergoing alignment. Fine tuning 
will be completed by October 1, 1974 at which time the 10830 A 
filter will be redelivered to Caltech for acceptance testing 
at BBSO. However, due to the unfortunate and unsuspected 
substrate materials problem, no 1083 OA observations were obtained 
during the Skylab/ATM mission. Results obtained through the 
1083 OA filter will be reported on as they become available, 

7 . Controls 

Controls have been essentially totally reworked from the 
original telescope due to changes and additions to drive functions, 
camera controls, filter controls, offset capability, flip mirror 
positioning and the like. A controller has been built for the 
65 cm FVU system that will allow programming, in the proper 
sequence, flip mirrors and cameras once a preset frame rate 
and bench holding time have been selected. The controls and 
intercabling work have not been completed even though the 
telescope systems are operational. Modifications, improvements, 
and additions will continue to change control needs. Computer 
control, while temporarily shelved, will eventually be 
implemented, possibly to simulate control functions for astronauts 
on early shuttle/pallet missions. The controls, like the 
instrumentation will always be dynamically changing as improve- 
ment in filters, cameras, electronic circuitry, etc. occur. 
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8. Installation 


To allow installation by a truck mounted crane (rather than a 
barge which would be far more costly) the causeway at BBSO was 
widened and improved prior to teardown of the old telescope. 
^®9-^down of the old telescope was a cooperative project among 
scientists, staff and summer students. It required approximately 
two days. Shipment from Caltech to BBSO and crane installation 
'• of the major new telescope and mount mechanical elements on the 
existing declination bearings required just one very long day. 
Mechanical assembly and alignment was completed within two -months 
although some redesign was necessary along the way. The 
installation of cabling and controls as well as diligently 
checking out the equipment during this process has never 
really ended. The 8.6" Singer Link Refractor and the 10" 
vacuum refractor were operational during the final Skylab manned 
mission. Although data was obtained from the 65 cm FVU during 
this time frame, seeing problems coupled with frustrating 
delays in obtaining proper alignment and focu s of all optical 
elements prevented obtaining really high resolution imagery 
until late spring and early summer of 1974, nearly a year behind 
schedule. Imagery obtained now makes it all seem worthwhile but 
there were many times in between where progress seemed agonizingly 
slow to all of the staff who devoted long hard hours to this 
impressive task. The accompanying photographs will testify to 
what was accomplished and the qualxty of data presently 
obtainable. 
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PREFACE 


This is the final report for contract NAS 8-28021, "Filter Studies in Hel and 
Call," submitted in accordance with the terms of the contract. Herein, we 
describe ground-based observational programs conducted during Skylab using 
filters constructed for the line of Hel and the infrared Call line at 

S3k2k. 
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D3 events associated closely with rapid spot 
changes in a growing region, 5 August 1973. 

On the next day, 6 August 1973, several small 
spots suddenly appear, with an Imp. I flare. 

Rapid spot growth occurs again later in the day 
with a second D3 event. 
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Figure 
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XI 

XII 

XIIIA. 

XIIIB 

XIIIC 

XIIID 

XIV 
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Small spots suddenly appear and grow at either end of 
a event. 

Under more questionable seeing conditions, a small spot 
appears in the penumbra of a larger spot, near the 
terminal of a event. 

A small filamentary appendage accompanies a dying spot, 
while nearby a spot grows, with a D 3 event. 

The Dg chromosphere appears as an elevated shell 
(White, 1963 ) very much brighter above an active 
region. 

Call active regions seen near the west limb. 

A flare-associated wave in Call (85422) and Hi. 

Bottom row is Hi + 1.22. 

Beginning of Imp. IB flare on 5 September I 973 with 
eruption of a filament. 

Erupting filament associated with flare on 5 September 

1973. 

Surge developing with the flare and erupting filament 
on 5 September 1973* 

Slowing and reversal of surge direction - after the 
end of the flare on 5 September 1973. 

Flare in Call showing suggestion of fine structure. 

Dark borders are seen around plage of all active 
regions in Call at 8542 a. 

Development of an active region at 85422. 



SECTION 1 


the Hel FILTER OBSERVATIONS 


1.1 


Summary of the Hel Rm p ^ram 


The D3 Hel eyetem began flrat operation In April, 1973, and hae fnnotloned 
eueeesafully, with aooe changea and improvements added between manned Skylab 
« seions, to the present. April, 197U. We obtained over 12,000 feet of 35 „m 

time-lapse films in P3 at mostly 15 sec. intervals, with emphasis in obeerr 
on during the three manned periods of Skylab. 

The first D3 filter completed was a M (half-width) blrefrlngent filter 

in *b r ™ve plates, prepared and mounted 

in “ t==tal9u. aesorlbed by Title (l97h). This filter began operation in 

sZ^ b r first and second manned 

missions. After the second mission, the D3 filter was completely re- 

n a new design utilising alternating inefficient polaroids, developed 
and described by Title flQYU^ Tnm-Sc 4 

increase of betL Z ‘ ‘ransmlsslon 

trast. At this time the filter was also modified to provide a continuous 

ref at ,37. t .8X. Per the final Skylab missio:, we obtaild . 

taneous pictures at 15 sec. Intervals at D, and D, t .68. 

some technical difficulties were encountered during missions, but in general 

he D program has been exceptionally successful It has • a 

observafifin^i ,4 .. 4. successi-ul. It has yielded a wealth of 

Observational data for analysis and correlation with Skylab observations. 

For the final Skylab mission we also set up a full limb eclipsing system in 
Dq record prominence activity in D wi+h wv 4. 

active prominences. ^ especially ascending and 


144 



1.2 


Hours of D. 


Observations 


On the following page is a list of days with hours of observation with the 
large-scale system during the three manned Sl^lab missions. 


1.3 Preliminary Results 
1.3.1 Visible Activity in 

A detailed comparison of all events with Hx has not yet been made. Such 
comparison was deliberately omitted in hopes of more effectively assessing 
the nature of visible activity in without bias from the more visible Hv 
flares . 


For this reason we use the term events on activity except in the obvious 
cases of major flares. 

Emission 

Less than 10^ of all events observed with a .4.X filter against the solar 
disk are visible at any time in emission. Intense emission well above con- 
tinuum intensity is even less frequent and limited to only major flares 
classed invariably as bright in Hx . 

" Bright emission ” in can be described as four distinct features of flares: 

(1) Bright localized centers. These abruptly with flare start and occur 
adjacent to sunspots (lMSC/695791) (Harvey and Milkey, 1972) with life- 
times on the order of 20 minutes. 

(2) Moving ribbons. Two bright parallel ribbons start close together and 
move apart, completely analogous to many flares observed in Hy (Dodson, 
i960) (Morton and Malville, 1963). In Figure lA, the D^ flare ribbons 
expand or move apart at a velocity of 8 km/sec., consistent with Hh 
flares observed. 
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11 APRIL 1973 


K— 1 XIO^KM. -)l 


1852:00 


1845:10 


1847:55 


1853:00 


1856:00 


1904:30 


^ Figure lA.. Two bright ribbons of emission in the Imp. IB flare of 11 April 1973, 

8 km/sec. D 3 absorption obscures underlying sunspots. 


move apart at 
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2013:00 2022:30 2032:00 



2040:10 2105:00 2132:00 


Figure IB. Classed as linp. IIB, the second flare of 11 April 1973, has extensive ejecta (2022:30) 
and is visible entirely in absorption. The two later sub-flares at 2105:00 and 2132:00 
are noticeably darker in D 3 absorption than the second major flare. 


11 APRIL 1973 W — 1 X 10 > KM. — M 



2030 30 2032 00 2033 30 

Figure IC. Enlargement of the Imp. IIB flare, 11 April 1973> especially at 2032:00, resolves fine 
parallel filaments connecting dark strings of points. One string, the lower, moves as 
a sharply defined front, away from the other string. 


VD 


(3) Multiple arch condensations. A localized parallel system of bright 

arches appears in the post maximum phase of a major flare and coincides 
almost certainly to the tops of a complex "coronal loop system .” 

(Figure II) 

(^) Short- lived emission centers. These low intensity emission centers are 
observed adjacent to spots at the initial phase of major flares and also 
with the beginning of emission loops (Figure II ) limb. 


The visibility of "faint emission " in flares may increase as a function of 
proximity to the limb. Several sub-flares near the limb were seen in emission 
indicating that a sizeable number of events are invisible against the 
brighter center of the solar disk. This may be the case in Figure III. The 
red— shifted bright feature in the Hi? spectrum (D)j typical of many Ry flares j 
in not visible at + .sS (A). Note the brightest part of the Hx flare (C) 
does not coincide spatially with the darkest absorption (b). 


^3 


Absorption 


Almost all Har flares of any consequence probably have a period of limited 
visibility in in absorption when observing with a filter of .5^ or less 
and under good seeing conditions. We judge this by the large numbers of D 
events observed with the improved filter of .4^ half-width. 


3 


Besides long-lived absorption features in such as bright plage or active 
filaments, there are five distinguishable types of distinct, transient 
activity. 


Arch condensations. The majority of small events are observed as arch 
shape condensations, with one or both ends terminating at sunspots. 
With many events observed simultaneously at and 
absorption is best observed at line center, and is rarely visible at 
+ .8& (Figure IIIA and B). This differs from flares which 
exhibit a marked preference for the red wing. 


150 



151 





1932:45 


1934:30 


1958:00 


1958:45 


1936:45 


2000:00 


2001:00 


1948:45 


2002:10 2003:15 2007:00 2039:00 

Figure II. A pre-flare filament darkening (193^230) abrupt disappearance with an Imp. IIB 
flare. Bright arches appearing abruptly in the post-maximum phase (2003:15) are 
probably the tops of a "complex coronal loop system. " 







M- 1 X1Q5 KM.-M 



2249:00 2252:00 2256:00 

Figure III. At Dj + .8X (A) sub-flares are rarely visible. They are best observed at "line center" 
(B). The Hy flare is red-shifted (d) at a place where Do is darkest at line 

center. The brightest Hr (C) and darkest do "not" coincide. 






Moving ribbons < With major riares ve observe large areas of slowly 
changing absorption features frequently as two parallel ribbons moving 
apart. Such features coincide extremely well with the emission, 
Figures I and II. 

( 3 ) Fine threads > In high resolution, Figure IC, at 2032 U.T., many absorp- 
tion features in a major flare resolve into fine parallel filaments 
connecting dark strings of points# One string of dark points (foot^ 
points) advances as a sharply defined front moving away from the other. 
This agrees well with some Hx flares observed in high resolution, with 
the difference that displays the very fine filament structures 
resembling coronal condensations . 

FoQt -points# The darkest observed absorption features in single 
pictures, easy to confuse with sunspots, are frequently the apparent 
foot-points of condensations. These foot -points frequently occur 
adjacent to sunspots. 

( 5 ) Erupting filaments. The "pre-flare filament darkening" generally ob- 
served in the blue wing of Hy (Smith emd Ramsey, 1964) (Martin and 
Ramsey, 1973) bas been observed with three major flares in Exactly 
as in Hy we observe it as a pronounced deirkening, expanding arch, and 
abrupt disappearance, Figure II. In one of three flares, we observed 
the breakup of the filament and ejection in a transverse direction at 
high velocity. 


1*3.2 Pre-flare Filament Darkening and the Brightening of Ascending 
Prominences in D* 


The physical mechanism for observing a brightening of ascending prominences 
in Ky presented by Ifyder, 1970 j is one of resonance scattering due to the 
increased photospheric radiation that a Doppler shifted feature sees when it 
is experiencing increasing radial velocity (Doppler brightening). In D^, we 
have no dense absorbing chromosphere and thus no means of producing increased 
resonance scattering by radial motion. Consequently, an ascending prominence 
in D^ should not brighten as much as Hy if resonance scattering is the primary 
mechanism. 
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If we interpret the pre-flare filament darkening against the disk as bright- 
ening of an ascending prominence (Smith and Ramsey, 19^4), then does in 
fact, at least qualitatively, behave very much like Hy (Figure H). 

To further test the hypothesis of Doppler brightening we set up a full limb 
eclipsing telescope to observe D^ prominences simultaneously with an Hy 
prominence patrol. Several ascending prominences were observed, which were 
judged to be moderately slow, and with these very little brightening was 
observed in either D^ or Hx. With all active, rapidly moving prominences, 
both D^ and Hx became significantly brighter during periods of increased 
prominence motion. With the one classic ’’rapidly” ascending prominence, 
Figure IV, observed to date, the prominence brightens in a way which 
appears identical to Hv. We plan to micro-densitometer this event to more 
carefully determine their relative brightness increases. 

The obvious weakness in these observations, beyond exposure difference, is 
that the Hx system has a iR filter and cannot contain such strong Doppler 
shifts as the SR D^ filter system. 

We plan to rectify this weakness in a new system using two eclipsing tele- 
scopes and essentially identical D^ and Hy filters. 

1 . 3.3 Sunspot Changes and D^ Activity 

Except for the possible temporary obscuration of sunspots by D^ transient 
events, the scarcity of the D^ chromosphere offers a convenient means of 
relating D^ events with photospheric sunspots and their changes. As anti- 
cipated, we find a strong correlation between the frequency of D^ events and 
changing sunspots. This is especially true with the appearance and growth of 
new spots (growing or emerging active regions). 

During the second Sltylab mission we observed the rapid growth of such a 
region, associated with numerous D^ events. The most interesting aspect of 
the activity is the suggestion of individual flares associating with rapid 
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190 1 



1908 1917 




Figure IV. A rapidly ascending prominence observed in Do (upper), brightens dramatically in a 
manner qualitatively consistent with that observed in Hv (lower). 



changes in small umbrae. We have observed almost continuous activity 
adjacent to and over a rapidly growing sunspot for a period of several hours. 
The suggestion is, however, that sunspot growth may occur in spurts associated 
with the individual events. On August 5, 1973 > Figure V, we point out 
several small spots changing rapidly, appearing to coincide with individual 
events. With the first event we observe one spot decrease in size while 
another grows, with the event spanning between. This is a rapidly growing 
region, and on the next day, August 6, 1973, Figure VI, we observed several 
small spots appear and grow rapidly with a large event. Rapid spot growth 

did not take place again until later in the day when we observed a second D 
event . 


In Figure VII, October 15, 1973 j we observe again rapid appearance of small 
umbrae associated with a event. 

In Figure VIII, September 1, 1973, under more questionable seeing conditions, 
we observe the appearance of a small spot in the penumbra of a large spot, 

CO incident with a event . 


Occasionally, small 
Figure IX. 


spots can both grow and die coincident with D 


3 


activity. 


We have shown several selected events with excellent conditions to show rapid 
changes in tiny umbrae, exactly coincident with activity. Much more 
frequently , even with major flares and excellent conditions, no spot changes 
could be detected on a time scale relating directly to the flares. 

Our few, but interesting rapid changes in tiny umbrae coincident with 
activity , point to the need for observing the "strong" sunspots magnetic 
fields with the high spatial and time resolution necessary to monitor such 
changes. Present magnetograph techniques do not sufficiently emphasize such 
strong fields of small spatial extent. With the recent introduction of the 
diode array to raagnetograph techniques, it is evident that such observations 
will soon be more practical to obtain. 
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Figure V. 


5 AUG. 1973 


1 X 105 km. 




Dj events associated closely with rapid spot changes in a growing region, 5 August 1973. 


6 AUG. 1973 K — 1 X lo® km. — w 



1718 1724 1732 1742 1752 1810 1951 


Figure VI. On the next day, 6 August 1973j several small spots suddenly appear, with an Imp. 

flare. Rapid spot growth occurs again later in the day with a second event. 


15 OCT. 1973 I 30,000 km I 




1745:00 1811:00 1815:00 



1824:00 


1837:00 


1858:00 


VO 


Figure VII. Small spots suddenly appear and grow at either end of a event 




1655:00 1722:00 1737:00 



1742:00 1747:00 1752:00 


Figure VIII. Under more questionable seeing conditions, a small spot appears in the penumbra of 
a larger spot, near the terminal of a D3 event. 


29 SEPT. 1973 K — 1 X 10® KM. 



H Figure IX . A small filamentary appendage accompanies a dying spot, while nearby a spot grows, 

I-* with a event. 


1.3.4 The Dj Chromosphere at the Limb 

One of the most interesting aspects of when observed with a narrow band 
filter at the solar limb, under good seeing conditions is the chromosphere. 
It is concentrated in a raised shell above the photospheric limb. We measured 
the center of this shell to be 1450 km above the photosphere, confirming 
measurements of this approximate height made by White, I 963 . In Figure X, 
we picture the raised shell and how it brightens markedly above an active 
region directly on the limb. Visually, the raised shell is evident 
around the entire limb, but brightens as expected over active regions. The 
increase in brightness over active regions is much greater than observed for 
Hr. If this increased brightness of over active regions corresponds to an 
increased number of Hel atoms and ions this is a possible source for the in- 
creased He abundance observed in space with particles from major solar flares 
(Hirshberg et al., 197 I). 
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Figure X. The chromosphere appears as an elevated shell (White, I963) very much brighter above 
an active region. 


SECTION 2 


THE Call observations AT 8542^ 


2.1 Objectives 

A I/4S filter for the infrared lines of Call was built and used for ground- 
based observations during the Skylab missions. This was an experimental pro- 
gram in that this is the first filter that has ever been constructed for the 
infrared Call lines. Our original goal was to take time-lapse observations 
in the Call line at 85428. This objective was achieved during Sky lab III 
although not during Sky lab II or IV. Both the successes and difficulties of 
this program are described below. 

2.2 Description of the Filter 

The functioning of this Call filter built by Spectra Optics, surpasses our 
original objectives in two respects: (l) The filter is tunable over several 

hundred angstroms. Thus, it is possible to observe all three of the infrared 
Call lines as well as other interesting lines nearby such as the Fel line at 
84688. (2) Two of the infrared Can lines (85428 and 84988) can be observed 

simultaneously or in rapid succession by simply rotating the 84988 and 8542S 
prefilters alternately into the beam. During Sky lab in, photographs were 
only obtained at 85428. However, visual assessment of 84988 and 85428 were 
made. No differences could be seen in the two lines, however, it should be 
noted that visual observations are limited by the quality of image obtained 
on the TV monitor which shows an interference pattern on the face of the tube. 
The narrowest passband that the filter can deliver is I/48. However, the fil- 
ter is versatile in that it can be used as a .378 filter or a .58 filter sim- 
ply by removing one or both polarizers on each end of the filter. The filter 
can also be used as a Doppler filter by tuning the 1/48 element out of phase 
with the rest of the elements. This yields two passbands slightly narrower 
than 1/4 angstrom, separated by a I/48. If the two passbands are placed 
symmetrically around a line, one can then observe both ascending and descend- 
ing mass by use of a polarizing prism behind the filter. 
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During the construction of the filter we requested two design alterations of 
the manufacturer, Spectra Optics. Our request was granted with no increase 
in cost but it did result in a delay in completing the filter. The changes 
were to: (l) Change the thicknesses of the calcite to provide two passbands 

centered as close as possible to the lines at 8498^ &nd 85428 rather than just 
the single passband at 85428. (2) To incorporate polarizing beamsplitters 

supplied by us into the filter in place of the sheet polaroid usually used 
in birefringent filters. This was an important change in order to increase 
the transmission of the filter. The polarizing coatings were purchased from 
Perkin-Elmer Corporation using a design specified by Dr. Alan Title, 

2.3 The Observational Program 

Testing of the completed filter was begun in June, 1973. Two problems were 
encountered. The filter would not tune properly and the polarizing prisms 
were found to scatter an excessive amount of light. Spectra Optics took the 
filter back to solve the tuning problem and Perkin-Elmer agreed to replace 
the defective polarizing coatings. 

Routine observations were begun on August 8, I973 using the I/28 passband. 
Useful observations were taken until after the end of the second manned 
mission (Skylab III) using the original set of polarizing prisms. 

The replacement prisms were delivered in early December 1973. However, it 
was not discovered until January 1974, when the new prisms were to be installed 
in the filter, that they did not meet the size specification and were too 
large to fit into the filter . This posed an additional problem. The polar- 
izing coatings are soluble in water. Therefore they had to be cut and re- 
polished in unconventional glass working vehicles such as oil or kerosene. 

The reworking of the prisms was accomplished after considerable experimen- 
tation. The final polish was achieved using ultra-low viscosity silicone oil 
as an abrasive vehicle . The filter was finally reassembled after the end of 
Skylab XV. 
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Photographs taken through the filter with the second set of polarizing prisms 
are significantly higher in contrast and apparent resolution. Figure XI 
shows two of the most recent photographs taken through the filter using the 
. 37 S passband centered on 8542^ line. These photographs were printed on the 
lowest contrast paper available, Agfa Rapid Print No. 1. In comparison, the 
earlier photographs , taken with the original set of polarizing prisms in the 
filter, have to be printed on paper 3 grades higher in contrast to achieve a 
similar appearance. 

The photographs in Figure XI were taken with the filter in a telescope having 
an excellent objective lens 7 inches in diameter. In contrast to Hel and 
Ha photographs taken through the same telescope, the Call photographs show 
less sharply defined fine structure. To test the visual sharpness of sun- 
spots through the Call filter, the 8542 X blocking filter was temporarily re- 
placed by an Ha blocking filter. The Call filter delievered images equal in 
sharpness to the observations made on the same day. Thus, the lack of 
sharply defined fine structure in the Call photographs relative to Ha is 
probably not due either to the filter or the telescope. The differences can 
be attributed to the differences between the IV N emulsion used for the Call 
photographs and the ultra fine grain, high contrast, S0392 emulsion used for 
Ha and to the relative intensity gradient between adjacent structures in the 
two lines . 

IV N, a special order emulsion made by Kbdak,was selected because it is the 
finest grain emulsion available for infrared wavelengths. I N and commercial 
high speed infrared films were also tested. Although these emulsions are 
much faster , the grain is so coarse as to inhibit high resolution photography 
of whole active regions on 35 film. During testing of the films, we were 
advised that the speed of the IV N emulsion could be increased by hypersen- 
sitlzing the film in a bath of water at approximately 3°C for a few minutes. 
We followed this recommendation and found that the speed increase was greater 
than a factor of 4 without appreciable increase in grain size. This reduces 
the normal exposure through filter with the l/4X passband from 1 sec. to less 
than 1/4 sec. 
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16 APRIL 1974 



17 APRIL 



Figure XI. Call active regions seen near the west limb. 
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2.4 Preliminaiy Results of the Call Observations 

On September 5j 1973 j an exceptionally good solar flaxe was recorded on board 
Skylab III - simultaneously with our ground-based observations on the Call 
filter, our Hy Doppler filter and multi-slit spectrograph. This is probably 
the best flare event recorded during all three Skylab missions. Although 
the flare was not exceptionally large (Class IB) it was associated with a 
filament eruption, a flare-associated shock wave and a large surge. 

The Call observations show no sign of the filament eruption and only a small 
part of the surge. The wave, however, is as bright or brighter than the 
bright wave seen in Haf . The arrows in Figure XU point to the curved wave 
front as it moves out from the flare. The apparent velocity of the wave in 
the Call and Hy slit jaw photographs is ~ 400 km/sec. The wave is not visible 
in the Har Doppler observations at ±1.2^ in the wings of the iix in Figures XII 
and XIII, nor is it readily discernible in the multiple array of Ifa spectra 
(Figures XII and XIIl). To our knowledge no other flare -as so dated waves 
were recorded during any of the Skylab missions. Several other small flares 
were recorded in the Call filter but without any sign of associated events. 

Another flare recorded with the Call filter was under-exposed, and is the only 
flare in which we can see that the flare has apparent fine structure. This 
flare, shown in Figure XIV is composed of several small elements. Each 
element has a relatively short lifetime. 

An interesting feature of active regions in the 85428 line is that the bright 
plage is surrounded by a region about 30,000 km wide which is darker than the 
Call quiet chromosphere. Several examples eire shown in Figure XV. The dark 
band is visible even for relatively young active regions and persists through- 
out the life of the active region. However, for the oldest, most spread out 
regions, the dark band seems to be broader and less pronounced. The band 
apparently corresponds to the outer parts of active regions in where the 
fibrils appear to lie horizontal to the solar surface. The dark band appears 
to increase in width with the develoiment of an active region (Figure XV). 
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1827:45 


1829:30 


1830:45 


183130 


Figure XII. A flare-associated wave in Call (8542X) and . Bottom row is Ifa + I.2X. 



5 SEPT.1973 



I828 M5 U.T. 1829:30 I830'0b 


Figure XIIIA. Beginning of Imp. IB flare on 5 September 1973 with the eruption 
of a filament. 
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I830'45 I83|:45 I833:|5 

Figure XIIIB. Erupting filament associated with flare on 5 September 1973. 


171 







5 SEPT. ’973 



1835:30 1837:45 I843:|5 


Figure XIIIC. Surge developing with the flare and erupting flJameiit on 5 
September 1973 • 
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Figure XIIID. 


Slowing and reversal of surge direction 
flare on 5 September 1973 • 


- after the end of the 
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8 AUG 1973 |_,o» km - I 
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2351:25 2359:25 2410:55 


Figure XIV. Flare in Call showing siiggestion of fine structure. 
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Figure XV. Dark borders are seen around plage of all active regions in Call 
at 8542S. 
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The Call at 854-2S observations reveal the existence of umbral flashes. On 
normal exposures, however, they are at the limit of visibility. Much longer 
exposures will be req.uired to make these flashes readily visible. Doppler 
observations of the flashes would be worthwhile for comparison with Hy umbral 
oscillations. 

The development of active regions in Call is illustrated in Figure XVI. A 
sunspot in an extensive region appears to change very little in the 7 day 
period. In contrast a new region, initially seen on 2? August, undergoes 
great change. This region, probably new on 26 August, is already losing its 
sunspots and decaying by 30 August. Another new region appears on 30 August 
and exhibits a rapid increase in size with the development of sunspots. By 
31 August the region begins to take on the appearance of the background net- 
work, revealing the presence of super granule cells which change very slowly. 
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Figure XVI. Development of an active region at 8542X. 


177 






REFERENCES 


Hlrshberg, J. , Asbridge, J.R. and Robbins, D.E.: I97I, Solar Physics 18. 313 
Ityder, C.L, ; 1970, Solar Riysics 1^, 147 

Martin, S.F. and Ramsey, H,E,: 1972, Solar Activity Observations and 

Predictions . ed. Murray and McIntosh, p 371 

Smith, S.F., and Ramsey, H.E.: 1964, Z, Astrophvs. 60. 1 

Title, A.: 1974, Submitted to Solar Riysics, "Partial Bolaroidn in Rirp- 

fringent Filters” 

Title, A.: 1973, Solar Physics 33 , 521 

White, O.R.: 1963, Ap, J. 138, 1316 


178 



Appendix E 
UCSD Rep6tt 

Direct Observation of Temperature Amplitude of 
Solar 300-Second Oscillations 
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Direct ObBervatlon of Ten^rature An|>Iitude of 
Solar 300-Sec Oec Illations 


H. S. Hudson and C. A. Lindsey 

Department of Physics 
University of California, San Diego 
La Jolla, California 9^037 


ABSTRACT 

nie 300-second oscillations form the dominant source of variability 
of the solar Infrared continuum. We have observed them at 20 g. with an 
an^jlltude <= 3.0 °K over an area with an effective diameter of 33". 

This new mode of observation of the 300-sec osc Illations should make 
possible a fundamental inqprovement In our knowledge of their nature and 
origin. 
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Intensity variations in the Rayleigh- Jeans portion of the solar 
continuum spectrum, longvard of 1 ji, permit direct observations of the 
ten5»erature structure of the solar atmosphere. Table 1 indicates the 
source heights in the HSRA. model atmosphere (Glngerlch, Hoyes, Kalkofen 
and Cuny 1971) for several of the infrared windows available for ground- 
based observations. The existence of 300-sec oscillations in the infra- 
red continuum provides a powerful new tool for the exploration of this 
phenomenon: variation with altitude, spatial structure, relative phase 

of temperatures and velocity variations, etc. 

We report here the Initial observations at 20 ii, made with a two- 
beam photometer at the UCSD-Unlversity of Minnesota 60” (152 cm) Infrared 
telescope of the Mt. Lemmon Observatory (altitude 28ol+ m). The high level 
of sky noise and atmospheric transparency fluctuations at Infrared wave- 
lengths req.ulre the use of high, dry site, and a two-beam photometric 
system to effect sky cancellation. The observations reported here, 
obtained 1973 June 4 near local noon, utilized circular beams of 23" 
diameter, separated by ~ 2', and roughly centered at NIO EI 9 . An aperture 
filter of black polyethylene helped reduce heating of the telescope. An 
interference filter (OCLI L19420-8) plus a KES-5 Dewar entrance window 
(long -wavelength cutoff ~ 4o |x), plus the Earth's atmosphere, determined 
the spectral passband as approximately 5 ^ centered on 22 p.. Subsequent 
observations under inferior conditions required the use of a KCl window 
(cutoff ~ 22 ^) to suppress sky noise. The rms aitqplltude of the variations 
of the difference signal between the two beams on the solar disk amounted 
to about 0 . 06 ^ of the signal obtained with one beam on the disk and one 
off, for a 1 Hz bandpass. This limb reading calibrates the signal strength 
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at AT =6.3 °K for an aBsumed brightnefis temperature of U9OO ®K. The 

rms 

correepondlng peak-to-peak difference signal vlth both beams off the Sun 
amounted to 0 . 006 $, a noise level presumably corresponding to sky noise 
or to detector noise. Most of this noise occurs at higher frequencies 
than the 300-sec oscillations. Spatial scans of the tengperature structure 
of the solar atmosphere showed a high degree of reproducibility for short 
time Intervals, so we have little doubt of the solar origin of the bulk ^ 
of the long-period variations. We assume that the variations In the two 
beams add independently and that the beam size exceeds the size of an 
oscillating element, so that the effective beam diameter becomes x 23" = 33". 

The peak frequency of the oscillatory feature in Figure 1 lies at 
3.*»- X lO”'^ Hz, which corresponds closely to the 300-sec period expected 
from observations (e.g. Noyes 1967) of the velocity field. The area under 
the peak In Figure 1 yields an rms ten^ierature fluctuation of 3.0 ''K, 
after integration between 2.5 euid 4.1 x 10"^ Hz. The rms an5>lltude does 
not have much quantitative significance because of the Intrinsic varia- 
bility of the phenomenon. However, simultaneous temperature and Doppler 
measurements should enable us to determine the relative anplltude and 
phase of individual oscillating elements. We have observed the 300-sec 
oscillations at other wavelengths, ranging from 8.6 p, to 350 pi. 
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TABLE 1 


Wavelength 

(microns) 


HSRA Atmosphere 


Optical Depth Height 

(at 5000 A) (kilometers) 


0.5 

1.00 

10.0 

0.12 

20 

3.5 X 

35 

1.3 X 

350 

1.6 X 

1200 

7.2 X 



0 


130 

10-2 

210 

K)-2 

270 

lo-** 

530 

10"^ 

930 
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FIGURE CAPTION 


Figure 1 Pover spectrum of an 83-minute data Baii^>le at 20 p, obtained 
1973 June U, computed from an autocorrelation function U1.5 
minutes in length. The observations show ten5>erature fluctuations 
in an effective aperture of 33" diameter, with an rms amplitude 
of 1.5 ‘’k between 2.5 and 4.1 x lO”^ Hz. < 
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ABSTRACT 


A theoretical stu^ Is made of the visible and UV line radiation of 
He I atoms and He II ions from a plane -par el^lel model flare layer. Codes have 
been developed for the solution of the statistically steady state equations 
for a 30 level He I - II III model, and the line and continuum transport 
equations* These codes are described and documented in the report along with 
sample solutions* Optical depths and some line intensities are presented for 
a 1000 km thick layer • Solutions of the steady state equations are presented 
for electron temperatures 10 ^ > 5 x 10 ^ °K and electron densities 

,ao -5 

10 - 10 cm . 
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HELXU^j MISSION FROM ACTIVE SOLAR REGIONS 
I. INTRODUCTION 

The purpose of this program has been to develop codes for the 
simultaneous calculation of He I and II resonance line and He I line 
Intensities from model flare regions. These lines vere chosen because 
of the spectral ranges of the Sky lab high resolution spectrographs and 
because of the planned program to obtain flltergrams on a patrol 
basis at the Lockheed I^e Canyon Observatory during the ASM mission. 

The NHL spectrohellograms incorporate simultaneous measurements of the 
He I and II resonance lines. 

A plane-parallel layer Irradiated on erne side by the photospherlc 
radiation field vas chosen as the geometric model* A statistically 
steady state and uniform electron temperature and density vith position 
vere assumed. The energy level model consists of all terms through prin- 
cipal quantum number 4. Our study has been confined to conditions we 
believe characteristic of flare regions, namely electron temperatures be- 
tween 10*^ and 5 x 10^ °K and electron densities between 10^® and lO^^cm”^. 
An extensive compilation of electron Impact excitation rates has been 
made as part of this study. The results were published In Solar Physics 
(Benson and Kulander, 1972). 

The statistically steetdy state level populations of model He I 
atoms have been calculated by a number of Investigators for temperatures 
and densities characteristic of the outer solar atmosphere. Almost none 
of these authors has considered a evifflciently detailed energy level 
structure In the model atom to ekcourately obtain the line emission. 
Jefferies (19^^) e.g. treats the 2a and 2p levels as a single level. 

De Jager and de Groot (19??) consider the 2s and 2p terms sepenrately 
but the term structure of higher levels Is Ignored. This higher term 
structure Is also Ignored by Athay end Johnson (i960). They also neglect 
the effect of the He II Ion processes by using other values for the 
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He l/He II equilibrium. Athay and Jdinson arrive at the conclusion that 
In the temperature range 40,000 - 50,000°K, the D3 line will appear In 
emission for n^ > 10^ almost Independently of T^. 

Zlrln (1956) assumes In his calc ul ations that transitions between 
terms of a given level are of negligible importance in determining the 
occupation numbers. This is known to be a poor assumption, Shklovsky 
and Kbnonovltch (1958) have calculated the D3 line Intensity but have made 
a number of unrealistic physical assumptions. More recently Hearn 
(1969) has calculated the occupation numbers of a 41 level He I atom and 
one level He II Ion but he only j«*osents results for the resonance line 
intensities. We shall demonstrate that more levels axe required In He II 
to obtain correct line Intensities. 

Jefferies (1957) has calculated the D 3 line Intensity from a layer 
assumed to be optically thick in the D3 line. The transport equation was 
solved assuming Incoherent scattering with no photospheric radiation In 
the line. Jefferies' results are very qualitative since It Is known that 
the D3 line Is probably not optically thick. 

To obtain accurate line Intensities, simultaneous solution of the 
line and continuum radiative transfer equations and the steady state 
populations Is required. To accomplish this, we have developed three 
separate codes. The first code (Code l) solves the statistical equili- 
brium equations for a 30 level He I-II-IXI system given the appropriate 
rates. The basic equations and sample solutions are given In Section II. 
The code Is described In Section V and a listing Is given In Appendix A. 

The second code (Code 2) represents a numerical solution of the 
line transport equations for a finite layer. The solution is of the 
Integral form of the transport equation by expansion of the source func- 
tion In terms of a finite sum. The mathematical method used is summarized 
by Avrett and Loeser (1969). Complete frequency redistribution and a 
Gaussian absorption profile are assumed. The basic equations a nd sample 
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solution are given in Section III. The code is described in Section V 
and a listing is given In Appendix B. 

The third code (Code 3) solves the continuum transport equation 
by expansion of the source function in a very similar manner to the line 
transport equation. The equations are given in Section IV. The code 
la described in Section V and a listing is given in Appendix C. 
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II. THE STEADY STATE EQUATIONS 

A. Energy Level Model 

The 50 assumed energy levels for the system of ions He I - He III are 
given in Table II-l. There are 19 levels for He I, 10 for He II and 1 for 
He III. These levels are shown in Figures II -1 and II-2. The levels are 
numbered 1 - 50 in order of increasing energy. Levels 1, 20 and 50 are the 
ground states of He I, II and III^ respectively. This model was chosen to 
provide accurate solution for the first two resonance lines in He I and II 
and the D5 and IO 850 lines of He I. Many other lines are Included but were 
not the primary lines under consideration. The allowed transitions included 
in the model are listed in Table II. 2 together with f numbers and wavelengths. 
The inclusion of the 4S,P,D and F levels separately is necessary because at 
the lower electron densities considered radiative de-excitation rates can 
become larger than collision rates between the n ■ levels. At higher electron 
densities the collision rates between the n ■ ^ terms dominate all other rates 
in or out of these terms and the relative populations are Boltzmann. 


Population Equations 

The rate equation describing the population of the bound or continuum 
state i is 


(R - R ) = 0 


II. 1 


where ^yji ^ij ^re the total transition rates/cm^ by process y from state 
j to i and from i to j respectively. The sum S represents a sum over discrete 
states and an Integration over continuum states. We shall assume the particle 
translational dlstrlbution"functions to be Maxwellian and the external 
continuum radiation field to be Planckian> in which case it is possible to 
integrate over the continuum and replace it by one additional term in the discrete 
sum. The atomic transition processes are radiative excitation and ionization, 
excitation and ionization by atoms and electrons and their inverses. 
Because of their higher velocities, electron collisions generally dominate the 
rates and hence only electron Inelastic rates will be considered. 
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TABLE n.l 

He I, II Energy Levels 


He I 

- i = 

1 

... Energy (ev) 

Wave Nos. 

st.wt 

1 

2 

Is 

A 

0 

0 

1 

2 

ls 2 s 


19.821 

159850 

3 

3 

ls 2 s 

^8 

20.618 

166272 

1 

4 

ls2p 


20.966 

169081 

9 

5 

ls2p 

V 

21.220 

171129 

3 

6 

ls 3 s 


22.721 

183231 

3 

7 

ls 3 s 


22.923 

184859 

1 

8 

ls 3 p 

V 

23.009 

185559 

9 

9 

ls 3 d 


23.076 

IB6096 

15 

10 

ls 3 d 

h ) 

23.076 

186099 

5 

11 

ls 3 p 

IpO 

23.089 

186204 

3 

12 

1b4b 

^8 

23.596 

190292 

3 

13 

1s48 


23.676 

190935 

1 

14 

ls 4 p 

3pO 

23.710 

191211 

9 

15 

ls 4 d 


23.738 

191439 

15 

l6 

ls 4 d 


23.739 

191441 

5 

17 

ls 4 f 


23.739 

191447 

27 

18 

ls 4 f 


23.739 

191447 

7 

19 

ls 4 p 

IpO 

23 . 7 ^*^ 

191487 

3 

He II 

- i = 

2 




1 

Is 

2 

0 

0 

2 

2 

2s 

2b 

40.8099 

329179.57 

2 

3 

2p 

2 p 0 

40.8091 

329182.02 

6 

4 

3 s 

2s 

48.3662 

390140.76 

2 

5 

3 p 

2 p 0 

48.3664 

390141.49 

6 

6 

3 d 

2d 

48 . 3665 

390242.64 

10 

7 

4 s 

2s 

51.0113 

411476.98 

2 

8 

4 p 

2 pO 

51.0114 

411477.28 

6 

9 

4d 

2d 

51.0115 

411477.77 

10 

10 

4 f 

2 p 0 

51.0117 

411477.95 

14 
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TABLE II. 2 


He II Lines 


He I 


Upper 

Level 

Lower 

Level 

Notation 

Ui) 

A 

( 10®/ sec ) 

f 

4 

2 

1 

10830 

.1022 

.5391 

5 

1 

2 UV 

584.4 

17.99 

.2762 

5 

3 


20582 

.01976 

.3764 

6 

4 

10 

7065 

.278 

.0693 

7 

5 

45 

7281 

.181 

.0480 

8 

2 

2 

3889 

.09478 

.06446 

8 

6 


4,30+4 

.0108 

.896 

9 

4 

11 

5876 

.706 

.609 

9 

8 


1.86+5 

1.28-4 

.111 

10 

5 

’ 46 

6678 

.638 

.711 

11 

1 

3 UV 

537.1 

5.66 

.0734 

11 

3 

4 

5016 

.1338 

.1514 

11 

7 


7.43+4 

.00253 

.629 

12 

4' 

12 

4713 

.106 

.0118 

12 

8 


21120 

.0652 

.145 

13 

5 

47 

5049 

.0655 

.00834 

13 

11 


21132 

.0459 

.103 

14 

2 

3 

3188 

.0505 

.0231 

14 

6 


12538 

.00608 

.0429 

14 

9 


19543 

.00597 

.0205 

14 

12 


1.09+5 

.0505 

1.21 

15 

4 

14 

4472 

.251 

.025 

15 

8 


17002 

.0668 

.482 

15 

14 


4. 39+5 

4.16-5 

.200 

16 

5 

48 

4922 

.202 

.122 

16 

11 


19089 

.0711 

.647 

17 

9 


18688 

.139 

1.02 

17 

.15 


1.43+7 

6.01-10 

.0033 

18 

10 


18699 

.138 

1.01 

18 

16 


1.67+7 

4.34-10 

.00253 

19 

1 

4 UV 

522.2 

2.46 

.030 

19 

3 

5 

3965 

.0717 

.0507 

19 

7 


15088 

.0137 

.140 
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TABLE II. 2 (Continued) 


He I 


Upper 

Level 

Lower 

Level 

Notation \( k ) 

A 

(l0®/sec) 

f 

19 

10 

18560 

.00277 

.00858 

19 

13 

1 . 81+5 

5.78-4 

.853 

19 

16 

2.17+€ 

5.65-7 

.024 


He II 


3 

1 

3o3.8o 

100. 

.4162 

4 

3 

1640.5 

1.01 

.01359 

5 

1 

256.3 

26.8 

.07910 

5 

2 

1640.4 

3.59 : 

.4349 

6 

3 

1640,4 

10.3 

.6958 

7 

3 

1215.1 

.413 

3.045-3 

7 

5 

4687.0 

.294 

.03225 

8 

1 

243.03 

10.9 

.02899 

8 

2 

1215.1 

1.55 

.1028 

8 

4 

4686.8 

.491 

.4847 

8 

6 

4687.2 

.056 

.01099 

9 

3 

1215.1 

3.30 

.1218 

9 

5 

4686.9 

1.13 

.6183 

10 

6 

4687.1 

2.21 

1.018 
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We may write the total rate from 1 to j more specifically as, 

t 


R . = n . P . . 

yij 1 IJ 




II.2 


where A* ^ and C. . are the radiative and collisonal transition rates/particle 
ij ij 

from the ith to the Jth state. In the statistically steady state, equation 


II. 1 reduces to 


S 





j 


0 ; 




II. 5 


The system of equations represented by II. 3 can easily be solved for the 
populations n^ when all the known. Ihis is the case when the gas 

is totally optically thin (in all lines), since the external radiation field 
is specified and the internal radiation field does not produce a significant 
upward transition rate. We can characterize the system of linear equations 
(II.5) by a matrix whose coefficients a^^j are equal to The diagonal 

elements are the P^^^. 


Bie general solution of equations (II. 5) is given by White (1961) 


n 




II. 4 


mi ^ 

where P is the co -factor of the coefficient of n. in the mth equation 

J X 

(i.e. the matrix element total number of He particles/cm^. 

We denote the discrete level corresponding to the continuum, i.e. an ionization 
or recombination, by c. We may then write the rate equation for level i of 
ion q as. 


where J and J ' refer to other bound levels in the qth ion and in other stages 
of ionization, respectively, P^. is the total excitation or de-excitation rate 
In ion q from i-J; «-nd are the total ionization rates from i in 

q to J ' in q + 1 and to 1 in q from J ' in q - 1, respectively; P^^^ , and 


200 



c*o+l 

total reconiblnation rates from i in q to j’ In q -i and to i 
in q from J’ in q+1, respectively. Specifically we may write the P*s as follows: 

bound -bound 


1> J 
1 < J 


to A J , ij 14 r e IJ 


■14 


^14®i4®r “e^l4 


II.6 


bound -free (from initial state q,i) 


ionization . 
recombination 


P<1 

IcJ 


W A 


A <1 




Iji ■ ij. 

Vij' * IJ. 


n.7 


bound -free (to initial state q, i) 
ionization pJ"J^ ° 

recombination P®'J+^ = ff J+J; 


II. 8 


O 

vhere n^ la the electron denslty/cm ; B^ ia the Planck function at temperature 
T^; W is the dilution factor; is a free parameter, A^y B^. are the Elnateln 

tranaltlon probabllltlea; A is the photoionization rate/lon; H J are 

the colllBlonal tranaltlon^^ratea/electron per ion for bound-bound and bound-free 
processea, respectively; njj is the colllsional transition rate/electron^ per 
ion for free-bound recombination; and a^y are the recoznbination, and 
stimulated recombination coefficients, respectively. 


The basic quantities required for the evaluation of the radiative rates 
are the oscillator strengths and photoionization cross sections; for the coUi- 
slonal rates the excitation and ionization cross sections. The Inverse cross 
sections and Einstein coefficients can be obtained from the usual detailed balance 
relations. The quantities n, 0, 0 and a generally have the form 

,r v<a(v)/(v) dv, 
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where v Is the electron velocity, Q(v) the cross section and /(v) the electron 
translational dlstrlhutlon function. A may be written in terms of the photo- 
lonlzatlon cross section a(v) as 

4 „ J . B(v) dv. II. 10 

p is expressed similarly in terms of the stimulated recombination cross section 
b(v). 


Solutions can be obtained for exxy arbitrary line radiation field by 
suitable choice of the parameter Y which is closely related to the net radia- 
tive bracket (NRB). Y varies between 0 and 1 being 0 for a thin layer (no 
external radiation) and 1 for a very thick layer where the line radiation field 
is Plancklan. For a thin line with external photospherlc radiation Y = l/2 B(Tg). 
Various continuum radiation fields can be considered by suitable choices of 
A. 
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Iteaetlnr 

* In the NSS compilation by WlMe ,t «1 rioSfil 

«>d r*^'“ H»rtree-rock calculations of Stewart 

for the 1 f ’T ^ »*' «» I «"■» oalculatlona by Peach (1967) 

we« 1 , : ^ ^ .11 other level, croa. aectiona 

calculated uelng the quantum-defect method of Burgee, and Seaton (l96o). 

Colll.lonal lonlmatlon rate, were obtained from meaeurement. of 

JtToT^^jL*^? ^PP jnted by Kleffer «.d Bunn (1966) for 1 ^ from meaeure- 

^nte of long (1967) for 2 ^ S and from calculation, of Bolder, et . 1 , (X 96 i) 

for the, He II ground state. For other levels the lon^»-f^ * 

Allen (1961). ionization rates were taken from 


of th ®' ..citation croee .action, are crucial to the eolutlon 

the .teady atate equation.. For thl. reaeon we have made an erten.lve 
tabulation and etudy of rate, from many eourc... A paper entitled "aectron 
impact Encltatlon Bate, for Helium" de.crlblng there rate. wa. publlahed In 

the Becember 19 TC l.aue of Solar Phyelc. (Beneon and Kulander, 1972). Hxcltatlon 


exp (-ax ) 
o 


n » at” 

Where ’ ''' " “”**“*'• ^ «' ^ «>. temperature range 

conatdered wa, bOOO-50,000% and for H. 11, lo"* - lo5>K. Bate, between aU 

evel. Of the model of Section II.A were calculated, lie Inver.e rate, for both 

— '-Ulbrlum 

impact rtLrregCr«ri^rHrr«ter ==ncernmg the electron 
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Generally both forbidden and allowed rates are In better agreement 
for higher temperatures than for lower temperatures. The allowed rates are 
In much better agreement than the forbidden rates. The allowed rates from 
ground state or level 2 generally show differences as large as a factor of 
10 for low temperatures and as large as 5 for high temperatures. Between 
other levels (n s 3) the differences are only as much as factors of 3 at 
low temperatures and 50-10(D^ at high temperatures. 


In He I the forbidden rates generally differ by as much as a factor 
of 100 at low temperatures, a factor of 20-50 at higher temperatures. Differ- 
ences as high as a factor of 10^ are noted In a few cases, differences of 
factors of 10^ are not uncommon. If one does not consider the hipest and 
the lowest rate value for each transition the differences are generally 
reduced to factors of about 10 at lower temperatures and 5 at higher tempera- 
tures. Forbidden transitions with lower level n = 2 show differences 2-3 
times less than those with n s 1. 

For transitions between the higher levels only the very approximate 
cross sections of Green (1966) and Allen (1963) ore available for forbidden 
transitions and for allowed transitions those of Seaton (I962), Saraph (196^) 
and Mlhalas and Stone (I968). Comparing the two approximate results for 
forbidden transitions with the other measured and calculated values we find 
the approximate results generally lower by factors of 2-10. For higher levels 
the Green cross sections are generally higher than those of Allen. With re- 
gard to the Beaton and Saraph approximations. It Is not conclusive which is 
more often closer to the other values. The Seaton approximation gives values 
generally higher than the Saraph approximation for lower temperatures while 
the reverse situation holds for higher temperatures. 

For He II there are relatively few cross sections available. The 
largest differences, being about a factor of 5» uMCh lees than for He I. 
For both He I and He II experimentally determined cross seotiona give lower 
rates than the calculated values. 
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The actual colllslonal excitation rates chosen represent mean values 
of those given in the paper, taken from various sources listed there. To 
detemlne the sensitivity of the solution to these rates the statistical equi- 
librium equations were solved after increasing the Individual values of 
the collision rates per electron f) between each state by a factor of 100. This 
was done to determine the sensitivity of the solution to the collision exci- 
tation rates. We chose the optically thin case with T » gOtOOG*^, n « lo^® 
for this test. We Illustrate in Figure II. 3 the effect on the population of 
the levels 1, 2, 4, 9, 11 and 20 of each perturbation In collision rate. Of 
course in many instances increasing a particular by 100 ha,d little effect 
on the populations of the levels mentioned. In Figure II. 3 we show the 10 
transitions 1-J which have the largest effect on each of the above levels. 

The ratio of the population after increasing ^00 to that before the 

perturbation is shown. We note that changes of 10 to ^0 in these Important 
populations result from an uncertainty of a factor of 100 in the rates. 

All the triplet levels of He 1 had about the same response to the change 
in fl. When the rate from the He I ground state (G.&) to triplet levels increased 
the triplet populations increased. When the triplet to singlet rates increased 
the triplet populations decreased as would be expected. The He I G.S, population 
is very sensitive to the rate from the O.S. to the metastable 2 level ind 
2 level). The 1 level population was decreased by about 20 when 
n(l - 2 ^S) was increased by 100. At the same time the He II G.S. pupulation 

was increased. An increased rate from the triplets to singlets results in 

1 2 
increased 1 S 'pbpulation and a corresponding decrease in 1 S population. 

Ihe populations of the upper levels of the 58^ and 537 JL lines are very sensitive 

1 1 

to the rate from the 1 3 to 2 S level increasing as this rate Increases. 
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A code which we shall designate as code 1 has been developed to solve 
the steady state population equations given in Section II. B. This code is 
described in Section V. We discuss in this section some sample solutions of 
the steady state population equations, 

^e gas Is assumed Irradiated over 2n ster by a BlacKbody spectrum 
at 6000°K representing the photospherlc radiation field. The gas Is assumed 
to be optically thin for all lines and contlnuua unless otherwise specified. 

1. Ionization Equilibrium 

In this section we Illustrate results for the Ionization equilibrium. 
Hg.n-4 shows the ratio n^/ug^ with T^ - 30,000°K at various values of n . 

For an atomic model with only levels 1 and 20 the ratio n^/ng^ is approximately 
constant with n^ since both the colllslonal Ionization and reconib Inatlon rates 
are proportional to n^. The optically thin solution for two levels is shown. 
The results are generally within a factor of 2. For very low n^, all of the 
levels of He I except 5, 11 and 19 represent additional paths frXra level 1 
to level 20. This Is so because the photolonlzatlon rates from these levels 
exceed almost all colllslonal rates, n^/ng^ Is hence lower at n^ = 10^° than 
the two level solution. For somewhat higher electron densities the colllslonal 
rates between the singlets and triplets exceed the photolonlzatlon rates. 

Hence, the recombination to Intermediate levels results In conversion to the 
5, 11 or 19 level and thence to the 1 level by spontaneous emission. The 
rate from 20 to 1 is now enhanced sjhd n^/ng^ exceeds the two level solution. 

At still higher electron density the colllslonal ionization rates begin to be- 
come greater than either the spontaneous emission or colllsional deexcltatlon 
rates. Now all Intermediate levels represent paths from level 1 to 20 and 

“l^“20 lower than the two level solution. This is seen to be the case 

at n « 10^^. 
e 


We turn to the solution when the resonance lines are optically thick. 
The result for detailed balance (db) In the first resonance line is shown 
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by the curve labeled = 0 (p is the net radiative brackett). Also shown are 
results for db in the second and third resonance lines. Increasing the radiative 
excitation rate increases the effective ionization rate and the amount of n^Q by 
several orders of magnitude. We thus note the important result that the ioniza- 
tion equilibrium is strongly influenced by the optical thickness of the resonance 
lines. The LTE result is also shown and is labeled LTB. 

2 . Level Populations 

Many sample solutions have been obtained for the level populations or non- 
equilibrium parameters, b, where b => n/n 

' equll' 

We have completed a parametric study of the effect of optical thickness in the 
resonance lines and continuua upon the level populations. Level populations were 
obtained for various physically meaningful combinations' of resonance lines and 
continuua being optically thick or thin. That is to say each resonance line or 
continuua was assumed to have a net radiative brackett of either 0 (completely 
thlck)^or 1 (thin). The results cover the temperature range T - 10,000 to 
50,000 K and electron density range n^ - 10^° -10^^ cm"^. Due to space limita- 
tions we cannot present results for aU 30 levels of the model. Hence levels 

I, 3» 9, ao and 2? were chosen to Illustrate the results. Referring to section 

II. A we see that levels 1 and 20 are the ground states of He I and II, 

level 5 the upper level of the 584l\ line, level 9 the upper level of the D3 line 
and level 2? the iqjper level of the k6Q6k . line of He II. Figures II. 5-29 show 
results for these 5 levels for electron temperatures 10^, 2 x 10^, 

3 X 10 , 4 X 10 and 5 X 10 °K. In each figure the rat^io of the actual popu- 
lation to the equilibrium population, b, is given as a function of n . 

e 

At each temperature calculations were made for combinations of net radiative 
brackets corresponding to layers of varying total thickness. There are 5 
resonance lines in the model for both He I and He II. These 6 lines together with 
the Lyman continuua for each ion are allowed to become optically thick in 
our calculation. Thus, there are eight lines and continuua which can be 
optically thick or thin depending upon the physical thickness of the layer. 

Fac’h figure shows a completely optically thin solution labeled T and 8 
ou.rr solutions labeled 1 - 8 which correspond to the combinations of optical 
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Figure II -16 b. 
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Figure 11-22 
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thicknesses given in Table II. 3. On this table (l) refers to the 584^, line, (2) 
refers to the line, (3) refers to the 522 A line and (4) refers to 

the Lyman continuum of He I. Similarly the numbers 5-8 refer to the 304, 

'.‘'/j, lines and Lyman continuum of He II. When a number appears in 

Table II. 3 the corresponding transition has been assumed optically thick 
(i.e. in radiative detailed balance) in obtaining rate coefficients for the 
level population solutions. Transitions not appearing in the table are 
assumed optically thin. It is noted that progressing from case 1 to 8 corres- 
ponds in general to the layer becoming thicker. Case 8 always corresponds to 
all 8 lines and continuua becoming optically thick. For example, with 

“ 30,000°K case 4 represents the NRB = 0 in the first 3 resonance lines of 
He I and the first resonance line of He II. 

We shall discuss briefly the solutions presented In p'igures 11.5-30. For 
4o 

T^ = 10 Kraost of the He is He I. The He I ground state population is not 

affected by optical thickness except slightly at n = 10^^. Similarily Figure II. 6 

£> ® 

shows that once the $84 A line becomes thick the radiation field in other lines 

and continuua do not affect it. The D3 line upper level population shown in Pi- 

& 

gure II. 7 is a strong function of the 584 A line and He I Lyman continuum optical 

o 

depths and to a lesser extent dependent on the 537 and 522 A lines. Increasing 
the singlet populations hence increases the D5 emission. The optical thicknesses 
of the He II lines and continuua do not affect the D3 line. From Figure II. 0 ;i (, 
sooi that the He II ground state exhibits a similar behavior. At 20,000°K the 
He I and He II densities are very roughly the same. The thin solution shows 
that He I is greatly overpopulated and He II slightly underpopulated with respect 
to LTE, Increasing the optical depth in the He I resonance lines and continuua 
now increases the ionization by absorption of photospheric radiation from upper 
resonance line levels and thus increases the He II density and decreases the 
He I density. The D3 line upper level population increases with increasing 
population of the singlet levels but decreases when the He I Lyman continuum 
becomes thick duo to decreasing total He I density. As before the He II resonance 
line and continuum radiation fields do not affect the He I level populations. 

At 30 , 000 ‘^’K both He I and II ground states arc overpopula1;ed at the oeponso 
of He III. The He II ground state population is sensitive only l-o tfu; Lyman 
coni.inuum of lie II while the He I population depends upon both continuua and thie 
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Table Il-J 

OPnCALLY THICK UKSS CHOSEN 


Case Number 

123 A b 678 



He I resonance lines. The D3 line Is sensitive only to the optical thickness In 

the two I<jnnan contlnuua. Hie level 27 population depends upon the He II resonance 

lines and contlnuua only. At T ■ A0|00CPK all the levels have a strong dependence 

o 

on the optical thickness of the He II I^rman continuum. When this continuum Is thick 
the population of He IH Is greatly increased while those of Ife I and II are decreased; 
then decreasing all of the level populations. For levels 9 c^nd 20 this continuum 
has the largest Influence. Levels 1 and 5 are also dependent on the He I resonance 
line and continuum optical depths. At 5C>00CP K the level populations have the 
saim.! beJuivlor as at ^K);00 cPk. 
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j. Optical Thickness and Line Intensity 

Figures II.5“^-9 have illustrated the effect of optical thickness in 
ihe rr:r;onarice lines and continuua on certain level populations. Wo now sof'k 
t.o dr-tormirie t-?io approximate optical depths of a flare Fur purj^GW’S 

<)V i Llusi ratiiui we* consider a 1000 km thick layer;. For ea(‘h of Uk; sarrit* rmuic 
/^i veii in Ta'bJ^* II, > w*- ?iave obiaincd the line: ct'ritc'r optii^al thirkner.r. ui l.he 
O' V; ^ , D o l-O^l'^OA IJ tujs and at the liirv-shoid u:r l.iif' Lyiti/ui (‘ont.I.riuua, of 

li<- ;i arjd II, Of>l-|j:al depths art: not. shovrri for tiie ^LIin.'s. 

fi 

ij';eaus<- llr‘.y ai-e- always simply cons/tant fraction cd' llu.' 'jiVi and 

11 m: ofjticial Ui.LckneruiC?s , namely 'r( 537 ) " -27 t( 922 ) - .llT(^8t), 

t(256) - ,19 t( 30 ^ 0 , and - ,070 t( 304 ). 

TLe optical thicknesses are given in Tables 11,4 - II. 8 The He 
density is assunied to be one tenth of the electron density. The number 
rollowing cacti entry is the power of 10 hy which t>ie erilay is maltipli(-:d . 

W^‘ nolc t.hal bbt* 10830 and D 3 3 ,ines can become thick fur high elci*tr‘un den- 
rUU.er, c*verj al. 90/X)0^K, These .Lines do nut be'e<;me Uiif-k at. 10 , 0 (.)(/^K. TIku'c 
ar<- un.ny eas'-s i.rj whle.fi a nu/nber «.>f .Lines and eoritinuua are (^pkiealLy tlijck. 
This doeTJ noL mean >r.»wever l.hai. simultaneous transport o'quationr. niusl. i)e solv;-d 
l'(jr Uif.Gw; liner, and cuntlnuua. Which, line and eont.lnuuifi radiation i'if.'lds mnrl 
l)e (d)Lalned mulLaneously depends upon the level population being sui^/lil as 
weJJ, as the temperature^ density and layer optical thickness. To dete.nnine the 
c-ffec‘i. of vai"ious lines and continuua on levels 1, 9 ? 20 and 27 we can rei’er 

(,o .b‘:i./njrf‘S II. ’> * Py. For exam]>le, IVible II. 9 shows Lhr approximatr.' maximum 
pe rerai f/i/-*/' r rrevr, cneounterc'ti in t>ie Dp line uppe^r levc'l ])Of>ulal.i('in by solving 
only Tor i.fie lints: n.nd (‘onlinnua /^iven in eac:h box. Tiio el'I.Va’t on :n\v U'Vt'i 
'•an In: fdd.aJ rie'd Vtrm (lie oui.piil. from Code 1 . .figures II. !3 - ^'9 il.Lui; l.i-a(.e I lie 
< rrer:l. t)i' varying Uie NRL bel-ween 0 and 1. In firacticc: when ilio 'l.aye?’ Iiecornes 
eplrirstjly l.h.iek j i- dot's no(, ijsu;j..lly (lecome so t]i:iek (-hat Llie .NRP> 0 , Jn ro’dr-r 
\‘i\r l.iie l r’s.d i a (-:i fjM Pie.Ld in become sa(.u raid'd wc' requ i. r‘<* 

> ie Tl)"^ , 
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Table II -U 

OPTICAL THICKNESS - 1000 KM LAYER 






T 0 
e 

10,000 

“K 











Case 





TRANSITION 

n 

& 

T 

1 

2 

3 

4 

5 

6 

7 

8 


0 

1—1 

0 

H 










584 


3.7+4 

= 

= 

B 

1.3+4 

B 

B 

= 

3 

504 


7.6 


B 

J-S 

4.1 

B 

B 

3 

3 

5876 


2 . 8-8 

1 . 1-5 

2 . 1-5 

2.7-5 

8 . 0-4 

» 

» 

3 

3 

10850 


2.6-7 

1 . 0-4 

1 . 9-4 

2.5-4 

7.4-3 

a 

B 

B 

B 

304 


1.5-2 

9.6 



8 . 6+2 

B 

B 

S 

S 

228 


1 . 4-6 

9.3-4 

1 . 8-3 

2 • 

8.3-2 

B 

S 

3 

3 

584 


3.7+5 

= 

B 

B 

M 

B 

s 

a 

3 

5 o 4 


76 

b; 

EB 

0 

B 

B 

a 

a 

3 

5876 


1.3-6 

1.7-4 

2 . 2-4 

3.7-4 

6.7-3 


tm 

3 

3 

10830 


1.2-5 

1 . 2-5 

2 . 1-3 

3 . 4-3 

6 . 2-2 

B 

B 

3 

S 

5 o 4 


.046 




1 . 6+5 

B 

B 

3 

3 

228 

10 ^ 

4 . 4-6 

1 . 0-5 

1.9-3 

3.1-3 

.15 

= 

B 

= 

3 

584 


3.7+6 



= 

= 

B 

B 

s 

S 

5o4 


7.6+2 

S3 

B 

B 


B 

B 

3 

3 

5876 


1.5-5 

1 . 1-3 

2 . 4-3 

4.7-3 

1.5-2 

B 

a 

B 

3 

10830 


1 . 2-4 

8 . 6-5 

2 . 0-2 

3.7-2 

0.12 

BS 

s 

B 

B 

304 


0.10 




1.3+3 

B 

a 

a 

3 

228 

10^5 

9.3-6 

1 . 0-3 

2 . 8-3 

6 . 0-3 

1 . 2-1 

B 

B 

3 

a 

584 


3.7+7 

s 

B 

c 

B 

B 

B 

a 

» 

VJl 

0 


7.6+5 

= 

a 

B 

B 

B 

a 

a 

B 

5876 


2 . 1-4 

1.0-2 

3.1-2 

4 . 1-2 

6 . 2-2 

a 

a 

3 

3 

10850 


1 . 4-3 

4.7-2 

0.13 

0.17 

0.26 

B 

3 

B 

a 

3o4 


0.56 


1.5+22 

2.5+2 

1 . 1+3 

=> 

a 

3 

a 

228 


5.4-5 

2.7-3 

1.3-2 

2 . 4-2 

1 . 0-1 

B 

a 

3 

a 

58)1 


3.718 



B 


B 

3 

r; 

3 

504 


7.6+4 

ss 

B 

B 

BS 

B 

3 

3 

s 

5876 


5.6-3 

0.15 

0.24 

0.26 

0.29 

B 

3 

a 

a 

10830 


1 . 4-2 

0.37 

0.54 

0.59 

0.66 

B 

3 

3 

a 

304 



1.7+2 

5.0+2 

6 . 1+2 

9.3+2 

a 

B 

t= 

a 

228 


4 . 5-4 

1 . 6-2 

4 . 8-2 

5.9-2 

9 . 0-2 

a 

3 

B 

a 
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T = 2 X 10^ ® K 


Case 


ANSITION 

n 

p 

T 

1 

2 

5 

4 

5 

6 

7 

8 













10 










584 


1.7+5 

5.4 

0.92 

0.52 

= 

= 

= 

7.0-5 

— 

504 


0.56 

7.2-4 

1.9-4 

1 . 1-4 

= 

= , 

= 

!• 5“6 

= 

5876 


1.7-5 

1.8-5 

1.7-5 

= 

= 

- 

= 

9 . 6-4 

ss 

10850 


1 . 6-2 

1.7-2 

1 . 6-2 

= 

= 

= 

= 

9.5-5 

s= 

504 


1.5+5 

2.8+5 

= 

= 


= . .. 


2.0+5 

= 

228 


0.10 

0.19 

r= 


S5 


= 

0.15 



10 










584 


2 . 0+4 

2.2+2 

61. 

26. 

- 

= 

= 

1.0 


504 


5.7 

6.5-2 

1.7-2 

7.4-5 



= 

2.9-4 


5876 


4.5-2 

0.11 

= 

= 


- 

= 

5.5-2 

— 

10850 


0.58 

1.0 

0.86 

0.95 

• = 

= 


2.9-4 

— 

504 


5 - 1+5 

2.1+4 

= 

= 


r: 

= 

= 

— 

228 


0.45. 

1.8 

= 


- 






10 










584 


2.5+5 

1.0+4 

2.0+5 

5.6+2 

= 


= 

100 . 


504 


65. 

2.8 

0.56 

0.16 


= . 

- 

2 . 8-2 

= 

5876 


0.90 

•4.5 

4.0 

5.1 

= 

= 


1.0 

= 

10850 


5.5 

24 . 

22 . 

17. 



= 


“ 

5o4 


5.5+4 

2 . 0+5 

= 



= 

= 

— 


228 

10^5 

2.7 

16. 








584 


2 . 0+6 

1 . 2+5 

1.7+4 

8.5+5 

= • 

= 

= 

5.2+5 

~ 

504 


6 . 1+2 

57. 

5.2 

2.6 


= 

= 

1.6 

IS 

5876 


18. 

65. 

48 . 

58. 

- 


= 

28. 

= 

30850 


57. 

1.5+2 

90. 

69. 


= 

:= 

52. 


504 


4 . 0+5 

1.9+6 

= 

- 

= 

= 

— 

~ 

— 

228 

i 4 

56 . 

170. 

:= 


— 

= 





3 0 










'384 


2. OH 7 

1 .2h6 

5.8+5 

5 . 0 H -5 

r- 



2.6+5 


'3<)4 



V3<'. 

1 10. 

87 . 

: ‘ 



75 . 

*T 




7. 

'..4 (2 

H 2 



-T 

4 . 8 +;.’ 

.-t: 




; ' , ' 1 . * 

•-.812 

',.4 K.' 



- 

5 . 1+2 

r: 

.v'4 


I.U) 

1 .'MY 


- 

= 


= 

== 

= 

2: '8 


V/h 

I.. 7+5 

= 

= 

= 

= 

= 

= 

= 
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Table IZ-6 

OPTICAL THICKNESS - 1000 KM LAYER 

k « 

T - 3 X 10 ® K 
e 


Case 


TRANSITION 

n 

e 

T 

1 

2 

3 

4 

5 

6 

7 

8 

584 

10^° 

23. 

a 

22. 

21 . 

2.5-2 

5.2-5 

2.7-3 

1.7-5 

1 . 6^5 

504 


8 . 4-5 

B 

a 

8 . 2-3 

1 . 0-6 

2 . 0-7 

1 . 0-7 

5.9-6 

6.0-11 

5876 


1.3-3 

m 

B 

1.2-5 

' j . j-k 

7.0-4 

8 

4.3-6 

5.2-6 

10830 


1 . 2-2 

a 

B 

1 . 1-2 

7.2-3 

6.5-3 

6 . 4-3 

4 . 0-5 

3.0-5 

504 


1 . 6+3 

at 

a 

= 

B 

a 

a 

10 . 

B 

228 


0.18 

B 

a 

0.17 


3 

a 

1 . 1-3 

« 

584 


6.5+2 

e 

B 

a 

2.0 

0.43 

0.15 

8.8-5 

2 . 6 - 4 • 

504 


0.19 

a 

a 

B 

7.2-4 

1 . 5 -*^ 

5.0-5 

3.2-6 

1.0-7 

5876 


6 . 8-2 

a 

a 

a 

5.2-2 

4.5-2 

= 

2.9-3 

8.2-5 

10830 


0.57 

a 

a 

a 

0.43 

0.57 

a 

2 . 4-2 

7.0-5 

504 


1 . 6+4 

a 

a 

B 

a 

s 

a 

1 . 0+5 

B 

228 


1.7 

a 

a 

a 

a 

a 

B 

0.10 

B 

584 


7.3+3 

a 

. 

a 

62 . 

9.6 

2.5 

1.0 

0.13 

504 


2.6 

m 

B 

IS 

.22 

5 .i *-3 

8.9-4 

5.6-4 

4 . 6-5 

5876 


1.5 

a 

a 

» 

1.4 

B 

1.1 

0.42 

0.15 

10830 


6.8 

a 

B 

a 

CVJ 

• 

VO 

5.9 

^.7 

1.8 

0.65 

304 


1 . 6+5 

a 

a 

B 

a 

B 

a 

6 . 4+4 

a 

228 


18. 

a 

a 

B 

» 

a 

- 

7.0 

m 

584 

10^5 

4 . 8+4 


a 

a 

6.5+2 

76. 

58. 

32. 

17. 

VJl 

0 


17. 

a 

B 

= 

0.23 

2.7-2 

1 . 4-2 

1 . 2-2 

6 . 4-3 

5876 


22 . 

= 

- 

- 

18. 

15. 

12 . 

10 . 

7.3 

10830 


52 . 

- 

- 

m 

' 25 . 

19. 

15. 

15. 

10 . 

5 o 4 


1 . 6+6 

- 

m 

a 

- 

- 

a 

1 , 4+6 

- 

228 


1 . 8+2 

a 

m 


m 

. 

B 

1.5+2 

a 

584 

10’“* 

3.1+5 

3.3+5 

5.7+5 

. 

5.9+3 

1.8+5 

1.4+5 


1.1+3 

5 o 4 


1 . 2+2 

\ c 

1.5+2 

m 

2.3 

0.69 

0.48 

- 

0.38 

5876 


2.0+2 

= 

2 . 4+2 

m 

1 . 9+2 

1.4+2 

1.5+2 

B 

1 . 2+2 

10830 


2 . 6+2 

a 

2 . 1+2 

m 

1.7+2 

1.2+2 

1 . 1+2 

B 

1 . 0+2 

3 o 4 


1.5+7 

a 

B 

m 

1 . 6+7 

a 

a 

=» 

S 

228 


1 .& 4-5 

a 

a 

a 

B 

a 

a 

a 

= 

240 













Table II -7 

OPTICAL THECKMESS - 1000 KM LAYER 

T = 4 X 10^ ° K 
e 


ORIGINAL PAGE IS 
OF POOR QUALITY 



58U 


1.9 


1.5 

0.57 

1.9-4 

1 . 0-7 

1.7-8 

8.5-9 

5.5-11 

504 


7 . 6-4 

= 

4.5-4 

1.5-4 

7.6-8 

4 . 0-11 

6.8-12 

5.4-12 

2 . 2-14 

5876 


CO 

= 

6.4-4 

1 . 6-4 

1 . 6-8 

1 . 4-8 

1.2-8 

= 

6 . 6-9 

10850 


7.7-3 


5.9-5 

1.5-3 

7.7-7 

1 . 3-7 

1.0-7 


6 . 0-8 

504 


1.5+3 

= 

1.2+5 

2.9+2 

2.9-2 

c= 

= 

= 

= 

?. 2 Q 

10 ^^ 

0.18 

= 

0.11 

.035 

5.5-6 

= 


■ 


584 


4 o. 

58. 

35. 

15. 

6.7-3 

4.5-5 

8.0-6 

2 , 7-7 

4.7-9 

5o4 


1.7-2 


1.5-2 

7.0-3 

2 . 8-6 

1 . 8-8 

3.4-9 

1 . 1-10 

2 . 0-12 

5876 


0.52 

0.51 

0.27 

0.13 

5.6-6 

= 

5.8-6 

= 

1 . 2-6 

ID850 


0.26 

0.25 

0.21 

0.11 

4.5-5 


3.9-5 

= 

1.3-5 

504 


1 . 4+4 

= 

1.5+4 

7.6+5 

2.8 

= 

= 

= 

= 

228 

10 ^ 

1.8 

= 

1.6 

0.81 

5.6-4 

“ 




584 


4 . 2+2 ■ 

4 , 1+2 

3.7+2 

2.5+2 

0.82 

8 . 1-5 

1.1-3 

2.5-4 

3.5-5 

504 


0.17 

= 

0.15 

0.11 

5.4-4 

5-^-6 

4.7-7 

1 . 2-7 

1 . 4-8 

5876 


0.69 

0.67 

0.61 

o. 4 i 

1 . 4-5 

= 

- 

1 . 1-5 

4 . 0-4 

:i.o850 


2.6 

2.5 

2.5 

1.5 

5.0-5 

= 

4.9-3 

4 . 0-5 

1 . 4-3 

5o4 


1.5+5 

1 . 4+5 

1.3+5 

8.7+4 

2.9+2 

= 

= 

= 

= 

228 

10^5 

18. 

17. 

16. 

11. 

5.6-2 

V 




584 


2.7+3 

2.5+3 

2.3+3 

1 . 6+5 


0.80 

0.10 

4.5-2 

2 . 0-2 

504 


1.0 

0.92 

0.84 

0.68 

2 , 2-2 

2 . 1-4 

5.6-5 

1.8-5 

7.7-6 

5876 


9.1 

8.7 

8.0 

5.6 

0.19 

0.16 

0.15 

0.10 

8 . 0-2 

10850 


11 - 

10 . 

9.2 

6.5 

0.22 

0.19 

o.i 4 

0.12 

8.5-2 

504 


1,5+6 

1 . 4+6 

1 . 2+6 

8.8+5 

5.0+4 

= 

= 

= 

- 

;.'28 

! h 
10 

180. 

170. 

160, 

100 . 

5.6 

— 

— 

— 

' “ 

584 


2 . 1+4 

2 . 0+4 

1.7+4 

1 . 2+4 

3 - 1+3 




9.0 

5 ( )4 


8.4 

y .8 

6.4 

3.0 

0.77 


1 

0 

5.2-5 

2 . 2-5 

',.".70 


'I'l. 


8 V 

58. 

15. 


10 . 

0.0 

8.5 

lo8-'iO 


Y(-. 

72. 

W'. 

45. 

11 . 

JO. 

7-0 

6.0 

6 , 0 

•'io 4 


1 .4 ly 


.1 . 2 ty 

').o +6 

2 . 5-16 

- 


- 

- 



1 .8t 5 

1 . Y-i ') 

,1 . 4 1 5 

.1 .;m 5 

5 - 0 '2 

- 

=; 

- 



Tacie Il-a 

OPTICAL THICKIIBS8 - XOOO KM lAYER 

k 

T = 5 X 10 ® K 
e 

Ik 


Case 


TRANSITION 


T 

1 

2 

5 

— 

5 

6 

7.._ 

8 


,.10 

10 










584 


0.35 

0.21 

4.2-2 

3.5-3 

2.1-7 

• 

- 

- 

- 

504 


1.5-^^ 

9.0-5 

1.8-5 

1.5-6 

9.0-11 

- 

- 

- 

• 

5876 


5.4-4 

3.4-4 

6.5-5 

5.5-6 

5.2-10 

- 



- 

10850 


4.9-3 

3.0-3 

6.0-4 

5.0-5 

5.0-9 

- 


- 

- 

304 


1.5+5 

7.8+2 

1.6+2 

15. 

7.8-4 

St 

= 

X 

- 

228 


0.18 

0.11 

2.2-2 

1.8-3 

1.1-7 

3 

= 


8 

CO 

ITS 


6.6 

4.0 

1.3 

0.26 

4.0-5 

- 

- 

- 

- 

504 


5.0-3 

1.8-3 

5.9-4 

1.2-4 

1.8-8 

- 

- 

- 

mm 

5876 


1.8-2 

1.1-2 

3.6-3 

7.2-4 

1.1-7 

m 

- 


- 

10830 


o.i4 

8.4-2 

2.8-2 

5.6-3 

8.4-7 

- 

- 

m 

mm 

504 


1.3+4 

7.8+3 

2.6+3 

5.2+2 

7.8-2 

n 

- 

n 

s 

228 

10^ 

1.8 

1.1 

0.36 

7.2-2 

1.1-5 

m 

■ 

m 

• 

584 


62. 

37. 

12. 

5.1 

3.7-3 

- 

- 

- 


504 


5.0-2 

1.8-2 

5.8-5 

1.5-3 

1.8-6 

- 

- 

- 

- 

5876 


0.37 

0.22 

7.4-2 

1.8.2 

2.2-5 



- 

- 

10850 


1.2 

0.75 

0.24 

6.0-2 

7.2-5 

- 

- 

- 


3 o 4 


1.5+5 

7.8+4 

2.6+4 

6.5+5 

7.8 

D 

- 

0 


228 

10« 

18. 

11. 

3.6 

0.90 

1.1-3 

“ 

m 

8 

" 

584 


4.1+2 

2.4+2 

82. 

20. 

0.25 

4.0-3 

4.5-4 

2.3-4 

9.0-5 

504 


0.18 

0.11 

3.6-2 

8.8-3 

1.1-4 

1.8-6 

2.0-7 

1.0-7 

4.0-8 

5876 


4.8 

2.9 

0.96 

0.24 

2.9-3 

2.6-3 

2.2-3 

1.7-3 

1.3-5 

10830 


4.9 

2.9 

1.0 

0.25 

2.9-3 

2.6-3 

2.1-3 

1.7-3 

1.3-3 

304 


1.3+6 

7.&^5 

2.6+5 

6.5+4 

7.&*2 

- 

m 

- 


228 

10^“ 

180. 

110. 

36. 

6.8 

0.11 

■ 

■ 

m 

■m 

584 


3.<H3 

1.^3 

5.0*2 160. 

19. 

0.35 

0.11 

8.5-2 

5.9-2 

504 


1.4 

0.84 

0.22 

7.0-2 

8.9-3 

1.6-4 

5.1-5 

4.0-5 

2.8-5 

5876 


47. 

30. 

7.8 

2.4 

0.29 

0.25 

0.20 

8 


10830 


35. 

21. 

5.6 

1.7 

0.21 

0.18 

0.l4 

8 

= 

^o4 


\ . M 7 


2.1-»6 

6.5+5 

7.8+4 

= 

a 

X 

s 



l S 

1.1+3 


9.0+1 

11. 

= 

8 

8 

s 
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Table II -9 ERROR IH D5 Line Intensity 
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TABLE II. 10 


CHABAjCTTERISTIC VALUES OF e, Tl, C 


V°K) 


e 


1 

Line 

(i) 

Case 

# 

10*^ 

io “ 

4.1-9 

2.3-4 , 

2.9-18 

584 

T 



II 

It 

5.1-14 

ft 

3 



1.2-8 

1.2-8 

1.8-28 

304 

T 



II 

ft 

tt 

fl 

3 


lO^" 

4.0-5 

8.0-3 

4.0-14 

584 

T 



II 

II 

8.3-13 

tl 

3 



1.2-4 

4.8-5 

1.8-24 

304 

T 



fl 

ft 

II 

ft 

3 

4 

2 X 10 

10^ 

1.1-8 

2.3-4 

6.1-13 

584 

T 



II 

II 

3.8-3 

tl 

3 



8.2-9 

1.2-8 

2.5-18 

304 

T 



n 

ft 

tt 

ft 

3 


xo^ 

1.1-4 

2.2-2 

6.4-9 

584 

T 



n 

tl 

4.9-7 

tl 

3 



8.2-5 

5.8-5 

2.5-14 

304 

T 



n 

II 

tl 

II 

3 

3 X 10^ 

io “ 

1.9-8 

2.3-4 

4.3-11 

584 

T 



ii 

II 

tl 

n 

3 



6.7-9 

1.4-8 

6 . 1-15 

304 

T 



n 

II 

4.2-9 

11 

3 


lo"* 

6.7-5 

1.2-4 

5.8-11 

584 

T 

■ 


II 

II 

1.4-8 

II 

3 



1.9-4 

3.4-2 

4.0-7 

304 

T 



II 

II 

It 

II 

3 

4 X 10** 

io“ 

2.8-8 , 

2.3-4 

3.9-10 

584 

T 



If 

2.1-4 

2.9-10 

It 

3 



5.8-9 

1. 6-8 

3.1-13 

304 

T 



ir 

II 

^•7r7 

n 

3 
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TABLE II. 10 (Continued) 


Case 

# 


Te(°K) 


(om-3) 


10 


14 


5 X 10 


10 


10 


10 


14 


Line 

T) C (A) 


2.8-4 

4.5-2 

3.4-6 

584 

T 

If 

II 

It 

II 

3 

5.8-5 

1.9-4 

2.8-9 

304 

T 

fl 

II 

1.9-6 

tl 

3 

3.8-8 

2.3-4 

1.6-9 

584 

T 

<1 

2,0-4 

9.0-10 

fl 

3 

5.2-9 

1.9-8 

3.4-12 

304 

T 

tt 

tl 

8.1-6 

It 

3 

3.8-4 

2.7-4 

3.0-8 

304 

T 

If 

5.8-2 

If 

tl 

3 

5.2-5 

2.7-4 

3.0-8 

304 

T 

Tl 

tl 

2.7-5 

II 

3 
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for our assumed Doppler line profile. (For the continuum we require t > 
where C Is the ratio of colli slonal to radiative recoaiblnatlon.) 

Some representative values of the parameters e and T\ (and t) For the 

^84 and 304A lines are given In Table 11.10 . Values were obtained for each 

temperature at n « 10^*^ and 10^ for the optically thin case and case 3 
6 

given In Table II. 3. We can see from the optical depth tables that under many 
conditions the layer is effectively optically thin, l.e. 

1< T < (s + 

In this case the total energy emitted In the line from the surface of a layer 
of optical thickness t^ Is given approximately by 

E~2/W (s + !♦**) t^ 

where = L/B. 

The values of T) and t vary with optical depth and an Integration 
over depth Is required to obtain the total line intensity from the layer. 

Using the optical depth tables and estimates of 1) and L from code 1 
we have obtained approximate solutions for some line Intensities. The lines 
chosen are the ^37« 304^ Unas and the D3 11m. Figures II.3P-3i^ give 
results as a function of T. for various n. while Figures U.33-39 show the 
variation with n^ for various T^. We note a considerably different variation 
with T and n for the D3, 584 suad 304^ line emission. These different 
variations give some confidence that simultaneous emission meMurements 
of the three lines could yield unique temperatures and densities from the 
emitting region. 


i 
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Figure II -50 


TOTAL LINE EMISSION 







Figure 11-51 


TOTAL UM immiTi 


n - 10^^ 

0 
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Figure 11-33 TOTAL LINE EMISSION 
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k. Population Rates 

It is very instructive to look at the various processes involved 

in populating and depopulating the upper radiating levels. In other 

words, to determine the exact processes contributing significantly to 

the terms T| and L in the transport equation. One way to do this is to 

assume a 3 or 4 level atom and evaluate T) and /. from the analytic 

"tti 

solution choosing different levels for the 3 or 4 levels. This 
would give us an approximate answer. Another way is to use the 
complete solution of the statistical equilibrium equations for the 
populations to evaluate the rates directly into the upper and lower 
levels of the lines in question. 

Solutions of the full SSS equations were used in obtaining the 

relative rate processes shown in Figs. 11-40 and 11-41. Fig. 11-40 

Illustrates the processes populating the 564A line upper level at 

T = 40.000^K for various electron densities. Four processes always 
e ^ 

enter significantly. The largest of these is photoexcitation from 
the 2^S level by absorption of photoapheric radiation. Direct colli- 
sional excitation from the ground state is next, followed by colli sional 
excitation from the 2^S and radiative decay from the 3^D. Hence, both 
the photospherie radiation and coupling to the triplet levels are impor- 
tant. 


Fig. 11-41 shows relative processes populating the D3 line upper 

level for the same electron temperature and densities. A much stronger 

dependence on n^ is noted. Fhotoabsorption of photospherie radiation In 

the D3 line Itself Is the dominant mechanism at low n^. At high n^, 

this process is small and the collislonal rate from the 3^ is dominant 

3 3 

with collislonal excitation from the 3 P and 2 S levels also being 
significant. Thus, the triplet - singlet interaction is again very 
Important. 


I 
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Figure 


BEIATIVE PROCESSES POPUIATING 2*P LEVEL OP He I 
Te » 40,00 CPk 











III. LIME TRAMSPOHT SOLUTION 
A. Basic Equations 

The solution of the line transport equation that we have used in 
the code developed (Code 2) Is an Iterative solution of the Integral 
equation. The details of the method have "been developed by Avrett an d 
Loeser (I969) In a manner convenient for the simultaneous solution of both 
line and continuum transport equations. Prom Section II. D, we note that 
solution for the line and continuum radiation field simultaneously is 

necessary. We will only briefly summarize the method and equations In 
this report. 


In the statistically steady state the rate equation describing the 
population n^ of the state 1 is 


^11 - ‘■ly 111.1 

where is the total transition rate from 1 to J per second per particle 
in the 1 state. In general, P^^ = vhere and represent 

the radiative and collislonal transition rates respectively. We shall 
assume a Maxwellian distribution for the electrons and helium particles and 
since we also assume a known external radiation field, the transitions 
Involving the continuum can be represented by a single term In Eq. m. 1. We 
can characterize the system of linear equations m.i by a matrix whose co- 
efficients a^^j are equal to la representing matrix elements and 

co-factors thereof, we shall always let the first subscript refer to the 
row and the second to the column. Because of the definition of the transi- 
tion rate the subscripts of the P* s will be reversed when substituted for 
the matrix elements a. The general solution of Eq. m. i is 


n, = X P^ ; 
1 m ^ 


N 


m 


E P“<5 
J 


III. 2 
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where is the co-factor of the element a . = P. in the coefficient ma- 

mj jm 

trix represented by the i equations of type mi and N is the total number 

3 

of helium particles per cm . The matrix of coefficients has the property 
that the co-factor of all the elements in a column are equal, i.e. is 
independent of m. 


When the medium becomes optically non-thin for certain frequencies 
the radiation field producing Internal excitation and ionization for these 
frequencies is no longer merely the external radiation field but is partly 
dependent on the internal properties of the gas and must be determined from 
the radiative -transfer equation 


dl 



where cos u represents the angle between the direction of propagation and 
the outward normal z and t = (*k dz, S is the source function, k is 
the linear absorption coefficient and is the specific intensity of the 
radiation. In USE however in the non-LTE case must be specified 

in tenns of microscopic processes. In terms of such processes the trans- 
fer equation governing the spectral line between upper level u and lower 
level I may be written as 


- 43tu 


dl 

V 

dz 




I 

V 


III. 4 


n 




Where represents the continuum absorption coefficient at frequency v 

and S is the continuum source function. B , B and A ^ are the Einstein 
C XU ux ux 

transition probabilities for absorption, stimulated emission and spontaneous 
emission. and are the normalized emission, stimulated emission 

and absorption coefficients within the line defined such that 

CO OD a> 

I = f i|f^dv = r r du) = 1. 

o o 4rt o 
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The continuum abaorptlon coefficient Is generally very Rmai i compared with 
the line absorption coefficient near the line center and will be neglected 
in determining the source function within the line. Using the standard 
relations between the Einstein coefficients and assuming J = t = iit 
the source function becomes 

where g represents the statistical weight. The minus one term in the 
denominator represents stimulated emissions. 


In evaluating the radiative excitation rate R^j for transitions 
between bound levels the line radiation field enters as 

€0 

o 

where 

= sr I, 

1b the nieezi intensity ond du) represents the solid suigle* Jt Is thus con- 
venient to formulate the transfer equation In terms of J rather than I . 

V V 

It Is now convenient to separate those components involving the unknown 
radiation fields from the co-factors# This is done by expemdlng the 

determinant in terms of Its co-factors Thus 


a r p 

p'^ »= E P 

4u 


P + 
^lu^ * 



^ul<^' 


lu 


+ E P 

Vu^l 



III. 6 


Actually, may appear in many of the co-factors since the line u-1 nay 
fall in the ionization continuum of some other transition. The Influence 
well as the line radiation In generca on the bound-free radiative 
rates wl31. be neglected. Using the st an da r d relationship between the 
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Einstein coefficients, equation III. 6, and remembering that 
Q = Q and J J^$^dv the source function may be written as 


^ul ^ul pUl 


J J^i^dv + eB + 

ul 

1 + e + T) 


L 


III. 7 


where B is the Planck function 


^"^e^ " *^ul C 


lu 

ul 


e = 


"ul 

^ul 


2 hv' 


ul 


&, 


u 


L ■ "ul B /, '■ik*' 

A^Q k/uA 


kl 


t| ^ S p 

. ^ul , / vik^ 

A^iQ 


III, 8 


The terms entering the numerator of Eq. III. 7 each represent a method of 
populating the upper level from the lower level. The first term repre- 
sents direct radiative excitation, the second direct colllslonal excitation 
and the third any combination of radiative or colllslonal processes involv- 
ing one or more intermediate levels in going from the lower to the upper 
level. The denominator, on the other hand, consists of terms indicating 
transition paths from the upper to the lower level. All the terms are 
normalized with respect to A^^. The first term represents direct radiative 
de-excitation, the second direct colllslonal de-excitation and the third 
any Indirect process going from the upper to the lower state. 
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Equation III. 7 1 b solved using a discrete ordinate method for 
the frequency IntegretL In vhlch we assume 


f F(x)dx . III. 9 

o 


The coefficients are to be found from a given set of dimensionless frequency 
values Xj^, k = 1,2,..K; is the value of F at x *s Xj^. The solution 
for Is given In Section III.C. 


We shall assume t to be Gaussian, l.e. 


_ exp( -v^) . Av 

■ 00 ^^^ " 


D 


Avjj Is the Doppler half width given by 


_2 

c ^ M ' 


1/2 


III. 10 


III. 11 


where M Is the Helium particle mass. 


The source function Is obtained at N depth points (i = 1 — - N) 
within the assumed layer. The depth points are located at specified physical 
depths which do not change during the calculation. The values of are 
obtained from the matrix equation 


N 

E 

J=1 




III. 12 


where 1 and J refer to depth points. 


The coefficients are given by 





where 


ij 


1 , J “ 1 
0 , J / 1 


III. 13 
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are weighting functions, e^is a coupling parameter. 


C . are given by 
J 


C. = 
0 


1 " M'l c 

1+e". “i l+€^. '^ij 

1 1 j=l 




III. 14 


S . 


In these equations, is a coupling parameter. 


- (A) 

W., ' 
ij 


and 


w 


ij 


.S- r A w 

- /T, \ 


(A) 




"k "ijk 


?lkil 
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III. 16 


where k refers to a specific frequency, Aj^ are wei^ting functions, 

i, 

profile function is simply 


The 




ik 


= e 


The optical depth values T are given by 

IJv 


*^ik 




dT = k dz 

, L 


h.v 

^2 '"I ®12 
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The weighting coefficients represent an expansion of the mean 

intensity terms of the source function 




j=i 


ijk v' vj 
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Substituting the solution of the radiative transport equation for 


J we obtain 

V 


N 


t J E (|t-T I) S(t) dt * E W s(t.) 
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(A) 


Where the frequency subscript has been dropped* To evaluate W . ^ ^ we assume 

that is represented by linear segments between optical depth points n and 

n + 1 . S(t) in the interval n s t s T is given by 

n+x 


T ^ - t t - T 

S(t) a S (-Qii ;—) + S , ( ^ ) 

\+i - \ Vi “ \ 


III. 19 


Equations fcrW^j^^ are obtained by substituting Eq. ILI9 into Eq. II. I8. 

The so obtained depend only on the set of values chosen. Different 

expressions for s,re obtained for J > i ^ J < i, and J = 1 . ¥e give 

these expressions in order . 


A 


W depends only on the coefficient of S 
1 J n 




(t„,- t) 


^n-l ^ '•’n 


III . 20 
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For ri > 1, 


2 W. = [ 
in 




+Eo T -T. -E, T ,• 

3' n i' 3' n-1 


“ T -1 J 
n n-l'^ 




Ti(E2|t„-Ti|-E2|t„^1-tJ)+ e 


•e 




■" ®3l%+l-'ll-^3l'^n-'l 
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For n < i, 


(^i-"n-l> 


(■'n+l-^l> 


in C T 




“(^l“T ) -(T..-T ,) 

-f* *®3(VVl) -"3'^l-VV(V-n-l)- 
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T - T , 

n n-1 


'^n+l" 
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B. 


Evaluation of e^and 


C" , 


«' and B appearing in Eqs. III.IJ and are related to the source 

function parameters ,T)and t as follows: 


e'B = eB + C 

« ' as S + T1 


111.24 


where B is the Planck function. Ibe values of and b“ depend upon the 
levels Included in the model other than the upper and lower line levels. 

For two bound levels (l,2) and one continuum level (k) t* and are given 

by: 


^1 21 0 )^ 21 


C CD P 

c 'b® = (l - p) B + a - 2 ^ 


«2i 

2h 


where 


a 




•,2 

c 




^2 ^1 


-hVg^/kT 


For three bound levels and a continuum level 


III. 2 5 
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X - Y 


/ s 

f B - OY 


where 


1 - _ M ^L 

X = (c + P + P +P - iWi ^^1 

'^21 *^21 ^23 23 ^^2 


III. 27 


III. 28 
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^ ^ ^12 

21 


“12 


R. 


^2 




^12 ^52 ^52 


1 ^ 2 ^ = 


‘25 23 


M- — P ^P "^P *^P 

22 51 ^51 ^52 52 


«2 = ^3-^^15 
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For some applications it is convenient to solve for two line radiation 
fields simultaneously using only the foxir levels in Code 2. The two lines 
are then represented by the 3 ~ 1 and 2-1 transitions. The corresponding 
X and Y values for the 3 ■ 1 lines are 


1 _ _ M.o 

* "51 " '’52 * > ■ 


Y = 


!k. 


R 


A 


w 


1 ) 


where 


^11 ^21 ^21 ^23 ^23 


®i ^12 ^12 


III . 31 


III . 32 


/ s 

In the general case c and B are obtained from Code 1 at each 
optical depth. Code 1 can be run for arbitrary values of J for the lines 
or A for the continuum rates. 
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C. Evaluation of 

The evaluation of Aj^ follows that of Avrett and Loeser (I969) 
which is recommended for further details. is represented by 


K 




0 :£ X s; 




where 


f(x) = 1 , 0 s: X s 


when J = 1, and 


f(x) 


■! 


(1 - ~) (1 - y 0 <: X s 


0 , 


x^ s X Xj^ 


when J s 2,3,. . .K.ly Is an adjustable pareuueter. 


III. 33 


III. 34 


is given by 


where f 



is the Inverse 


of f(xj^)j and 




Xj^ f J ^ 1 

I Xj (1 - ^), J . 2,3 ...K. 


III. 35 


III. 36 
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D. 


Sample Solution - l\fo Linea 


The inputs to code 2 to begin the iterative solution are the 
n^/ni. Approximate starting values are obtained from code 1. We shall 
Illustrate the solution by discussing the populations n^n^ and 

rij^/n^^, where 1, 2 and 3 designate the lower level and two upper line 
levels for which the transport equations are solved and K the continuum 
level, and are obtained as described in Section III.B, n and 
the optical depths at the geometric depth points are calculated by the 
code from the given population ratios and For the initial run, 

the upward radiative rates are those given by the optically thin case. 

Sj values are calculated from Eq. III. 12. is found from the S and 
the new radiative excitation rate based on is used to obtain new 
values of n^/n^ through code 1 or a simpler 3 or 4 level solution. The 

new values of and J are used as input to code 2 (or 3) to obtain 

the other line intensity. 

We Illustrate the solution for two lines in Fig. III-l for T = 

2 X 10 °K and n^ = 10^^ cm’^. The Figure shows results for ^ 

and as a function of the number of Iterations. The Initial 

values of n^/n^ and n^/n^ are 10 and is 10^. The solutions 

have essentially converged after 2-3 iterations. Solutions at the surface 
of the layer are labeled -0; solutions at the center of the layer are 
labeled -C. The codes have been tested for a number of simultaneous line 
and continuum transfer problems. There have been no convergence difficul- 
ties provided the lines and contlnuua chosen for solution do significantly 
effect each other. 
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IV. CONTINUUM TRANSPORT SOLUTION 

A, Basic Equations 


The statistical equilibrium equation for one bound level (n^^) and 
one contiutm level (i^) Is 





[ 




J d V + C 1 
h.v *'v Ik-* 



-hVKiT 

e 



+ Jy) dv + J, 


IV. 1 


where the * refers to the LTE population at a given electron temperature and 
density. It Is convenient to express the population ratio n^/n^ in terms of 
the non-equilibrium parameter b^. 


IV. 2 

Ttie radiative transpoi*t equation for v ^ can be written 



where neglecting stimulated recombination (which is Justified for He for the 
temperatures considered) 


B 



-hv/kT 

e 


dT = k dz 
V V 


k 


V 



IV. 4 
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The values of at depth 1 are found by solving the set of equations 


N 

S M, 
J=1 


IJ 


IJ 


1, 2 N 


IV. 5 


where 


M 


IJ “IJ 


e® Ri 

®1 "l 


^ 0. J (< 1 


K 

Z 

kol 


A' W 


IJk 


(A-1) fk' 


-y. 6, 
e k J 


IV.6 


and 


"i 




IV. 7 


is the optical depth at the continuum threshold at depth 1, 0^ = h.v^/kT^, 
yj^ = k = 1, 2, ...K, a set of dimensionless frequency values and g|^ 

are values such that y^“^ represents the frequency dependence of the photo- 
ionlzatlon cross section. The optical depth is given by 

Sv I 


Ik g 


1 yi 




IV.8 


8f 

e. and e. are coupling parameters whose calculation is described In the next 

“ related to as follows, 


section. The values of w. 


w. 


IJk 


(A-1) 


V 


V 


IJk 

. ( 

IJk 


(A) 

'(A) 


- 1, S 


, J 1. 
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The mean intensity of the radiation field Is given by 


'Ik 


N 

Z 

J"1 


w. 


ijk 


(A) 


B* 

fjk_ 
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Q, "b 

B. Evaluation of _and_e^ 


for h^, 


and follow from the statistical equilibrium equation solution 




b. 


where 


1 + e 


00 a 
''l 
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IV. 12 


For two 


bound levels (1^2} and a continuum level kj and €.^ are 


'*’= (<ik* ) 


R 


Ik 


P + P 
^21 ^ *^^2k 


,a _ _1 (n , ("2 /”l ) *21 ("k/”2)'’k2 

' ' Ik P„. + P, 


R 


Ik 


21 "2k 


IV. 13 


For an N level model the corresponding equations are 


b 1 


N 


® ~ ^ (Cti, + 1,T £ 


R ''^Ik ' ^iTT *“ ^1 a 


a 1 


N 




^kA Q r* 
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Or b 

e and e can easily be obtained from program 1. 
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C. EvaLuatlQn^0f_^J^ 


! 

The coefficients ore defined such that 


P(y)dy . S a; 

1 kal 


IV. 15 


where is the value of F at y = yj^. Here, 1 <; y^^ s y^. With unity 
instead of zero as the lower integration limit, we detennine the 
coefficients as follows. 


Let f (y) =1, 1 s y s yj^, 


when j <= 1, and 


f (y)j - 


I 

( 0 , 


(1- (I y* y. 


IV. 16 


^ ^ 


when i = 2,3, ...K. Let f^“^ be the inverse of f'(yT,)j. 


Then 




IV. IT 


where 


6' 


1 5 


yjj”li j ■ 1 


IV, 18 


5 (y^-l)(l- •j); J • 2,3,...N. 


Y Is an adjustable paremeter. 
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V. CODES 

Cotig 1 ~ Solutloja of* 1 V steady State Pnpulflt.l^j^ 

EuuaUoaa 

1 . Main Program 

In our own Internal notation, this code is designated as P48. The 
main program reads f numbers, wavelengths (lambda), electron and radiation 
temperatures, the energy level model, recombination coefficients, con- 
verts energy differences from wave numbers to electron volts, and writes 
all the above before entering subroutine CONSTS. The subroutine CONSTS 
as well as other subroutines are described in the next section. The 
radiation temperature is that characterizing the external Blackbody 
radiation field Incident upon the layer. The indices on the f numbers 
refer to the ion and the upper and lower line levels respectively. 


On return to the main program, more quantities are read in and 
printed out, such as: the total line center optical depth of the layer, 

the dilution factor for the external radiation, the electron density, 
the indeces of the chosen thick line, the depth points at which 
solutions are to be obtained, a constant divisor for the elements of 
the main matrix to prevent overflow, Y values for thick lines, and 
Various control numbers. These quantities are entered here in order 
to be able to run a series of different solutions by changing any or 
of them without changing the electron or radiation temperatures 
and, thereby, having to recalculate the many reaction rates. 

After the above have been printed and the ratio of cal- 

culated the subroutine AMAT is entered to calculate the elements^and 
cofactors of the A matrix for the optically thin case. The rate equa- 
tion matrix represented by Eq. II-3 is printed out. The matrix Is 
30 X 30 corresponding to the number of energy levels. Selected cofactors 
are printed out. The accuracy of the solution can be checked by compari- 
son of cofactors with the same J values. 


277 



nie solution for the non-equlllbrlujn populations n Is obtained 
in MAT. The solutions are printed for the optlealUr thin case along 
with ME solutions for the same temperature and density. The optical 
depth at the center of each line Is calculated and printed out for a 
layer of optical thickness Tj^ In the line specified in the input. 

Some further quantities such as the physical thickness of the layer H 
are calculated and the suhroutlne ELIM Is entered to compute e, 11 and C 

The program now repeats the ca;.culatlon with the specified set 
of 1 walues read in. Next, a Do loop on depth la set up. In the list- 
ing given there Is no depth variable. In calling for one to obtain 
values of 1) and c for example, it would enter Into the calculation of 
the quantity rp, idilch Is used to alter the Y value for the thick line. 
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CQl 


2 . 


Description of Sub-Programs 


ubroutinc AMAtCKK) 

AM/\T (;a]-culate.s elements of the A matrix using quantities determined 

b;/ CONSiTj. Onc;e all the elements are found, they are divided by a constant. 

to prevent, overflow. Cofactors are calculated, summed, and used to determine 

I'Yom this the n. are found. Next the n . are found and the ratio n./n^ . 

1 eqi r t'qj 

Last, before returning to the main program, the optical depths are calculated 
and printed for each f value. 

KK = 0 , optically thin solution 

KK = 1 , optically thick solution for any Y ^ 0.0. 

I^inction COFACT (HR. NC. NE. P) 

COh'ACT finds the cofactor of matrix D for row NR, column NC. NE Ifj 
the number of rows and columns in D. 

Function Cg>LI. ('f. I, d. K. N) 

tiuiction C0LL calculates certain electron collisional excitation and 
ionization coefficients for special points. 

T = electron temperature (°K) 

I, J, K = indeces of coefficient 

N = 1 , excitation coefficient calculated 
N = 2 , ionization coefficient calculated 


Ou brout.ine CONC’.I.’C 


Initially, if the radiation tempo-'rature is different from the elfM:tron 

l.r.mperature, CONC'Kl rcad.s in new value.s of P.. .. caiculatod a1. 'I' . It then 

'J. J k r 

(tn I r*!! <M*t/ron coll Iniorial. ilonizati on and rocombJ naU on coon*;] c i < Tit;; , 

Uif radi.al>ivo :Lori;] znd/i t;n rut.o, Uio li. .. and B*. idir- ol(‘<?t;ron (‘oil :i M,i onu.l 

IJK. .I.K.J 

t'xo.i i on oc k'I’IM o.i ord., \.Yw Kin:H,e:in sponLancomi transition probabilit.y^ find 
l.ho I'li n;d.o i on Invinn.i l.*i on pr*i)babil:l ty multiplied by the Planck function. 

Tho Luibiic r;l pi, i ivroru to tlio ion, the second subscript the beginning level 
and the last uubecrlpt the ending level of the trajisition. 


2 


origin ALPAGBK 

OF POOR QUALH^ 
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Subroutine ELIM (RHO. AP) 


AP = Einstein spontaneous transition probability for the thick 
line (K7, K8, K9) 

RHO = hv5/c^, where v is calculated at (K7, K8, K9) 

This computes double cofactors for use in calculating eta and lota. 

It is called from the main program and uses the function EZ^. 

litinction me (IRA. ICA_. KRA^ KCA. BMAT. Z) 

IRA s row index of 1st row to be eliminated 

IGA s column index of 1st column to be eliminated 

KRA = row index of second row to be eliminated 

KCA « column index of second column to be eliminated 

RMAT «= contains reduced matrix 
Z «= sign of cofactor (Z = ± l), 

ELM2 eliminates two rows and two columns from the A matrix and places 
the reduced matrix in RMAT, 


3 
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B. 


Radiative Transfer Codes 


Two codes have been developed under this program. The first 
(code 2 ) represents a numerical solution of the line transport equation 
given in Section III. The second (code 3) represents a numerical solu- 
tion of the continuum transport equation. Listings of these codes appear 
in Appendices B and C. The simultaneous solution of several line and/or 
continuum transport eqs. with the population eqs. is accomplished by 
iteration. Initial values for population ratios are estimated from 
various solutions obtained from code 1 (see Section II). After obtain- 
ing the radiation intensities from codes 2 and 3 the new radiative 
rates are used as input to code 1 to obtain new population ratios to be 
used in the next iteration. 


Code 2 requires the specification of certain quantities involving 
and while code 3 requires 6^ and There are two methods provided 

for determining these quantities. One option is to calculate them in 
code 1 and simply read them in. (Code 1 is now set up to calculate the 
quantities T] and which are read in and are related in a simple manner -- 
see Eq. III. 24. The quantities e®’ and are easily related to the co- 

factors generated by code 1. Code 1 will generate arbitrary cofactors 
depending upon input parameters.) 


A more approximate method of calculating these quantities is pro- 
vided directly in the two codes. Code 2 will obtain the solution for 

» o 

(C and B'" from a 4 level model (3 bound levels and 1 continuum level). 
Code bus the option of solving for and from a 3 level model 

(2 bound levels and 1 continuum level). These approximate solutions can 
be very useful for some purposes • 


A symmetrical (about the center of the layer) geometrical grid 
(<'ni;t(.>d y, in thn program) is set up given the total geometric depth and 
miml'ci- .<r decades. Ciauviitly , the maximuni number of decades is 

set at giving a total number of points of 3l. This calls for solution 
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of a 3l X 31 matrix and may create underflow or overflow problems. 

The number of decades required for convergence must be tested for each 
problem. 


Optical depth points at each geometric depth are next calculated 
from an integration over geometric depth at a given frequency. 

Steps in frequency are now calculated using the quantities 

V- V 

X £ bound level cases and y = v/v^^ In the continuum. Step 

' Av^ 

sizes are input quantities and are constant in x and y. The maximum x 

is found by choosing the first value of x such that e"^ t ^ 1.0. 

3 inax * 

and the maximum y such that t /y i 1.0. 

max''' 

Certain quantities uniquely dependent on frequency are now calcu- 
lated, such as f(y), g(y), A(y), f’(x), g'(x), and A'(x). A loop over 

frequency is set up with Index k to compute the and W (hy 

ij^ IJk 

means of the subroutine WMAT) and then sums over frequency are taken within 
the loop to arrive at the coefficients in the final equations: in 

the bound case and In the continuum. 

Once the integration over frequency is completed, the program 
sets up a matrix, called EM, for solution of the simultaneous equations 
to get Sj or l/by. These are then punched in cards etlong with the pop- 
ulation ratios to be used as inputs to the next step In the Iteration. 

^IJ used to determine new line or continuum radiative 

rates which are then used as input to code 1. New values for the popula- 
tions and the parameters B®, and cj are obtained and are used in 
beginning the next step of the iteration. 
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DESCRIPTION OF SUBPROGRAMS 


SUBROUTINE WMAT (NZ, TOUT) 

WMAT calculates ^ those areas of the matrices 

where Is less than TOUT. 

FUNCTION NOSONI (A, X, L. LMAX) 

NOSONI is a matrix inversion routine using the method outlined on page 
4 34 of Hildebrand, Introduction to Numerical Analysis (New York, 19 56 ). 

A io a matrix of order L with column dimension IMAX. Its elements are 
assumed to be stored columnwise in the usual Fortran manner, X is work- 
ing storage of length L. The inverse of A will replace A. Upon return, 
NOSONI = 1 if inversion went properly, = 0 if a divisor is zero, in which 
case, A may contain garbage. 

FUNCTION ESB (X, N) 

This routine uses the function EXIN to obtain Ej^(X), the exponential integral 
of order 1 in X, It then calculates Ej^(X^, where N = 2, 3 or 4' from the rec 
recurrence relation i [e"^- ZEj^(x)] 

FUNCTION EXIN (Y) 

liDCIN obtains E^(y) from polynomial approximations shown in sections 5.1.53 
and 5- 1-56 of the Handbook of Mathematical Functions , National Bureau of 
Standards , 19^4 . 
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Introduction 


This Final Report summarizes the research .'erforined undnr 
NASA grant NCR 21-001-024. from 1 Oct.. 1972 thru 3i March. .974. 
This grant allowed the continuation of solar radio observat lon-^ 
which had previously been funded with "In-house" research funds. 

The research was performed under the direction of B. L. Go two , 
Antenna control, circuit construction, and routine cataloging of 
data were performed by R. J. Sneeringer. 


ObaervatLons and Data Reduction 

ObaervntJons commenced In November, 1972 (after a month 
of receiver Improvement), and continued until the end of the last 
manned Skylab mission in February, 1974. As explained in our 
original proposal, it was not possible to obtain unlnterruped 
sunrise to sunset coverage, due to the sharing of the 60 ft. antenna 
with other laboratory programs. Despite this limitation, 1408 hours 
of observations were obtained. 

The spectra were recorded in real time, both on film and 
magnetic tape. The filming was performed with a co.-.tlnuous motion 
camera running at 0.8 in/min. This rate of film travel w'll only 
allow 0.2 H to be re.solved, but this Is sufficient .c inu-.vals 

which require further study with the full time resolution of 0.01s 
that Is preserved on the magnetic tapes. High speed replays of all 
bursts were subsequently made at times which caused the least 
Interruption to the observing program, l.e., during periods of 
exceptionally low solar activity, during the unmanned portions of 
the Skylab flight, etc. 

A catalog of the observations is given in Appendix I. A 
similar catalog covering the period May, 1973 - February, 1974, 
has been submitted to World Data Center - A (NOAA) , for inclusion 
in a catalog of ground based observations in support of Skylab. 
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Equipment Modification 


Our CRT display was modified so as to greatly improve the 
contrast obtainable on weak solar bursts, at the expense of decreasing the 
dynamic range displayed to 10 dB. The display was run in this high 
contrast mode throughout the period November, 1972 - February, 1974, 

For non-saturated filming of very strong bursts, the CRT controls are 
easily readjusted and the magnetic tape replay through the spec- 
trograph display. 

An automatic intensity calibration scheme (hourly) was 
constructed and Installed In November, 1973, Up until this 
time Intensity calibrations were performed manually once or 
twice a day. 


Research Results 


Preliminary reduction of our high time resolution observations 
has revealed the fact that there is often considerable curvature 
present at the low frequency extremity of the fast-drift bursts in 
our frequency range (see Figure 1). The entire burst occurs on a 
time scale of 0.5 s, so the detection of this curvature would have 
been impossible with the 0.2 s or greater scanning period of former 
studies. Stimulated by this new finding we are currently considering 
the following hypothesis: The majority of fast-drift declmetric 

wavelength bursts are the result of streams of electrons that are 
guided along closed magnetic field lines. This hypothesis appears 
to be capable of explaining the surprisingly loose correlation 
between type III bursts at decimeter and meter wavelengths. It can 
also account for the observation^ that shorter lived decimeter 
wavelength type III bursts exhibit higher drift rates. 

On the theoretical side, we have studied pulsating type IV 
2 

solar radio bursts. The most interesting result of this study 

i. Young, C. W. , Spacer, C. L. Moreton, G. E. , and Roberts, 

J. A.: 1961, Astrophys. J. . 133 . 243. 

r. G'liwols, B. L. : 1973, Solar Phys » , 33 , ^75 ireprintnd In 

Appendix II). 
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was the way in which the Razln effect enhances the depth of modulation 
of the pulsations. Also of interest Is the fact that when significant 
synchrotron self-absorption is present, the pulsations break up into 
two distinct bands which pulsate 180* out of phase with each other. 

Colloboration with Other Groups 

Our Data has been supplied directly to two of the investigators 
on Skylab. These are; Dr. E. B. Mayfield of the Aerospace Corporation 
(experiment S-056), and Dr. A. S. Krieger of AS & E (experiment S-54). 
Data has also been furnished to Dr. J. C. Brown who is currently visiting 
the Astronomical Institute at Utrecht. 

A catalog of our observations during the first manned Skylab fllglit 
was included as part of a NOAA publication^. A complete catalog 
covering the entire Skylab flight has recently been submitted through 
the same channels. 

Scientific Meetings and Publications 

Gotwols, B. L. : "Pulsating Type IV Solar Radio Bursts", 

presented at the 140th meeting of the AAS, June 27, 1973, 

Columbus, Ohio. Abstract - BAAS . _5, 340, 1973 

Gotwols, B. L. ; "Solar Radio Pulsations", presented at 

lAU Symposium no. 57 on The Solar Corona, September 13, 1973, 

Surfers Paradise, Australia. Abstract - to be published in 
the Proceedings. 

Gotwols, B, L. : 1973, Solar Phys . 33 , 475 (reprinted in Appendix II). 


3. Coffey, H. : 1973, "Preliminary Catalog of Ground-Based Skylab- 

Coordinated Solar Observing Programs for the Period May 28 to 
July 26, 1973", World Data Center - A for Solar - Terrestrial 
Physics, NOAA, Boulder, Colorado. 
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z 
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Fig. i Solar bursts recorded on 21 Nov. 197^i: (a) high time resolution dynamic spectra of fast drift 

bursts, many of which exhibit a low frequency turnover. Much detail has been lost in the process 
of reproduction, (b) schematic tracing of the original film from which Figure 2a was produced. 
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APPENDIX I 


CATALOG OF SOLAR RADIO OBSERVATIONS 
TAKEN AT THE APPLIED PHYSICS LABORATORY 
November, 1972 - February, 1974 


ORIGINAL 
OF POOR QUALITY 
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THE M0PKIN« UNiVCRSITY 

applied physics laboratory 


A Guide for Classification of Solar Radio Bursts 
Observed with the APL Spectrograph^ (565-1000 MHz) 


TABLE I - Description of the various types of spectra 


SPECTRAL 

TYPE 

APL 

SYMBOL 

DESCRIPTION AND COMMENTS 1 

1 

T- 

1 

Storm bursts 


2 

Slow drift bursts 

i III 

3 

Fast drift bursts; v> 100 MHz/sec 

1 

! IV 

4 

1 ^ 

Prolonged continuum j 

V 

i 5 ' 

1 

1 Brief continuum (normally following type III bursts)! 

- 

* 6 

Intermediate drift bursts; v 30-100 MHz/sec * 

UNCLF 

1 UNCLE j 

1 

Unclassified activity. 


TABLE II 

- Symbols appended to the spectral type 

SYMBOL 

DESCRIPTION 

P 

Pulsations 


Small group (< 10) of bursts 

j GG 

1 Large group 10) of bursts 

I c 

j Underlying continuum 

i u 

' U-shaped burst of type III j 

i RS 

t 

' Reverse-slope burst ! 

1 DP 
1 

Drifting pair j 

1 



Intermittent activity in this period 1 


TABLE 

III - Intensity Scale 



1 

SYMBOL i 

FLUX DENSITY XI o"^^ Wm"^ Hz"^ 


r 

i 

1 

25 - 65 


i 

j 

2 

65 - 650 


i 

1 

3 

> 650 


1 


1 Gotwols, B. L. and Phipps J. , 1972, Solar Phys. 26, 386. 
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JHU/APU 










SOLAR 

RADIO 

EMISSION 




spectral observations (565-1000 MHZ) 






NOV. 

1972 


times 

OF 


HURSTS 



DBS. 

(UT) 





date 

start 

emd 

ST ART (UT) 

FN('(UT) 

INT 

. TYPE remarks 

IS 

2039 

2054 





16 

1329 

1 447 






1510 

2052 





17 

132fl 

1 725 






1727 

2M6 






1520 

1932 

1933.2 

1933.4 

1 

3(j 


1836 

2 '53 





21 

13A0 

2053 

1347.8 

1348.1 

2 

36 




1 349,5 

1349,7 

1 

36 




1351.0 

1351,1 

1 

36 




1432.9 

1435,5 

3 

366 




1519,0 

1519.4 

2 

366 




1525,2 

1525.3 

1 

3 




1529,0 

1536.4 

2 

366 




1615.0 

1616.2 

3 

366 




1617.5 

1619,7 

3 

366 




1640.4 

1640.9 

2 

36 




1642.2 

1 642.4 

2 

36 




1728.0 

1728.1 

1 

36 




1730,1 

1732.7 

3 

366 




1 756.3 

1 758.2 

3 

366 




IHl 7.5 

1817,6 

1 

36 




1833.2 

1833.4 

2 

36 




1850,1 

1850.2 

2 

36 




1902.8 

1903,9 

3 

366 




1905.2 

1905,8 

2 

366 

2? 

13A2 

1 ^55 

1422-6 

1423,0 

1 

36 


2016 

2051 

1631 .1 

1631 ,6 

1 

UNCLE 

27 

1311 

1456 






1525 

20 50 





?M 

1314 

2^13 





2P 

1305 

1‘='39 






1623 

2o51 
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JHU/APL 


SOLAR RADIO EMISSION 
SPECTRAL OBSERVATIONS (565-lOrtn MHZ) 


NO'/, 1972 


TIMES OF •BURSTS 

OHS. (UT) 


date 

start 

E^!D 

STart(UT) 

ENO(UT) 

INT, 

TYPE 

3/^ 

1301 

2046 

1522.0 

1523.6 

2 

UNCLE 




1711.7 

1711.9 

1 

3G 




171S.5 

1715.5 

1 

36 




1R16.4 

1917.3 

2 

366 


remarks 
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JHU/APL 


Solar radio emission 
spectral observations <565-1000 MHZ) 


DEC* 1972 



times 

OBS, 

OF 

<UT> 


BURSTS 


DATE 

START 

End 

START (UT) 

ENO(UT) 

INT. type 

ns 

1259 

1431 

I35ti 

1606 




n« 

1953 

20 49 




1 1 

1304 

1723 

1634 

1 q 1 0 




1? 

130R 

1950 

1914 

2006 

1615,0 

1742.2 

1615.2 
1 742.6 

1 36 

3 30RS 

n 

1327 

1459 

1659 

2010 

1456 
1*^16 
190 0 
21 13 




14 

1257 

1901 

1R57 

2033 




1 s 

1301 
150 3 
IbO? 

1427 

1S39 

2009 




1 s 

1257 

151« 

1 447 

2 26 




19 

1257 

1414 

14 1(. 

1443 




2 1 

1331 

1711 

1956 

2104 

1617 
1924 
2040 
21 14 




21 

1257 

1530 

1924 

1 S25 
1956 
2050 




22 

1304 

1P55 

1321 .6 

1321.8 

1 3G 

2^ 

1301 

1529 

IsOO 

1945 

2015. 1 

2016.1 

3 36 


remarks 


LOW FREQ FND 
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JHU/APL 


SOLAR RADIO EMISSION 
SPECTRAL observations (565-1000 MHZ) 

DEC. 1972 



TIMES 

OF 


OHS, 

<UT) 

Date 

START 

Eno 

26 

2015 

2054 

27 

1301 

1^07 


1703 

1915 


2013 

2042 

29 

1307 

1522 


1643 

1«19 


1845 

2''48 


BURSTS 


FNoUJT) INT. TYPE 


REMARKS 
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JHU/APU 


SOLAR RADIO EMISSION 
SPECTRAL observations <565-1000 MHZ) 


JAN. 1973 


times of hursts 

08S, (UT) 

DATE START EMO START(UT) ENi^(UT) INT. TYPE REMARKS 

04 1306 le;O0 

1519 1H35 

1906 2051 

05 1311 1*34 195S.8 19S7.2 3 3G 

1920 2058 

OH 1304 1617 

1625 1634 

no 1537 1607 

17 1323 1546 

1559 1710 

2016 2n43 

IH 1321 IPIO 
1925 2055 

19 1307 1735 

1819 1947 

2? 1414 1446 

1454 1608 

1645 1752 

1826 1933 

2006 2126 

23 1308 1725 

1850 2050 

24 1307 1656 

1725 2003 

2034 2n54 

25 1309 1937 

26 1307 1438 

1613 1732 

1806 2000 
2032 2n51 


297 



JHU/APL 


SOLAR RADIO EMISSION 
spectral observations (S65-1000 MHZ) 


JAN. 1973 


times of 'URSTS 

OHS, (UT) 

date start End start(ut> eni*<uti int. type 

29 1306 1745 

1817 1937 

30 1339 1 355 

1949 2000 

31 1337 1646 

1715 1732 


REMARKS 
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JHU/APL 


SOLAR RADIO EMISSION 
spectral observations (565-1000 MHZ) 


PER. 1973 



times 

OBS. 

or 

(UT) 


HURSTS 



date 

START 

End 

STaRT(UT) 

ENn (UT) 

INT. 

type 

01 

1304 

1546 

1646 

1813 

1508 
1616 
1734 
2 '38 

2017.6 

2017.7 

2 

3U 

0? 

1305 
1618 
1 758 
1914 
1950 

1514 

1727 

1834 

1933 

2053 





ns 

131 3 
1627 
1843 

1 538 

1 7 37 

2 <53 





06 

1305 

1735 

1945 

1702 

1912 

1953 





09 

1638 

1943 

1713 

2050 





1? 

1424 

1615 

1703 

1549 

1640 

2136 

2038.3 

2038.3 

1 

3 

1 3 

1306 

1544 

1S17 

2nS0 





1 A 

1257 

1656 

1 627 

2 56 





IS ’ 

1302 

1626 

1^15 
2 54 





} 6 

1422 

1557 

1 527 
2106 





22 

1511 

1656 

1847 

1 ^46 
1018 
2t53 






remarks 



JHU/APU 


OATE 

P3 

2ft 

2 « 


SOLAR RADIO EMISSION 
spectral observations (565-1000 MHZ) 

FEB. 1973 


times 

OMS. 

OP 

(UT) 


>URSTS 



STAPT 

End 

start (UT) 

FN. (UT) 

INT, 

TYPE 

1321 

H45 

1546 

1921 

2010 

1410 

1517 

1847 

1943 

2055 





1307 

1603 

173? 

1528 

1702 

2o54 

?01?.l 

2013.3 

1 

UNCLF 

1305 

1459 

1650 

1839 

2041 

1427 
lft07 
1805 
2 0 35 
2050 

1733.0 

1735.4 

1733,0 

1735,4 

2 

2 

unclf 

30 


remarks 
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JHU/APL 


SOLAR RADIO EMISSION 
SPEOTRAL OBSERVATIONS (S6S-1000 MHZ) 


MAR. 1973 



times 

OHS, 

OF 

(UT) 


HURSTS 



date 

start 

Ek'O 

STaRT(UT) 

ENi. (UT) 

INT. 

TYPE 

01 

1324 

1429 

1602 

1708 

1359 

1501 

1635 

1849 


. 



0? 

1253 

1633 

1432 

1743 

• 




0*^ 

1314 

1S52 

1711 

150B 

1650 

2n56 





OR 

1301 

1419 

1611 

1724 

1348 
1532 
1 647 
2'^47 





09 

1518 

1636 

1603 

1805 





I't 

1553 

1651 

1644 

2030 





11 

1306 

1351 





17 

1308 

2009 





IH 

131? 

2 23 





IR 

130 1 
1346 
162] 

1321 
1459 
2 44 

18?0.0 

1820.0 

1 

3G 

?? 

I3iil 

1356 

1 t27 

2 51 





23 

1326 

1504 

1659 

1 4411 
1634 
2041 

1 404.8 
1919,2 
2040.8 

1404,8 

1919.4 

2041.2 

3 

1 

3 

3'^ 

uncle 

3GG 

24 

1304 

162« 

1314.1 

1314,5 

2 

3G 


remarks 



JHU/APL 


Solar radio emission 
SPECTRAL OBSERVATIONS <565-1000 MHZ) 


MAR. 1973 



times 

OhS, 

OF 

(UT) 


BURSTS 




date 

start 

End 

start <UT) 

EN'MUT) 

INT. 

TYPE 

remarks 

?5 

1310 

2013 

1616. B 

1610.2 

1 

3GG 


30 

1338 

152A 

1453 

2o42 

143R.2 

1526.0 
1613.6 

1717.0 

1 439.4 

1530.1 
16)4,2 

1718.2 

1 

3 

2 

1 

UNCLE 

3G6 

3G6 

UNCLE 



302 



JHU/APL 


SOLAR RADIO EMISSION 
spectral observations (565-1000 MHZ) 


APR* 1973 


times Of worsts 

OBS, (UT) 


date start 

End 

start (UT ) 

ENMUT) 

INT. 

TYPE 

REMARKS 

01 1740 

IQOO 

1 747.3 

1 7^f7,8 

3 

30RS 




1R4?. 1 

1 fl42.3 

2 

3<ilJ . 
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JHU/APL 


SOLAR RAOIO EMISSION 
spectral observations (565-1000 MHZ) 


may 1973 



times 

OHS, 

OF 

(UT) 


‘BURSTS 



date 

start 

EmD 

start (UT) 

ENn (UT) 

INT. 

TYPE 

17 

1749 

2a26 

1749, or 

1025.3 
1027.9 

1915.3 

1750.5 
1027.1 

1828.5 
1917.9 

2 

2 

1 

3 

3GG 

3GG 

3GG 

3GG 

IH 

1203 

1958 

1527.4 

1531 .6 

2 

4P 

19 

1210 

1026 





20 

1211 

1936 

1654.3 

1 654,4 

1 

3(i 

21 

1211 

1714 

1754 

1055 

1653 

1749 

1052 

1955 





24 

1345 

1424 

1601 

1415 

1550 

1957 





2S 

1201 

1639 

1823 

1635 

1H16 

1959 





?.h 

1155 

1322 

1639 

1725 

1310 

1619 

1718 

1927 





27 

1205 

1550 

1543 

1031 

133) .2 
1606.2 

1331.6 

1607.0 

1 

2 

3G 

3GG 

?R 

1224 

1029 





29 

1630 

1956 





3) 

1206 

1957 






REMARKS 
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JHU/APL 


Solar radio emission 
SPETTRAL observations (56S-iOOO MHZ) 

JUN. 1973 


times of BURSTS 

O0S. (UT) 


DATE 

START 

Emo 

start (UT ) 

ENf^ (UT) 

INT, 

TYPE 

remarks 

01 

1157 

1955 






03 

1231 

1906 






04 

1234 

1459 

1822 

1453 

1751 

2'il9 






09 

1231 

1926 






1 0 

115ft 

1918 






1 1 

1219 

1955 







1245 

1523 

2004 

1 S 16 

2000 

2047 






13 

1201 

1950 

1322.0 

1322,0 

2 

UNCLE 

NARROW BANDWIDTH 

14 

1201 

1503 

1432 

1953 






IS 

120 B 

1954 

1409,3 

14 09,4 

2 

3G 


1 ft 

1224 

1539 

1S35 

1926 

1422.0 

1427.9 

1426.0 
1 428.4 

.3 

2 

4*-^ 

3G 


1 7 

1203 

1926 






IB 

1220 

1510 

1440 

1952 






21 

1226 

1510 

1507 

1956 






2 ? 

1159 

1930 







23 

1238 

1B52 

25 

1214 

1356 


14 33 

1 746 
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JHU/APL 


SOLAR RADIO EMISSION 
spectral observations (S65-1000 MH7) 


JUN, 1973 


times of ^URSTS 

OMS, (UT) 


oatE 

start 

En') 

start (UT) 

EN'i(UT) 

INT, 

TYPE 


1850 

1020 

1 R5P*9 

1059.7 

2 

3GG 

20 

1202 

I'^IO 

131 0,0 

1311.5 

2 

3G 


1855 

1031 

1515.4 

1517.4 

2 

4P 




loop, 6 

19 )9,5 

1 

4 

30 

12U 

1018 

1516,3 

1518.0 

3 

366 


remarks 
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JHU/APL 


SoLftR RADIO EMISSION 
SPECTRAL OBSERVATIONS (565-1000 MHZ) 


JUL, 1973 



times 

or 



OBS. 

(UT) 


DATE 

START 

End 

start (UT) 

01 

1210 

1925 


D? 

123H 

l^^SO 



1500 

1953 


0 5 

1155 

1920 


0 7 

1229 

1449 



1456 

1930 


OR 

1217 

1930 


09 

1211 

1 752 

1650,3 


1820 

1938 


12 

1154 

1251 



1307 

1432 



1502 

1H05 



1832 

1916 



1938 

1954 


1 3 

1156 

1918 


1 4 

1139 

1845 


15 

1131 

1 B49 


15 

1243 

1334 



1357 

1 549 



180 9 

1947 


21 

1129 

1725 


22 

1126 

I?37 



1250 

1Q55 


23 

1226 

1314 



135? 

1 944 



1547 

1731 



1825 

2000 



•BURSTS 


EN(KUT) INT. TYPE 


1650.3 2 311 


25 


REMARKS 


1209 


145B 
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JHU/APL 


SOLAR radio EMISSION! 

SPECTRAL OBSERVATIONS (S65-1000 MHZ) 



times 

OF 



OMS, 

(un 


iATE 

start 

End 

STflRT(UT) e 

29 

1120 

1P56 

1309 1 

30 

1219 

UOO 



1650 

1726 



1751 

1B29 



1901 

1055 



JUL. 1973 
«UR5TS 

Ni.(UT) INT. TYPE REMARKS 

3S3 2 A 
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JHU/apl 


SOLAR RADIO EMISSION 
SPECTRAL ORSFRVATIOmS (56S-1000 MH7) 

AUG. 1973 



times 

OF 


HURSTS 




OBS, 

(UT) 





hate 

start 

EmO 

START <UT1 

ENi>(UT) 

INT. 

TYPE 

01 

11S5 

1446 






1729 

1011 






1911 

1Q55 





f)? 

120 3 

1 329 






15S4 

1732 






IHOl . 

1955 





0 3 

1244 

1423 





0 4 

1130 

1«42 





OS 

1132 

1416 






1423 

1«45 





OFi 

1214 

1712 






1 740 

1920 





oo 

1 isi 

1450 

1 151 .OF 

1151.0 

3 

3GII 


1456 

1727 

1153.4 

1 154.1 

1 

UNCLF 


1R44 

1945 

1414.9 

1415.9 

2 

3GU 




1550.4 

1552.6 

2 

36GU 

n 

1156 

1657 






1723 

1748 






1005 

1026 






1055 

1 943 





1 1 

1134 

1 045 





1? 

1 1 37 

1 047 





1 3 

1210 

1 4 06 






141? 

1 700 






100 3 

1941 





\f^ 

1203 

1>^3J 






1603 

1952 





1 7 

11S7 

1429 






1436 

1 643 






1820 

1950 






REMARKS 
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JHU/APL 


SOLAR RADIO EMiSSIOfg 
spectral OBSERVATIOMS (565-lOOn MH7> 

AUG, 1973 


times Of 

OOS. (UT> 


Hate 

start 

EmO 

1 A 

1 149 

1045 

IR 

1151 

1254 

1?44 
1 844 

20 

1217 

1318 

1637 

l?47 

1449 

1925 

23 

1 156 
1450 
1550 
1842 

1418 

1524 

1012 

1954 

24 

1201 

1519 

1454 

1535 

30 

1213 

1055 

31 

1204 

1225 

1219 

1941 


-URSTS 


STart<UT) ENi(UT) INT, TYPE 


1503.5 15<.3,B 1 3G 


remarks 
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JHU/APL 


SOLAR RADIO EMISSION 
spectral OBSERVATIOMS (S65-1000 MH7) 


SEP. 1973 


times of BURSTS 

OHS, (UT) 


date 

start 

Emu 

ni 

1204 

1 A26 


1747 

1925 

0 2 

1209 

1?19 


1422 

1S35 

06 

1215 

1226 


1346 

1359 


1A22 

1458 


1622 

1647 


1 746 

I P44 


1930 

1946 

0 7 

1211 

1320 


1519 

1937 

OR 

1217 

1918 

1 0 

1225 

1443 


1510 

1 748 

13 

1215 

1338 


1341 

1457 


1615 

101 0 


IBlfl 

1944 

14 

1208 

1423 


1 454 

1938 

IS 

1206 

1046 

16 

1512 

1 727 

17 

1217 

1 422 


150? 

1(442 


1657 

1 950 

20 

1522 

16 38 

21 

1206 

14 13 


1442 

I'^IU 


1601 

2nl7 


STaRT<UT) EN.i(UT) 


IPIA.9 12P3.4 

ib;>3.3 

U3S.7 1635,7 


INT, TYPE REMARKS 


2 30G 

1 UNCLF 


2 3 
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JHU/APL 


Solar radio emission 

SPECTRAL OBSERVATIONS {565*1000 MHZ) 


SEP, 1973 


times of ^UmSTS 

ORS. (UT) 


DATE 

start 

EmD 

start (UT) 

EM ) (tj D 

INT. 

TYPE 

remarks 

22 

1216 

1A24 

1404 

1924 






23 

1220 

1019 

1749.9 

17SO.O 

2 

36 


2A 

1339 

1A51 

1730 

1425 

1708 

1R13 

1759.1 

1759. ? 

1 

36 





/ 
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JHU/APL 


Solar radio emission 
SPECTRAL observations (565-1000 MHZ) 


OCT. 1973 



times 

OBS. 

OF 

(UT) 


HURSTS 



date 

start 

E^|^ 

ST ART (UT ) 

FN( (UT) 

I NT. 

TYPE 

01 

1242 

1404 





OS 

150? 

1850 

1 R46 
1946 





12 

1224 

1944 

1 4S 0 . 8 

1450.9 

1 

3 

13 

1217 

1927 





1 4 

1216 

1925 





IS 

1631 

1829 

1 755 
2000 





IK 

1251 

1955 





19 

1219 

1405 

1 332 
1930 






1227 

1923 





?1 

1219 

1R30 





2 ? 

1243 

1 945 





2S 

1232 

1352 

154 0 

1318 
1S05 
1 948 

1 5S2.8 

1 553,0 

1 

3G 

26 

1224 

1450 

) 447 
1945 





27 

1224 

1 742 

1 SSO. 0 
ISSS.5 
1 5 '7 

1550.0 

1559. 1 
1613 

2 

2 

2 

UNCLF 

6GRS 

4 

?K 

1227 
1 84S 
1950 

1 039 

1 93H 
2056 






REMARKS 



JHU/APL 


SOLAR RADIO EMISSION 
SPECTRAL observations (565-lOOfl MHZ) 


NOv. 1973 


times of 
Ohs. (UT) 


Date 

start 

End 

If) 

1308 

lAOfl 

1511 

1342 

1428 

1600 

17 

1351 

2^23 

IH 

130 8 
1324 

1321 
1 336 

19 

1323 

1737 

1914 

1639 

l«3H 

2-26 

29 

1333 

lAOS 

1700 

2014 

1345 

1423 

1803 

2047 

2A 

1324 

1014 

1852 

1351 

1R47 

2020 

2S 

1303 

1340 

27 

1406 

1917 

1914 

2106 

2P 

1431 

2107 

29 

1308 

2031 

39 

1305 

1726 


-UPSTS 


START(UT) FN(MUT) int, type 


IT'sin.T 17S0.7 1 3 


remarks 


{ 


ORIGINAL PAGE IS 
OF PGCB QUALIFY 



JHU/APU 


SOLAR RADIO EMISSION 
SPETTRAL observations (S65-1000 MHZ) 


DATE 

01 

0? 

0 3 

04 

05 

n#s 

07 
OH 
OR 
10 
1 1 
1? 

1 3 

14 

IH 


DEC. 1973 


TIMES 

OHS. 

OF 

(IJT) 


hursts 




start 

Emu 

start (UT ) 

ENn (un 

INT. 

TTPE 

REMARKS 

1307 

1916 






131H 

2004 






1327 

163? 

1402 

2n40 






1319 

1623 

1431 

2050 

1349.0 
1949.6 

1950.0 

1 349.0 
1950.9 

1950.1 

2 

1 

1 

3 

3G 

3PS 



1323 1411 

H24 U12 

I6?n 1455 
1745 1941 

200fl 2040 

1325 1 M4 

1 7 1 fl 1930 
2016 2046 

1305 2020 

1313 2-15 

133? 2-14 

1325 2030 

1330 U32 

1512 1545 

193H 2110 

131M 1942 

1521 1907 

1910 2040 

1301 1400 

1603 1734 


125R 


1434 


31.5 



JHU/APL 


SOLAR RADIO EMISSION 
SPECTRAL observations <565-1000 MHZ) 


DEC. 1<)73 

timfs op -URSTS 

OHS, (UT) 

Date start End start (UT) enikuT) int. type remarks 

15 1820 2023 

16 1301 1411 

17 1619 IR23 

1927 ?020 

IR 1328 1949 

2020 2048 

1° 1322 1021 

20 1318 2030 

21 1314 1619 

1624 1935 

2? 1309 1922 

23 1303 1S03 

1542 1931 

24 1322 2110 1521.2 1521.3 3 36 

26 1327 2050 

27 1715 2057 

28 1311 1906 

30 1335 1516 

1553 2027 

31 1310 1540 

1545 1600 

1604 1959 


316 



JHU/APL 


SOLAR RADIO EMISSION 
SPEfTRAL OBSERVATIONS (565-1000 MMZ) 


JAN* 1974 


times of mursts 

OBS, (UT) 

Date start end start(ut) enc)(ut) int, type 


01 

1418 

2l?2 

0? 

1343 

2o52 

01 

132? 

1908 

04 

1316 

143? 

1620 

2009 

1422 

1605 

1«18 

2043 

0 5 

1318 

1945 

06 

1 330 

2027 

07 

1313 

1935 

oo 

1429 

1746 

10 

1504 

2031 

1957 

2046 

1 1 

1319 

1449 

1716 

ns4 

1712 

1P52 

1? 

13 04 
1826 

1922 

1947 

1 1 

1311 

1945 

14 

1312 

1819 

2007 

1726 
1 942 
2040 

15 

1307 

1759 

1733 
1 945 

1 7 

1305 

2100 

IH 

1 304 
1935 

1 321 
2 '-54 


REMARKS 


317 



JHU/APL 


Solar radio emission 
spectral observations (565-1000 MHZ> 

JAN, 197A 


times of hURSTS 

ORS, (UT) 


Date start Emd start (ut> eno(ut) int, type 

IR 1312 1RA9 

20 1320 1606 

1614 1929 

21 1307 n35 

13A1 U40 

1648 1748 

1829 1P41 

1859 1921 

22 1305 IPOO 

1835 1900 

1928 2048 

23 1344 1522 

1526 1702 

1712 1P32 

1929 1946 

1950 2055 

24 1310 n40 

1344 1425 

1505 1538 

1617 IPOO 
1829 2053 

25 1316 1411 

1417 1548 

1702 1729 

1758 2040 

26 1309 1937 

27 1305 1942 

28 1304 1738 

1809 2029 j 

30 1417 1q 09 i 

1912 2105 


REMARKS 
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JHD/APL 


SOLAR RADIO EMISSION 
SPECTRAL OBSERVATIONS (S6S-1000 MHZ) 


JAN, 1974 



times 

OF 

HURSTS 



OBS. 

<UT) 



DATE 

START 

End start(UT) 

ENn(UT) INT. TYPE 

REMARKS 

31 

1819 

1956 
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JHU/APL 


date 

01 

0 ? 

03 

OA 

o«; 

Oft 

07 


SOLAR radio emission 
spectral observations (565-1000 MHZ) 

FEB. 197A 


times of bursts 

OflS, (UT) 

start EmO STaRT(UT) ENn(UT) INT. TYPE REMARKS 

1303 1709 

181? 2<nB 

I30fl 1F52 
1558 2033 

1324 2034 

1303 1733 

1801 2029 

1302 1411 

1428 lft5<» 

1729 1845 

1849 2039 

1332 1439 

1444 1628 

1650 2039 

1305 IftOO 
1628 2040 
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